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Abstract 

The numerical prediction, theoretical analysis, and experimental verification of the 
phenomenon of wave packet revivals in quantum systems has flourished over the 
last decade and a half. Quantum revivals are characterized by initially localized 
quantum states which have a short-term, quasi-classical time evolution, which then 
can spread significantly over several orbits, only to reform later in the form of a 
quantum revival in which the spreading reverses itself, the wave packet relocalizes, 
and the semi-classical periodicity is once again evident. Relocalization of the initial 
wave packet into a number of smaller copies of the initial packet ('minipackets' or 
'clones') is also possible, giving rise to fractional revivals. Systems exhibiting such 
behavior are a fundamental realization of time-dependent interference phenomena 
for bound states with quantized energies in quantum mechanics and are therefore 
of wide interest in the physics and chemistry communities. 

We review the theoretical machinery of quantum wave packet construction leading 
to the existence of revivals and fractional revivals, in systems with one (or more) 
quantum number(s), as well as discussing how information on the classical period 
and revival time is encoded in the energy eigenvalue spectrum. We discuss a number 
of one-dimensional model systems which exhibit revival behavior, including the 
infinite well, the quantum bouncer, and others, as well as several two-dimensional 
integrable quantum billiard systems. Finally, we briefly review the experimental 
evidence for wave packet revivals in atomic, molecular, and other systems, and 
related revival phenomena in condensed matter and optical systems. 



PACS numbers: 03.65. Ge, 03.65.Sq 
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I. INTRODUCTION 



The study of localized, time-dependent solutions to bound state problems in quantum 
mechanics has been of interest to both researchers and students of the subject alike since the 
earliest days of the development of the field. Schrodinger (fl*) and others (0) - (0) discussed the 
connections between the quantum and classical descriptions of nature by exhibiting explicit 
wave packet solutions to many familiar problems, including the cases of the free-particle, 
uniform acceleration (constant electric or gravitational field), harmonic oscillator (forerunner 
of coherent and squeezed states), and uniform magnetic field. Many such examples then 
appeared in early textbooks © - (0) only a decade later, discussing both wave packet 
spreading and periodic time- dependence in a way which is easily accessible to students even 
today. 

Despite Schrodinger 's hope (P) that "...wafe groups can he constructed which move round 
highly quantised Kepler ellipses and are the representations by wave mechanics of the hydro- 
gen electron..." (without spreading, as with the constant width harmonic oscillator packet 
he derived), early investigators soon found (0), (0) that such dispersion was a natural feature 
of wave packets for the Coulomb problem. 

Attempts at constructing localized semi-classical solutions (of the coherent-state type) for 
the Coulomb problem, following up on Schrodinger's suggestions, continued with theoretical 
results - (jM ) appearing in the literature much later. It was, however, the development 
of modern experimental techniques, involving the laser-induced excitation of atomic Rvd- 
berg wave packets, including the use of the 'pump-probe' (jT3) or 'phase-modulation' (l^ 
techniques to produce, and then monitor the subsequent time-development of, such states 
which led to widespread interest in the physics of wave packets. (For reviews of the subject, 
see fflgf ) - (|23).) Updated proposals for the production of such states in the context of Ryd- 
berg atoms, where one could study the connections between localized quantum mechanical 
solutions and semi-classical notions of particle trajectories, led first to the creation of such 
spatially localized states (|2^. to experiments which observed their return to near the atomic 
core (25) . and then the observation of the classical Kepler periodicity (|26[). (|27h of Rydberg 
wave packets, over only a few cycles in the early experiments. 

However, this interest also led to the prediction of qualitatively new features in the long- 
term time development in such bound state systems, such as quantum wave packet revivals. 
Parker and Stroud (|28[ ) were the first to find evidence for this behavior in numerical studies 
of Rydberg atoms, while Yeazell and Stroud (jHO), pll ) and others ((12), ()33l ) soon confirmed 
their predictions experimentally. 

The phenomenon of wave packet revivals, which has now been observed in many experi- 
mental situations, arises when a well-localized wave packet is produced and initially exhibits 
a short-term time evolution with almost classical periodicity (Td) and then spreads signifi- 
cantly after a number of orbits, entering a so-called collapsed phase where the probability is 
spread (not uniformly) around the classical trajectory. On a much longer time scale after the 
initial excitation, however, called the revival time (with Trev » Td), the packet relocalizes, 
in the form of a quantum revival, in which the spreading reverses itself and the classical 
periodicity is once again apparent. Even more interestingly, many experiments have since 
observed additional temporal structures, with smaller periodicities (fractions of Td), found 
at times equal to rational fractions of the revival time [pTrev/q)- These have been elegantly 
interpreted pH) as the temporary formation of a number of 'mini-packets' or 'clones' of the 
original packet, often with 1/g of the total probability, exhibiting local periodicities Td/q, 
and have come to be known as fractional revivals. Observations of fractional revivals have 
been made in a number of atomic (|3ll l - (|33h and molecular (l35l ) systems. 

A simple picture (l3fih of the time-dependence of the quantum state leading to these 
behaviors, modeling the individual energy eigenstates and their exponential {exp{—iEnt / h)) 
time-dependent factors as an ensemble of runners or race-cars on a circular track, is often 
cited. The quantum mechanical spreading arises from the differences in speed, while the 
classical periodicity of the system is observable over a number of revolutions (or laps.) For 
longer times, however, the runners/race- cars spread out and no correlations (or dumpings) 
are obvious, while after the fastest participants have lapped their slower competitors (once 
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or many times), obvious patterns can return, including smaller 'packs' of racers, clumped 
together, which model fractional revivals. 

A different metaphor involves the (deterministic) shuffling of an initially highly ordered 
deck of playing cards. One shuffling method involves splitting the deck into two equal 
halves, and then alternately placing the bottom card from each half into a pile, reforming 
and reordering the deck. After only a few such shuffles, the original order is seemingly 
completely lost and the cards appear to have randomized. After only a few more turns, 
however, clear patterns of ordered subsets of suits and ranks appear, increasingly so until 
after only eight such shuffles the deck has returned to its original highly ordered state. Only 
the simplest of mathematical concepts is required to describes these analogs of fractional 
and full revivals. 

The Dynamics of wave packets of highly- excited states of atoms and molecules, including 
a discussion of wave packet revivals and fractional revivals, and descriptions of the exper- 
imental observations of these phenomena, were discussed in the excellent 1991 review by 
Averbukh and Perelman (j37f ). while a nicely accessible general discussion by Bluhm and 
Kostelecky (38) has also appeared. There have been developments in the field since then, 
and many of the basic quantum mechanical concepts behind revival behavior have also be- 
gun to appear in the pedagogical literature, so it seems appropriate to provide a review of 
some of the fundamental ideas behind the short- and long-term behavior of quantum wave 
packets, describing both the classical periodicity and revival behavior of wave packets in 
many model systems, and their experimental realizations. 

The theoretical machinery required to understand many aspects of revival behavior is 
sufficiently accessible, and potentially of enough general interest, that our review of the 
subject will contain many tutorial aspects, such as 

(i) the use of familiar model systems (such as the infinite well and others) as illustrative 
of the fundamental concepts, 

(ii) an emphasis on examples in which exact revival behavior is found (to be used as 
benchmarks for more realistic systems), 

(iii) a focus on semi-classical methods, including the WKB approximation, which are ap- 
propriate for many wave packet systems constructed from large n energy eigenstates, 

(iv) general discussions of how information on both the classical periodicity and quantum 
revival times are encoded in the energy eigenvalue spectrum, 

(v) and references not only to the original research literature, but to many of the peda- 
gogical papers appearing on the subject. 

In these areas, and others, we hope to extend the reviews in Refs. (fST*) and (|38h in useful 
ways. 

We start in SecOwith a general introduction to the theoretical tools required for under- 
standing revival behavior, including t he a utocorrelation function. We then turn to discus- 
sions of many model systems (in Sec. IIII|) which illustrate various aspects of the quantum 
mechanical time development of localized wave packets, including background material on 
unbound systems (free-particle and uniform acceleration cases) and for ones exhibiting only 
periodic behavior (the harmonic oscillator). The infinite well is then discussed in great de- 
tail, as are other familiar one-dimensional problems. Two-dimensional quantum systems, 
especially quantum billiard geometries, are studied in Sec. IIVI We then briefiy review 
experimental evidence for revival behavior in Sec. |V] in atomic (Coulomb) and molecular 
(vibrational, rotational) systems, as well as in situations where the quantum revivals are 
due to the quantized nature of the electromagnetic field in two-state atom-field systems, 
or in the excitation spectrum of Bose-Einstein condensates, and finally we discuss related 
revival phenomena in a variety of optical systems. 
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II. GENERAL ANALYSIS 



A. Autocorrelation function 



The study of the time-development of wave packet solutions of the Schrodinger equation 
has a long history and often makes use of the concept of the overlap {{iptlipo)) of the time- 
dependent quantum state (lipt)) with the initial state (iV'o))- For example, early work by 
Mandelstam and Tamm f(39^ on the time-development of isolated quantum systems led to 
the inequahty 

m\i^o)\' > cos^ (^^^ , valid for 0<t<^, (1) 

where AH = ^J {H"^) — {H)"^ is the uncertainty in the free-particle energy of the wave packet. 
These ideas were used to study energy-time uncertainty relations, as well as the minimum 
time required to reach an orthogonal quantum state. 

This overlap is most often referred to as the autocorrelation function and can be evaluated 
in either position- or momentum-space to give 

/+00 r+oo 
rix,t)i;ix,t)dx= / (2) 
oo J — oo 

which emphasizes that in order for \A{t)\ to be large, the wave function at later times must 
have significant overlap with the initial state in both x- and p-space. 

For one-dimensional bound state systems, where a wave packet is expanded in terms of 
energy eigenfunctions, Un{x), with quantized energy eigenvalues. En, in the form 



with 



n=0 

/ + 00 
[un{x)]* tlj{x,0)dx 
-oo 



oo 
n=0 



(3) 
(4) 

(5) 



the most useful form for the autocorrelation function ff40i) is 



and the evaluation of A{t) in this form for initially highly localized wave packets will be one 
of the main topics of this review. Besides being of obvious theoretical value in the analysis of 
time-dependent systems, more physically, the autocorrelation function is important because 
it is very directly related to the observable ionization signal (1^ , (|30[ ) in the pump-probe 
type experiments where such behavior is studied experimentally. 

Studies of the long-term time development of quantum states is also a topic with a long 
history, often being discussed from a more formal point of view. For example, Bocchieri 
and Loinger ()4ll ) produced statements of a ' quantum recurrence theorem^ which attempted 
to generalize Poincare's theorem to show that systems with discrete energy eigenvalues 
would eventually return arbitrarily closely to the initial state, in the sense that the norm 
iV't ~ ''/'oil could be made smaller than any arbitrarily small number. (Such ideas have also 
3een examined (l4^ from a pedagogical point of view.) Similar studies in the mathematical 
literature often deal with analyses of almost periodic functions (^S), P^- Such long-term 
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time behavior has also been compared (|45| ) to Fermi-Pasta-Ulam recurrences (l46l ) in the 
dynamical behavior of non-linearly coupled oscillators. Investigations of recurrence ( revival- 
like) phenomena in the context of non-stationary (time-dependent) Hamiltonians ()47l ) have 
resulted in similar theorems. 

In this context, Gutschick and Nieto (jlSl) noted that "...«/ one waits long enough any state 
which has significant overlap with a finite number of states will eventually return to "almost" 
its original shape... it amounts to having a totally dispersed wave packet suddenly regenerate 
itself after very long times. Peres (|48i) also pointed out that very long time scale recurrences 
of this type would occur in systems with hydrogenic energy levels. However, both analyses 
focused on time scales which are much longer than those ultimately observed in wave packet 
experiments. For the laser-induced wave packets produced in modern experiments, the 
structure of the expansion in eigenstates in Eqn. is such that a quite general systematic 
analysis of the time scales involved in the problem is possible and that is the topic of the 
next section. 



B. Time-dependence of one-dimensional localized wave packets 

In general, the time-dependence of an arbitrary time-dependent bound state wavefunc- 
tion, ip{x,t), with the expansion in eigenstates, Un{x), of the form in Eqn. Q can be quite 
complex. However, in many experimental realizations, a localized wave packet is excited 
with an energy spectrum which is tightly spread around a large central value of the quan- 
tum number, uq, so that uq >> An » 1. In that case, it is appropriate to expand the 
individual energy eigenvalues, E{n) = En, about this value, giving 

E{n) ^ E{no) + E'{n,){n - n,) + ^^{n - n^f + ^^^(n - n,f + ■ ■ ■ (6) 

I D 

where E'{nQ) = {dEn/dn)n=no and so forth. This gives the time-dependence of each individ- 
ual quantum eigenstate through the factors 



E{no)t + {n- no)E\no)t + ^{n - nofE"{no)t 
+ ]:{n-nQfE"'{n^)t + -- 



= exp [—iuJot — 2'ni{n — Uoji/Td — 2Txi{n — noYt/Trev (7) 
-2ni{n - Uoft/Tsuper H ) 

where each term in the expansion (after the first) defines an important characteristic time 
scale, via 

''~Wi^\ ' \E"{n,)\/2 ' Wr-- \E"'{n^)\/%- ^ ^ 

The first [ujq = E{no)/h) term is an unimportant, n-independent overall phase, common 
to all terms in the expansion, and which therefore induces no interference between them; 
it is similar to the time-dependent phase for a single stationary state solution and has no 
observable effect in \ip{x,t) ^. 

The second term in the expansion is familiar from correspondence principle arguments 
(jioi ) as being associated with the classical period of motion in the bound state, so that 

^ 27r^ 
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This connection is perhaps most easily seen using a semi-classical argument and the WKB 
quantization condition, familiar variations on classical action-angle methods (|50h - (jsil). 

For a particle of fixed energy ii^ in a ID bound state potential, V{x), we have E = 
mv{xY/2 + V{x) and the short time, dt, required to traverse a distance dx can be obtained 
from this and integrated over the range defined by the classical turning points; this then 
gives half the classical period as 

^^^^ ^E-V{X) Ja ^^\I^Ja ^E-V{X)' ^ ' 

The WKB quantization condition in this same potential (with the same classical turning 
points, a, h) can be written in the form 

l-b 

^En - V{x) rfx = (n + Cl + Cn)T^h (11) 

J a 

in terms of the matching coefficients Cl, Cr. We recall that the appropriate values of these 
constants (jS^) are given by C^^/j = 1/4 at 'linear' or 'smooth' turning points, where functions 
are matched smoothly onto Airy function solutions (and effectively penetrates roughly 
1/8 of a wavelength into the classically disallowed region (jHI)), while Cl,r = 1/2 at 'infinite 
wall' type boundaries where the wavefunction must vanish. This expression can then be 
differentiated implicitly with respect to the quantum number, n, to obtain 

I \dEnldn\dx ^ , , 



]a 2^E^-V{x) 

This, in turn, can be related to the classical period in Eqn. ()1()|1 to give 
nn I 1 , l-ah 

Tci = T = v2m / — = dx = , ,^ , , , . (13) 

Ja - V{x) \dEjdn\ ^ ' 

The most obvious example of such a connection is for the harmonic oscillator, where the 
WKB condition gives the exact eigenvalues. En = {n + l/2)hu!, and the classical period from 
Eqn. Q is 

27rh 2Tih 271 
\dE/dn\ riw uj 

as expected. For the special case of the oscillator, all wave packets are exactly periodic and 
all higher order derivatives {E",E"', ...) vanish, so no other longer time scales are present. 
Nauenberg (jic? ) has also provided an elegant argument connecting E'{no), the classical pe- 
riodicity, and the structure of the autocorrelation function, which we review in AppendixIBl 
For future reference, even in the presence of higher-order time scales, we define the 
classical component of the wave packet to be 

CO 

t)^Yl an«n(a:)e-2'^^("""")^"*/^ ^ ^^^^{A^-^^'''"^- (15) 

n=0 fc 

where we define k = n — n^. This component can be used to describe the short-term 
(t ~ << T^rev) time-development and is especially helpful in discussing fractional revivals. 
We note that here, and elsewhere, we will henceforth ignore the overall exp(— iii^(no)t//i) 
phase factor. 
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The next term in the expansion in Eqn. (jSJ is given by 



T 

J. rt 



2Tih 

|^"K)|/2 



(16) 



which will be associated with the quantum revival time scale. This time scale determines 
the relative importance of the (n — term in the exponent for t > 0. It is responsible for 
the long-term {t » T^i) spreading of the wave packet in the same way that the difference 
in the p-dependence of the px and p^t/2m terms in exp{i{px — p'^t /2m) / h) in the plane wave 
expansion of free particle wave packets gives rise to dispersion; it can be shown quite directly 
(Ref. (j^Cj) and Appendix EI) that the spreading time is proportional to Trev for t << Tj-ev 
More interestingly, for times of the order of T^-ev, the additional exp(27rz(n — nQ)H/Trev) 
phase terms all return to unity, giving the t ~ time-dependence in Eqn. (II 5|) and a return 
to approximate semi-classical behavior. 

The anharmonic oscillator provides an example of the calculation of the revival time of 
relevance to physically realizable systems such as the vibrational motion of molecules. For 
example, the addition of a cubic term, Xhwx^, to the usual oscillator Hamiltonian yields (in 
second-order perturbation theory (jS^)) the energy eigenvalues 



E„ = hw 



{n + 1/2) 



15A- 
4 



n ■ 



1/2)^ - ^A^' 



hw + n{l — a) — an^~\ (17) 



where 7 = 1/2 — llA^/8 and a = 15A^/4. The effect of the anharmonicity (for either sign 
of A) is to increase the classical period {Td = 27r/(l — a)uj), but more importantly it also 
introduces a finite revival time given by 



27T 



T, 



d 



— ^ — » Ta 

ijja a 



for a « 1 . 



(18) 



As another example of the hierarchy of relative magnitudes of T^ev and which is 
possible, we consider a family of (symmetric) one-dimensional power-law potentials 
(Hi defined by 

. k 



V^k){x) = Vo 



(19) 



This family includes both the harmonic oscillator (when k = 2 and Vo/L"^ = muj'^/2) and 
the infinite well of width 2L (when k +00, for any Vq). The same WKB approximation 
considered above (see Appendix |X1 for details) shows that the energy eigenvalues scale (for 
large n) as 

E.,{k) oc n^^/^^+^^ (20) 

which does indeed give the correct large n behavior for the oscillator [E^ oc n) and the 
infinite well {En oc n^); it even gives the appropriate dependence on quantum number for 
the Coulomb problem {En oc for k = —1) so that these simple arguments are useful 
for the more realistic case of Rydberg atoms. (Many similar scaling laws can be derived for 
truly three-dimensional systems in different ways, as shown in Ref. (59).) 

For this general class of potentials, we find that the ratio of the revival time to the 
classical period is given by 



(fc) 



\E'{no)\ \ 
27vh J 



k + 2 



k~2 



no . 



(21) 



The revival time diverges for the special case of the oscillator {k = 2), as do all of the 
other higher-order time scales, while for the general {k 7^ 2) case one can quite generally 
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have Trev >> for no >> 1, which is typical for the locahzed wave packets studied 
experimentally. We note that for the experimentally important case of the Coulomb potential 

{k = —1), the hierarchy is given by Trev°^''^ = {2no/3)T^F°^''\ 

We note that several important model systems, including the infinite well, and the 2D 
{Em oc m^) and 3D (Ei oc l{l + 1)) free rotors (38) have energy spectra with no higher 
than quadratic dependence on the quantum number, implying that higher order derivatives 
vanish and longer time scales are sent to infinity. The interaction energy of t he p articles in 
Bose-Einstein condensates can be modeled in the form En = Uon{n — l)/2 (|60() providing 
another highly realizable physical example with a quadratic energy spectrum, and revival 
behavior. 

Continuing in this fashion, we note that the next term in the expansion is typically an 
even longer time scale, called the superrevival time, defined by 



2TTh 



super 



\E"'{no)\/Q 



with 



(k) 



super 
-L rev 



3{k + 2) 



no 



(22) 



(23) 



for the power-law family of potentials; we discuss superrevivals in more detail in Sec. III. El 
Since this type of behavior is generic to localized wave packets, we will first exemplify 
the short-term semi-classical periodicity, then the approach to a collapsed (incoherent sum) 
state, as well as the structure of quantum revivals and fractional revivals using a simple 
model energy eigenvalue spectrum and standard Gaussian distribution of eigenstates. We 
will assume that the expansion coefficients are given by 



1 



which gives the required normalization 



exp 



(n 



no 



\2-\ 



(24) 



E 



(25) 



to exponential accuracy for rip >> An >> 1. For the examples in the rest of this section 
we use no = 400 which is typically larger than that found in experimental realizations, 
but useful for making the visualized examples we present more transparent. We often use 
An = 6, but also will examine the effect of varying An as well; several notable experiments 
on Rydberg atoms (30), (33) have made use of wave packets constructed from 5 — 10 states, 
corresponding to An 2 — 5 in our notation. 

For the energy eigenvalues, we choose, for simplicity, a generalized anharmonic oscillator 
spectrum of the form 

En = 27r(n - an^/2 + l3n^/6) (26) 
and use two sets of parameters for comparison, namely 



Case 


A 


B 


a 


1/800 


1/800 







2 ■ 10-6 




2 


1.515 


T 

rev 


1600 


4444.4 


Trev /Tel 


800 


2933.3 
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Table I. 



For Case A, with only linear and quadratic terms, we expect exact revival behavior, and 
this case is arranged to be even more special by having Tj-ev/Td be an integer and so can 
be considered as an ideal case. The introduction of an additional small {(3 7^ 0) term in 
Case B will exemplify systems containing higher order terms, longer time scales, and hence 
only approximate revival behavior. 



C. Classical periodicity and approach to the collapsed state 

Even in the absence of a specific physical system, one can still examine, in great detail, 
the time-dependence of a generic wave packet using the autocorrelation function in the form 

00 

= ^ |aj V^"*/^ (27) 

n=0 

Nauenberg pffl, for example, has shown how to elucidate the approximate periodicity ap- 
parent in \ A{t}\ for bound state systems in a quite general way, and we discuss his approach 
in Appendix El but we first present some numerical examples. 

Using the values for Case A in Table I, we examine the behavior of A{t) over the first 
100 classical periods in Fig. ^ corresponding in this case to 1/8 of the entire revival time. 
In that plot (upper half), the initial classical periodicity is clear, as is the effect of the wave 
packet spreading, shown by the decreasing value of A{t) at integral multiples of the classical 
period (also seen in the bottom half), as well as the increasing 'width' of the (decreasingly 
small) peaks in A{t). In more specific model systems where we can examine the dynamics, 
we can elucidate the nature of the loss of coherence to the collapsed phase in more detail, 
as in Sec. I111.D.4I 

For times longer than approximately 50Tc/ (in this case), the \A{t)\ oscillates (more or 
less rapidly) about a constant value, indicated by the horizontal dashed line. To understand 
this, we examine the general structure of |^(t)P and note that 



\A{t)\- 



Y,\ant + 2Y,K?\arn?cos[^^^-^^] . (28) 



For time scales of the order of the classical period, the cos((-E„ — Em)t/K) terms are still 
highly correlated and reproduce the approximate classical periodicity. For longer times, the 
oscillatory components become increasingly out of phase and can lead to high frequency 
(/ >> liijTci) excursions around the constant value given by the first term in Eqn. 
namely 



n=0 

which we will label the incoherent limit of |A(t)p. For the Gaussian expansion coefficients 
in Eqn. ()24|) . this sum can also be done approximately (to the same accuracy as the nor- 
malization) giving 

00 

\Ainc\^ = J2 = . r- ~ 0.047 for An = 6 (30) 

which is included in Fig. ^ as the dotted horizontal lines and indicated by the horizontal 
arrows. The larger the number of states contained in the expansion, the smaller the resulting 
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incoherent value during the collapsed state; because of the discrete nature of the bound state 
spectrum, the autocorrelation function does not asymptotically approach zero as it would 
in the free particle case, for example, but oscillates about this generic plateau value. 

In the bottom half of Fig. ^ we also show the values of A{t) at multiples of the classical 
period (squares) and at times a half period away (t = {n + l/2)Tci, stars) showing how the 
'in-phase' ('out-of-phase') components first shrink (grow) until they are of the same order 
as the incoherent value of lAj^cP? and then clearly exhibit highly correlated behavior, which 
is the first hint of fractional revivals; some evidence of this is the appearance of highly 
oscillatory behavior at fractional multiples of Trev (note the 1/9 - 1/11 labels in Fig. ^a).) 

We note, for future reference, that the form of the probability density, |'?/'(x, t)p during 
the collapsed or incoherent phase of the time-development can be similarly written in the 
form 



[X. 



t)\^ = J2\an\^\Unix)\' + 2^ 



J2 <a^wn(x)w;„(x)e^(^»-^™)*/^ 



(31) 



The probability density will then oscillate around a 'static' value determined by the incoher- 
ent sum of the probability densities for each eigenstate, |Mn(x)p. In the large n limit which 
is applicable here, WKB methods can be used to approximate the energy eigenstates, and 
one can write 

Pjx) = Mx)\' ^ (32) 



En - Vix) 

where En,Tn are the quantized WKB energies and classical periods in Eqns. (fTT|) and (fT^ . 
Similar statements can be made about the momentum-space probability densities. 



D. Revivals and fractional revivals 



1. Revival time 



We now turn our attention to the longer term behavior of the autocorrelation function 
over an entire revival time and plot |A(i:)| over the interval (0, Trev) in Fig.|21for three cases. 
The top two plots, (a) and (b), correspond to the parameters of Case A in Table I, showing 
the effect of changing An; for larger (smaller) values of An, the spread in momentum values 
is larger (smaller), so that the quantum mechanical spreading time, to, is smaller (larger) 
and the rate at which the initial classical periodicity is lost is corresponding faster (slower). 
The initial size of the wave packet (and, as we will see, the subsequent size of the 'mini' 
wave packet components at fractional revivals) is smaller (larger) for bigger (smaller) values 
of An and this correlation is also obvious in that finer details are apparent for the larger 
An case. 

In the bottom plot in Fig. Efc), we show the results corresponding to case B in Table 
I, where one includes contributions from terms of order (n — ng)^. The overall pattern is 
similar, and many of the same individual distinct features are seen, but without the obvious 
symmetry or near exact revivals as in the case of purely quadratic n-dependence where no 
higher order time scales are present. The rich structure of features at rational fractions 
of Trev which is apparent in both cases, and which are obviously highly correlated with 
fractional values of \A{t = 0)p = 1, is seemingly a robust feature of all the examples, and 
we systematically explore their general structure next. 

To examine the generic behavior of the wave packet near the revival time, T^ev, we use 
only the first two derivative terms in the expansion in Eqn. (jHl) and write 

i;{x, t ^ Trev) = afcMfc(x)e-2"^*/^^'e-2^**'' where k = n - uq (33) 

k 
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and we focus on the effect of the 0{t/Trev) or k'^ terms on the exponential time development. 
The additional phases arising from such terms all give unity and the wave packet is said to 
have revived since 

^(x, t ^ Tre,) = J2 «fcUfc(x)e-2"^'=*/^^' = M^, t) (34) 

k 

has returned to something like its initial form, exhibiting the classical periodicity. In the 
special case when T^ev/Td is an integer (as in Case A in Table I, and the top two plots in 
Fig. 12)) the revival occurs exactly in phase with the original time-development, and is exact 
(in that \A{t)\ returns to exactly unity). 

The robust prediction that the wave packet exhibits approximately the classical periodic- 
ity (in phase or not) near t = T^ev is also apparent in Fig. |21 (where we now plot |A(t)P near 
several full or fractional revival times) even for admixtures of higher order terms, such as 
the Case B values in Table I where \ A{t ~ Trev) I < 1, but the return to the initial periodicity 
is apparent. 



2. Fractional revivals 

We next examine ip{x,t) near half a revival time, and find that 

^(x, t ^ Trej2) = J2 afcM,(x)e-2-'=*/^-e--^'=' (35) 

k 

with additional exp(— ttz/c^) = ±1 factors. We can compare this to the semi-classical time 
evolution by noting that 

Mx, t + = J2 a.«.(x)e-2-^(*+^-/^)/^- = ^ a,«,(x)e-2-^*/^- e"-^ (36) 

k k 

and the additional ± phase factors in Eqn. (|35|) can be written in the same form since both 
the 

(Eqn. (jSlD) e-"*'^ = (-1)'=' = (-1)^ = e"'"''' (Eqn. (37) 
terms give the same result for even and odd values of k. We thus see that 

^{x, t ^ Trej2) = ^,;(x, t + Td/2) (38) 

and the wavepacket also reforms near the half revival time, with the original classical peri- 
odicity, but half a period of phase with the initial wave form (at least for integral Trev/Td) 
and we see this behavior in Fig. El for both the A and B cases. The approximate revival at 
Trev/ 2 seen for Case B is, in fact, somewhat better than at the full revival time Trev since the 
higher order anharmonicities have had less time to effect the phase structure of the revivals. 
We then urn our attention to the quarter-revival time where we can write 

ij{x, t ^ Trev/^) = afe«fc(x)e-2"''=*/^-e--^'='/2 (39) 

k 

and the additional phase factor beyond the classical terms is given by 

k = 2l even g'^'^^OVs _ ^-2m ^ ^ ^^q^-^ 

k = 21 + 1 odd e-i^{2l+l?/2 ^ ^^2mfi ^-2ml g-iV2 = _i (40b) 
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or 



^-..fcV2 ^ ^ +1 ^ ^^^^ . (41) 
-i for odd 



This factor can be written in the forms 



e 



2 2 ^ ^ V2 

and, using the resuh in Eqn. we can see that 

^(x,t^T,,,/4) = Y,a,u,{x)e-'-'''/^^^e--'''/' 



k 



(43) 



^ [e-*'^/^ t) + e+*-/^ t + T,,/2)] 



Thus, the wave packet near Trev/^ (and 3Tret,/4 as well) consists of two, high-correlated 
copies of the original packet, in- and out-of-phase, each containing half of the probability. 
This structure is similar to that proposed in early suggestions made for Schrodinger cat type 
states (j6l| ). and experimental work on Rydberg atom realizations of this state (p^ ) has been 
discussed in just that context. 

The signal for such behavior is the presence of peaks in |^(t)P with half the magnitude, 
and half the classical periodicity, as each 'mini-packet' or 'clone' or 'fractional revival' state 
now approaches the same location in phase space as the initial wave form twice during each 
classical period; this behavior is indeed seen in Fig. IHl with the exact revival structure and 
phase relation to the initial state for Case A (with a maximum value of |v4(t)p = (Iv^)^ = 
1/2 as shown by the horizontal dashed line), and in a more approximate manner as for Case 
B. 



3. General structure of fractional revivals 

A very general pattern of well-defined revivals, characterized by temporally localized 
structures in A{t), of local periodicity T^/q and with magnitude |A(t)P = l/g, at various 
rational multiples of the revival time given by t = pTrev/<l (where j9, q are mutually prime) 
is obvious from Fig. |21 and leads to the term fractional revivals. Averbukh and Perelman 
(|3^ were the first to analyze in detail the mathematical structure of the additional phase 
factors arising from the exp{—2iTik'^t /Tre^) terms at such times to discuss the ^^Universality 
in the long-term evolution of quantum wave packets beyond the correspondence principle 
limif and we reproduce here, in part,^their elegant arguments, for completeness. (We note 
that it has been pointed out (f63f ). (f64 ) that this problem, especially the calculation of the 
autocorrelation function at fractional revival times, is similar to that of the evaluation of 
Gauss sums (65) which has a long history in the mathematical literature (|66[).) The case 
of odd q is most straightforward and we consider the more general arguments for even q in 
Appendix O 

We are interested in the structure of the additional phase terms at times t = pT^ev/q of 
the form 

-2mk'^t/Trev _ -2TTi{pk'^/q) _ -2iTi0k (aa\ 
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and especially how to write these phases with k"^ in the exponential in terms of similar 
factors, but with linear dependence on k, as in Eqns. ^7\\ and We first note that the 
phase factor in Eqn. (j44j) will be periodic in k, namely 



or Pe = ek = ek+i = ^{k + lf (modi). (45) 

q q 



For this to occur, one requires that 

2pk + pP 







(mod 11 



(46) 



which, for odd q, is satisfied if / = g. One can then expand any periodic function (periodic 
in k) in terms of the basis states 



and write 



^-2m0k _ ^-2mpk^/q 



S = 0,1,...,/- 1 



s=0 



-2nisk/l 



(47) 



(48) 



which can be described as a generalized trigonometric identity, finite Fourier series, or dis- 
crete Fourier transform. This is already a useful result as we can then write the wavefunction 
near a fractional {p/q) revival in the form 



tp{x,t pTrev/q) = y^afcMfc(x)e" 

k 
l-l 



■2TTikt/T^l^-TTipk^/q 



l-l 

i!ciix,t + sTa/l) 



(49) 



s=0 



This form implies that there will be as many as / 'clones' of the original wave packet, with 
amplitude bs and probability l&sp, differing in phase from the original packet by fractions 
(Tci/l) of the classical period. 

The expansion in Eqn. ()48|1 can be inverted by multiplying both sides by exp{27Tirk/l), 
summing over / possible states, and using the periodicity of the 0^ to obtain 



l-l l-l 

^2TTirk/l—2TTipk'^ /q 



E 

k=0 



E*- 



s=0 



E- 

.fc=0 



^2iTik{r—s) /I 



l-l 



[ISr 



(50) 



or 



l-l 



1 ^2m{rk/l-pk^/q) 



(51) 



k=0 
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Using the arbitrariness of the summation index (due to the periodicity of the exponentials), 
we can formally relabel this relation using k ^ k — 1 and write 



,27ri{r{k-l)/l-p{k-lf /q 



-2nir/l —2TTip/g 



IE 



(rk/l+2kp/q-pk /q 



(52) 



^^-2mr/l g-27rjp/gj 



,2iTi(k/l{r+2pl/q)-pk /q 



or 



h^i = e 



2m{r/l+p/q) ^ 



where 



r + 



2pl 



(53) 



which recursively relates the expansion coefficients. This implies that all of the (non-zero) 
\br\ have the same magnitude and, for odd values of q (and hence I), it gives 



1 



for all r = 0, .., g — 1 



(54) 



and each of the / = q 'mini-packets' contains 1/g of the total probability, overlapping with 
the original t = packet in A{t) with a periodicity Td/q. 

As an example, consider the case of Tret,/3 where p = 1 and q = I = 3. In that case we 
use Eqn. (I53j) . once with r = and r' = 2 to obtain 



62 = e'"*/=^ 60 

and again with r = 2 and r' = 4 < — > r' = 1 (recall the periodicity) to find 

bi < > 64 = 62 • 

The explicit value of 60 is obtained by direct summation of Eqn. giving 



(55) 
(56) 



2 

60 = - V e-2-^'/3 = - [1 + e-''-'/' + e-S'^*/^] = 1 [1 + 2e-''-'/'] = . (57) 



k=0 



Thus, near t ^ Trev/3, we find that 

ijj{x,t ^Trev/3) = bo^ci{x,t) + bitpci{x,t + Td/3) + b2'4'ci{x,t + 2Tci/3) 

[M^, t) + e^^'l'' {V'd(x, t + T,z/3) + V'd(x, t + 2T,,/3)}] 



(58) 



I 

73 



as first observed by Averbukh and Perelman in Ref. (j34D . The autocorrelation function near 
Trev/3 is also shown in Fig. El for our two model systems, where the exact and approximate 
realizations of this behavior are apparent. 

The arguments for even q are similar, but differ slightly depending on whether g is a 
multiple of 4, as discussed in Appendix O The result however, is quite similar in that one 
finds r = q/2 copies of the classical wave packet, separated by multiples of Tc//r, each with 
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maximum I -Aft) P = 1/r. The cases of g = 2 and g = 4 are considere d ag ain explicitly in 
Appendix O and are shown to reproduce the results in Eqns. (|HH|) and ()43|) respectively. 

Averbukh and Perelman used their results to correctly analyze the numerical calculations 
of Parker and Stroud (|28h. identifying the revival at T^e^u/S and fractional revivals of order 
1/4,1/6,1/8 which had earlier been described only as "...a complex pattern of quantum 
beats." We reproduce the results of their analysis of the (simulated) data in Fig.l^ Fractional 
revival structures of order up to 1/7 have been observed in other systems (|33h. as shown in 
Fig. El 

The observability of higher-order revivals, which would have r distinct features in |v4(t)p 
separated by Td/r (for r = q (r = q/2) for odd (even) q) depends on the level of the 
incoherent 'background' in Eqn. (jSUI) . so that in general such features would not be observable 
if 

- = \A{t ^ pTrev/q)? < lAn^r = > ^ ^ or q> ^n2^ (59) 

so that larger An wave packets can 'resolve' higher-order fractional revivals, as noted in 
Fig. 121 The behavior at such times during the collapsed or incoherent phase is shown in 
Fig. 121 for our two model systems for a multiple of Trev/37 where the dashed horizontal line 
indicates the value of 1/g, while the |v4(t)p values are seen to oscillate instead about the 
incoherent value, shown by the dotted line. 

We note that a number of authors have considered generalizations of these ideas (|67h. 
applications to specific model systems (16S ) - (75), and extensions to Hamiltonians that are 
time-dependent through a slow change in a parameter (|76l). 



E. Superrevivals 

For systems with purely quadratic energy dependence on a single quantum number (such 
as the infinite well, rigid rotor, and others), there are no independent time scales longer 
than the revival time, and the pattern of fractional and full revivals will repeat itself in- 
definitely with the Trev time scale. For more realistic systems, with higher order terms in 
the expansion in Eqn. (jH}, the superrevival time, Tguper becomes important, and Bluhm and 
Kostelecky (ffzl ). have exhaustively analyzed the long-term time-dependence of wave packets, 
with emphasis on Rydberg atom applications. They find qualitatively new patterns of re- 
vival behavior, with periodicities in the motion of the packet characterized by periods which 
are fractions of Tj-ev, giving a self-similar structure to the auto-correlation function plots for 
t > Trev and t < Trev The wave packet behavior on the Tguper time scale is similar to that 
of the fractional revival structures seen on the Trev scale, with integral multiples of Tsuper/^ 
appearing prominently, and smaller time scale periodicities of "iTrev appearing, explicitly 
due to the presence of the third-derivative term in Eqn. (jSI). In performing their analvses, 
they necessarily had to generalize and extend the results of Averbukh and Perelman (34) to 

include the contributions of the ^-'^'^i-pi'^-^of h phase terms. Extensions to even more time 
scales (j63), (ffsh . with possible applications, have also been discussed. 



F. Revivals in systems with two or more quantum numbers 

The generalization of the discussion of classical periodicity and revivals (fractional or 
otherwise) to systems with more than one quantum number has been presented by several 
groups ga), jll), ®. 

For the case of two quantum numbers, for example, one first assumes a time-dependent 
quantum state of the form 

m = E «(ni,n.)«K,n.)e-^^("-"^)*/^ (60) 
ni,n2 
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where the coordinate labels have been suppressed, and could be either two position-space 
(e.g., {x,y), {r,9)) or two momentum-space ((p^.,py)) values. The energy eigenvalues are 
not assumed to factorize, but for localized wave packets, we do assume that they can be 
expanded about a central value {ni,n2), with the resulting expression (to second order) 
given by 



1 _ x2 -E'('"'l,'"'2)A 1/ _ n2 -E'(?^l, ^2 



+ ^(^i-^i)M T2 +i(^2-n2)M ^^"^r" (61) 
+ (ni - ni)(n2 - na) — - — H . 



dn\dn2 



(ni,n2) 



In the spirit of Eqn. (jSI), we can use this expansion to define important time scales. Two 
separate classical periods are given by 

- \dE/dn^\ - |ai?/9n2| ^ ' 

where we will henceforth suppress the subscripts indicating that the partial derivatives are 
evaluated at (ni,n2). The corresponding longer-term, revival times are defined by 

{l/2)\d^E{ni,n2)/dnl\ ' {l/2)\d^E{ni,n2) / dnl\ ^ ' 

and the mixed term 

\d^E{ni,n2)/dnidn2\' ^ ' 

There can be recognizable periodicities in the short-term semi-classical time-development, 
with the two classical periods beating against each other, with the most obvious case being 
when the two periods are commensurate with each other, namely when 

y(ni) 

= T (65) 

rp{n2) 5 ^ 

cl 

where a, h are relatively prime integers. (Similar results hold for action-angle variables in 
classical systems fjSOh . (|5ll).) A simple example of such behavior is for a two-dimensional 
oscillator with differing frequencies, namely 

V{x^y) = '^{u:y + uy) (66) 

with quantized energies 

E{n^, Uy) = {n^ + l/2)t\Wx + {riy + l/2)t\Wy (67) 
where the condition in Eqn. (|b5|) reduces to 



T-i(ni) 



T^^^) uj^ b 



(68) 
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which is the famihar condition for the existence of Lissajous figures. Similar analyses will be 
shown below to lead to the appropriate conditions on closed or periodic orbits in simp le two - 
dimensional quantum billiard system s such as the square ( or rect angular) case in Sec. ()1V.A|) . 
the equilateral triangle (Sec. llV.C|) . and circular (Sec. IIV.D|) footprints being the most 
familiar. Nauenberg (|40j) has also considered the patterns of classical periodicities in systems 
with several quantum numbers. We note that this connection of the quantized energy 
eigenvalue spectrum to the structure of classical recurring orbits is quite different from that 
found in another semiclassical approach, namely periodic orbit theory ffSl. ). (82) . 

The extension of these ideas to central potentials in three dimensions, where the energy 
eigenstates do not depend on the azimuthal quantum number {E{nr,l,m) = E{nr,l) only), 

implies that one of the corresponding classical periods, T^l^^ — oo, is irrelevant, related to 
the fact that classical orbits are planar for all times. In a similar context, certain special 3D 
central potentials depend only on special combinations of quantum numbers with important 
semi-classical connections. For example, for the 3D isotropic harmonic oscillator, defined by 

V{r) = V{r) = —-{x' + y' + z') = ^r\ (69) 

the quantized energies can be written in the form 

E{N) = {N + 3/2)huj (70) 

where can be written in terms of Cartesian quantum numbers in the form = n^+ny+Hz 
or in cylindrical {nr,m^,nz) or spherical language {nr,l) language. The fact that there is 
only one effective quantum number is related to the symmetry properties of the oscillator 
potential (and its factorizability in many coordinate systems) and implies that all classical 
periods are the same for commensurate). This fact, in turn, is related to the result of 
Bertrand's theorem (jsj) which states that the only 3D power-law potentials for which all 
orbits are closed are given by = 2 (the oscillator) and k = —1 (the Coulomb potential). 
For the latter problem, the quantized energies also exhibit special patterns of degeneracy 
(dependent on 5*0(4) symmetries beyond the simple SO (3) rotational invariance) so that the 
Coulomb energy spectrum (with Z the nuclear charge and fi the reduced mass of two-body 
system) 

depends only the principal quantum number, n, and not on the radial and angular quantum 
numbers, n^, /, separately. The factorizability of this problem in parabolic coordinates (jH3) 
leads to the same result for the quantized energies, and is useful for problems involving the 
addition of an external electric field (the Stark effect.) 

The general time-development of wave packets which depend on two independent quan- 
tum numbers, on longer time scales, is then determined by the interplay of the three revival 
times, and the corresponding commensurability condition required to observe revivals is 

7^(ni) — ( 2.\ 7^("2) — ( ^ 1 7^(ni,n2) ('7'y\ 

\d) \ f J 

where c, d and e, / are again pairs of relatively prime integers. Fractional revivals are also 
possible if hierarchies such as 

Tfra. = (^) T(-) = (^) T(-) = T(--) (73) 



exist and Bluhm, Kostelecky, and Tudose (jTOD have examined the structure of such fractional 
revivals, extending the results of Ref. (34) to the case of additional quantum numbers. We 
will provide examples of systems described by these time scales in the context of several 2D 
quantum billiard systems below. 
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III. MODEL SYSTEMS 



The time-dependence of localized quantum wave packets, including possible quantum 
revival behavior, has been discussed for a large number of pedagogically familiar, and phys- 
ically relevant, one-dimensional model systems. We briefly review several such cases, and 
then focus attention on the infinite well as a benchmark case where exact quantum revival 
behavior is found. 



A. Free particle wave packets 

The analysis of Gaussian free particle wave packets goes back at least to Darwin 
and is a staple of introductory textbooks. While such systems do not exhibit quantum 
revivals (or even classical periodicity), it is useful to briefly review the basic formalism of 
Gaussian wave packet solutions to the free-particle Schrodinger equation in one dimension. 
It helps establish notation for later use, as well as providing an example for comparison with 
more realistic systems, especially illustrating important aspects of wave packet spreading in 
explicit form. 

A general Gaussian free-particle wave packet, with arbitrary initial values of {x)q = Xq 
and {p)o = po, can be written in momentum space as 



^ g-a^(p-Po)^/2 g-ipxo/;i 



The various expectation values related to momentum are given by 

1 1 
(p)t = Po , kP^)t = Po + TT^ ' and Ap^ = Apo = — ^ (75) 

so that the momentum-space probability density does not disperse in time, and the expec- 
tation value of kinetic energy, given by (jp')t/2m, is obviously constant, as it should for a 
free particle solution. 

The corresponding solution in position-space is given by 

Jj(rji. — ^ pipo(x-xo)/h p-ip^t/2mh p-jx-XQ-pot/mf /2{ahf{l+it/to) /-^g\ 



where 



2 mh 2m(Axo)^ 



tn = mha = — -- — -rr = (77) 

2 (Apo)2 h ^ ' 

is the spreading time. The corresponding position-space probability density is 

P(x,t) = |V^(x,t)p = ^e-(— W)^/^? (78) 

where 



x{t) = xo + Pot/m and (3t = ah\/ 1 + t^tg (79) 
and the time- dependent expectation values of position are 



{x)t = x{t) = xo+pot/m, (x^)^ = [x(t)]2 + ^, and Axt = , (80) 



— , and Axt = , 
2 ' V2 
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which exhibits the classical constant rate of change for the expectation value (consistent 
with Ehrenfest's theorem) and the standard spreading. 

For comparison to later examples, we note that the autocorrelation function for this class 
of solutions is easily obtained in either p- or x-space from Eqn. and is given by 



exp 



2Ul 



or 



{t/2to) 



exp 



-2a'pl- 



tt/2to) 
{t/2t,f 



+ (t/2to) 



(81) 



(82) 



As expected, the only relevant time scale is the spreading time (actually 2to) and we note 
that while for times satisfying t « 2to there is an increasing exponential suppression of the 
overlap between ipt and ipo, the exponential factor does 'saturate' for long times, giving 



\A{t » 2to) 



2to 
t 



exp 



smce a 



V2Apo ■ 



(83) 



The asymptotic form of the exponential factor can perhaps be best understood by noting 
that the 'distance in position space' between the initial 'peak' at {x)q = xq, and that at 
later times when {x)t = Xq +pot/m, grows linearly with t, while for long times the position 
spread, 

Axt = Axo^/l + {t/t^Y Axo^, (84) 

increases in the same way. This leads to factors in the exponent of the form 



{<t)-xm' 



(pot/m) 



{Axoit/to)y 



Potp 
uiAxq 



(poa) 



smce 



ah 



We distinguish the exponential suppression factor in Eqn. (|H^ (which we can describe as 
'dynamical' as it depends on the initial wave packet parameter, po) from the more intrinsic 
pre-factor term (containing only the spreading time) which is due to the natural dispersion 
of the wave packet (which we can therefore describe as 'dispersive'.) 

We note here that the Gaussian packet in Eqns. (fTij) or (fTBj) clearly satisfies the general 
bound on A{t) in Eqn. (Q) in Sec. III.Al saturating it to order (9()f:^) for short times. In this 
case, one has 



2a^ 



and 



^ ' ^ 2m 



+ a2 + 4^4 



(86) 



which give 



{AHf 



1 
2m 



pI + 



4a2 



(87) 



B. Wave packets and the constant force or uniform acceleration problem 

The problem of a particle under the influence of a constant force is one of the most 
familiar in classical mechanics, but is less often treated in introductory quantum mechanics 
texts, especially in terms of time- dependent solutions, despite the fact that closed form 
solutions have been known since the time of Kennard (j2i). For that reason, we briefly review 
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the most straightforward momentum-space approach to this problem. This system also 
provides another example of the behavior of A{t), this time in a more dynamical system. 

In this case, where the potential is given by V{x) = —Fx, we can write the time-dependent 
Schrodinger equation in momentum-space as 



dp 



ct>{p,t) = ih^^^ (88) 



or 

,d(j){j>,t) d(j){j>,t)\ p2 



We note that the simple combination of derivatives guarantees that a function of the form 
— Ft) will make the left-hand side vanish, so we assume a solution of the form t) = 

$(p — Ft)(p{p), with arbitrary and to be determined. Using this form, Eqn. (jH^ 
reduces to 

with the solution 

= e-'P'/'^'"™. (91) 
We can then write the general solution as 

(f){p^ t) = (!){p-Ft)e-'P'/^"'^^ (92) 
or, using the arbitrariness of $(p), as 

(f)(p^ t) = 0o(p-Ft)e'«P-^*)'-P')/'''"^'') (93) 

where now (/)o(p) is some initial momentum distribution, since (/)(p, 0) = </'o(p)- (We note that 
because the p^ terms cancel in the exponential, we will be able to explicitly integrate Gaus- 
sian type initial momentum-space waveforms.) The momentum-space probability density 
now clearly satisfies 

\(^{p,t)\^ = \Mp-Ft)\^ (94) 

which demonstrates that the momentum distribution simply translates uniformly in time, 
with no change in shape. 

For any general initial 0o(p) we now have the time-dependent expectation values 

= {p)o + Ft (95a) 
(p')t = (p')o + 2(j9)oFt + {Ftf (95b) 



AP* = vAp^)^^^ = \/(P')o - (P)0 = (95c) 



{&)t = (x)o + ^ + f^ (95d) 
m 2m 



giving the expectation value 



W=(|;^ + "'W)=<0-F(f)„ (96) 

which, in turn, also agrees with a similar calculation of {E)t using E = ih{d/dt), all of which 
are consistent with a particle undergoing uniform acceleration. 
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Using the standard initial Gaussian momentum-space wavef unction, 0o(p) in Eqn. ()74p . 
as the (f)o{p — Ft) in Eqn. (|^. we can evaluate the position-space solution using the Fourier 
transform to obtain 



ip{x,t) 



^iFt(xo-Ft^/6m)/h ^i{po+Ft){x-xo-pot/2m)/h 



1 



a/ y^ah{l + it /to) 



X e 



-(x-(a;o+Pot/m+FtV2m.))V2(a^^)^(l+*t/*o) 



The corresponding probability density is given by 

1 



P(x,t) 



(x-x{t))V/3f 



where 



and 



x{t) = xo H h — 

m 2m 



{x)t = x{t) , {x )t = [x{t)] + , and 



Ax - ^* 



(97) 



(98) 



(99) 



(100) 



so that this accelerating wave packet spreads in the same manner as the free-particle Gaus- 
sian example. 

The calculation of the autocorrelation function can once again be done in either p- or 
x-space using Eqn. Q to give 



A{t) 



1 



exp 



{2tplt/mh-{aFt)\l + {t/2tof)) 
4(1 -it/2to) 



-iFt{xo-Ft^/6m)/h 



(101) 



and the same factors of 1 — it/2tQ as in Eqn. (jSip are obtained; this expression also reduces 
to that case in the free particle limit when F ^ 0, as it must. The modulus-squared is given 
by 



\Mt)\' 



exp 



2a'{p'o + {Ftoy{l + {t/2to 



v/1 + (t/2to)2 

and we note that this result can be obtained from Eqn. 



(t/2to) 



(102) 



pl^pl + {Fto)\l + {t/2tof: 



l + {t/2to)\ 
by the simple substitution 

(103) 



For this case of uniform acceleration, the wave packet spreading is identical (same Axt) as 
in the free-particle case, which can be understood by noting that the distance between two 
classical particles starting at the same initial location, undergoing the same force, but with 

slightly different initial velocities (or momenta, p[^^ — p'^^ = Apo) would be 



XA{t) - XB{t) = (xo + p\^h/m + Fi^/2m) - (xq + p'f't/m + Fi^ /2m) 



Apot 



m 



(104) 



which increases linearly with time, in exactly the same way as for the free-particle solutions 
(when F = 0). The 'distance' between the peaks in ipQ and ipt, however, eventually grows 
as t"^ so that the exponential ('dynamical') suppression in A(t) does not saturate, while the 
'dispersive' pre-factor is exactly the same as for the free-particle case. 
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C. Harmonic oscillator 



The harmonic oscillator provides the most straightforward example of a bound state 
system for which the periodic motion of wave packet solutions (especially Gaussian) is easily 
derivable. In this case, the initial value problem is perhaps most easily solved, especially for 
Gaussian wave packets, by propagator techniques (53) . In this approach, one writes 



+00 



where the propagator can be derived in a variety of ways 
form 



dx'^lj{x',0)K{x,x';t,0) (105) 

and can be written in the 



K{x,x';t,0) 



exp 



imuo 



2hsm{ujt) 



2TTihsm(ujt) 

Using the initial position-space wave function 

1 



{{x + {x') ) cos(ti;t) — 2xx') 



^P{x,0) 



^ipox/h^-(x-xo)^/2l3^ 



where (3 = ah, one can evaluate the time-dependent wave function in closed form as 



exp 



S{x,t) 
2(3L{t) 



where 



and 



L{t) = Pcos{ujt) + 



ih 

mujf3 



sin (cut) 



S{x, t) = —xl cos(co't) + 2xxo — x^ 



imuj(3'^ sin(u;t) 

n 



(106) 



(107) 



(108) 



(109) 



(110) 



2xQPQsm{ujt) 2if3'^pQX i/^^pg sin(ti;t) 



muo h mujh 

The corresponding position-space probability density can be written as 



(x,t)p 



1 



exp 



[x — Xq cos{ijjt) — po sin(u;t)/mu;) 



with 



and 



Ax. 



^\L{t)\ 

{x)t = Xq COs(co't) + 

\m\ 



Po sin (cut) 



V2 V2 



13"^ cos2(cut) + {h/mujl3fsm.\ujt) . 



(Ill) 

(112) 
(113) 



Thus, the expectation value evolves in accordance with classical expectations ()86[). while the 
width in position-space oscillates (from wide to narrow, or vice versa.) The momentum-space 
variables behave in a similarly correlated manner with 



{p)t = —muxosm{ut) + pocos^ut) 



W(V/5)^cos2(cut) + (mcu/5)2sin2(cjt) 



V2 



(114) 
(115) 
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For the special case of the 'minimum uncertainty' wave packet, where 



the width of the packet is fixed as 



/3^ = — = /3o^ (116) 



Axt = Axo = ^ (117) 



which is the same as the ground state oscillator energy eigenvalue state, but simply oscillates 
at the classical frequency, which is similar to the famous example first cited by Schrodinger 
^ or to later examples of coherent states (87). 

The total energy of the general solution can be written in the form 



iE)={^ + -.raco'4]+-^{^ + ^] ■ (US) 



Since this is also a bound state problem, one can expand the solution in Eqn. ()108|) in 
eigenstates as 

oo 

^(a;,t) = ^a„u„(x)e-*-^"*/^ where E„ = (n + 1/2) /icu . (119) 

n=0 

We can then note that 

(E) = ((n + 1/2))^ = ((n) + l/2)ncj (120) 

which for a;o,Po = and (3 = /Sq gives (n) = as expected, since in this limit we recover the 
ground state energy. The spread in n values is useful in discussions of wave packet behavior, 
and we note that 

{AEf = (E^) - {Ef = {{n + 1/2)2) _ ^ ^/2))2 = {Anf{nujf (121) 
and we find that 



(An)2 = \ 



2 / o\ "2 

pi r, f mujp ^ 



+ ^^ + ^-2 (122) 



which will be useful. In the xo,Po = and j3 = j3o limit, we find An = 0, as expected for an 
eigenstate. 

All wave packet solutions (Gaussian or not) of the harmonic oscillator can be shown 
(using, for example (l53f), the expansion in Eqn. (jll9|) ). to satisfy 

iP{x, t + mT,i) = (-1)™ ^(x, t) (123) 

with a similar result for 0(p, t) as well. This implies that the autocorrelation function will 
be periodic as well, with 

A{t + kT,i) = {-lfA{t). (124) 

The evaluation of |A(t)P for general values of (3, xq, and po is straightforward enough, 
but the resulting expressions are somewhat cumbersome, so we will focus on several special 
cases as illustrative. 



27 



Case I: Minimum uncertainty wave packets, /3 = Pq. 
In this case the evaluation of A{t) gives 



A{t) = a/ cos(c<jt) + i sm{ujt) 



exp 



^0 



+ 



2mujh 



[(1 — cos(co't)) — i sin(co't)] 



(125) 



where great simphfications have been made by noting that 

1 



cos(co't) — i sin (tut) 



cos (cut) + i sin (cut) . 



(126) 



(We note that a very similar expression arises in analyses of the macroscopic wavefunction 
for Bose-Einstein condensates (88) and the collapse and revival of the matter wave field 
has been observed experimentally (60).) Once again, the two important parameters appear 
together in quadrature, as in the uniform acceleration case, and we have 



\A{t)\' 



exp 



Xn 



+ 



pI 



/5n ' mujh 



{1 — cos(u;t)} 



(127) 



which clearly exhibits the expected periodicity. All of the suppression can be attributed 
to the 'dynamical' factors (those in the exponential, containing xq and po) as there is no 
'dispersive' pre-factor component for this constant width packet. 



Case II: Arbitrary (3, but xo,po = 0. For this case, the wave packet does not oscillate, but 
only 'pulsates' ff89.). and the time-dependent wave function simplifies to 



^{x,t) 



1 



exp 



— (x^(cos(co't) + {imu jS'^ / h) sin(ci;t)) 



y/{f3cos{ujt) + (yih/mujp) sin(cut)) L W^{.cos{ujt) + {ih/mu(3) sin(^t)) 
It is convenient to define the parameters 



(128) 



r = 



PI 



so that 



1 

r 



Pi 



(129) 



in terms of which the resulting autocorrelation function in this case has the very simple form 



or 



A{t) 



2 cos(a;t) — i(r + 1/r) sin(co't) 
1 

Vcos2(tut) + (r + l/r)2sin2(cut)/4 



(130) 



:i3ii 



all of which can be attributed to a 'dispersive' (but in this case periodic) pre-factor. 

We first note that in this case A{t) is invariant under the transformation r — > 1/r, in 
other words, the time-dependence is the same for both initially wide (/? > (3o) or narrow 
{P < Po) packets. Plots of |A(t)P over one classical period are shown in Fig. El where it is 
clear that the larger the deviation from the 'minimum uncertainty' wavepacket, the faster 
the wavepacket 'pulsates' away from its initial shape. It is also noteworthy that in this case 
|y4(Tc;/2)| = 1 so that the wave packet returns to its initial form (up to a constant complex 
phase) twice each classical period. This can be understood from the expansion of this wave 
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form in terms of energy eigenstates. In this case, where the parameters xq,po both vanish, 
one is expanding an even-parity function. In the eigenstate expansion, 

+00 +00 
ij{x, t)=J2 anM„(a;)e-^^"*/^ = e'^'^^/^ ^ a^^^(^)e— * , (132) 

n=0 n=0 

the Un{x) have parity (—1)", so that for an even-parity state, only the even (a2n) terms are 
nonvanishing and the n-dependent exponential factors in Eqn. ()132|) therefore oscillate twice 
as rapidly as in the general case. 

Finally, for the very special case where (3 = (3q [r = 1) as well, we recover the ground 
state energy eigenstate of the oscillator, with its trivial stationary-state time-dependence 
(ipoi^yt) = Uo{x) exp{—iEQt/h)) and Eqn. (jl3(J|) indeed reduces to 



/ = = e+»^*/2 (133) 

^ ^ y 2 cos{ujt) - 2i sin(a;t) ^ ^ 

as expected. 

Be cause of the specially symmetric nature of the potential, besides the relation in 
Eqn. p23j) . we also have 

i;{-x,t + T,i/2) = {-i)tP{x,t) and <l){-p,t + T,i/2) = {-t) <j){p,t) (134) 

so that half a period later, the wave-packet is also reproduced, but at the opposite 'co rner ' 
of phase space, namely with x ^ —x and p <-> —p; note that two applications of Eqn. ()134|) 
reproduce Eqn. ()123p . One can also show these connections using propagator techniques, 
provided one properly identifies the complex pre-factors, as described in detail in Ref. |53l). 

This type of 'mirror' behavior can be diagnosed using a variation on the standard auto- 
correlation function, namely 

A{t)= / tp*{-x,t)tp{x,0)dx= / (j)*{-p,t)(p{p,0)dp (135) 



00 



which measures the overlap of the initial state with th e 'acr oss-phase-space' version of itself 
at later times. Given the simple connections in Eqn. ()134p . we can immediately write, for 
Case I considered above, 



|i(t)|2 = exp 



(136) 



which is exponentially suppressed at integral multiples of T^, but unity at t = {2k + l)Tc//2. 
This type of anti- correlation function finds use in the study of wave packet revivals where 
quantum wave packets may reform near t = Trev/2, as in Eqn. (jHEI), but out of phase with 
the original packet. 

Finally, we note that these results can be extended to the case of the 'inverted' oscillator 
(corresponding to a particle in unstable equilibrium, i.e., at a local maximum, instead of a 
local minimum of the potential) defined by 

V{x) = -\m^'^x^ (137) 

One can make the substitutions 
ijp —up , to iuj , sm{cot) i smh.(ujt) , and cos(ti;t) — > cosh(a)t) (138) 
and we briefly discuss this in Appendix iDl 



29 



D. The infinite well 



1. General comments 



The one-dimensional infinite potential well is the most frequently encountered example 
of a ID bound state system, finding its way into every introductory quantum mechanics 
text. Given its important pedagogical role, and familiar solutions, it is not surprising that 
many aspects of wave packet propagation in gene ral, and quantum revivals in particular, in 
this system have been studied (jSSl). (jOOl ) - ()10Clf l. (Interestingly, the time developm ent o f 
quantum states in the infinite well was used as a debating point by Einstein and Born ()101l ). 
()102i ) in more general discussions about the nature of quantum mechanics.) 

In this review, we will examine both the short-term quasi-classical time-evolution of Gaus- 
sian wave packets, the spreading to a collapsed state, as well as the structure of the revivals, 
mirror or 'anti'-revivals, and fractional revivals using both the standard auto-correlation 

function, A{t) and a related 'anti-correlated' A{t). We also focus on visualizing the same 
phenomena through the time-development of expectation values (and uncertainties), in both 
position- and momentum-space, as well as emphasizing the visualization of such effects in a 
variety of ways, including the Wigner quasi-probability (phase-space) distribution and the 
use of what have been dubbed quantum carpets. It is worth noting that bound-state wave 
packets in this system were examined in detail by Segre and Sullivan (90) in 1976 where the 
existence of exact wave packets revivals was discussed very explicitly, with the authors also 
noting in a footnote that " We suspect that this almost periodic behavior of packet width is a 
general property of hound-state packets^ 

We begin by defining the familiar problem of a particle of mass m confined by the potential 

,0 for < X < L 
V{x) = { (139) 
oo otherwise 



which has energy eigenvalues and stationary state solutions given by 

pI^ h'^n'^r. 
2m 2mL 

and 



En = ^ = -^r—TT = n Eq where Pn = —r- for n = 1, 2, ... (140) 



2 /'mix\ 

Unix) = ^ -sin i^-^j . (141) 

We note that the energy eigenfunctions have a generalized parity property (about the mid- 
point of the well), namely 

Un{L-x) = y|sin = -y|cos(nvr)sin(^) = (-1)^^+^ «„(x) (142) 

so that for n odd (even), the eigenfunctions are even (odd) about the center {L/2) of the 

well. This connection suggests that we generalize the notion of the A{t) 'anti-correlation' 
function in Eqn. p35j) to the geometry of the standard infinite well by defining 

~A{t)= ! ip*{L- x,t)^{x,0)dx. (143) 



The classical period for this system, obtained from Eqn. Q, is 

2nh 2mL' 2L 2L _ 



Pr. 
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and Vn is the analog of the classical speed, giving a result in agreement with classical expec- 
tations. The revival time is then given by 



T = 



27Th 

\E''\/2 



2-Kh AmL"^ 



En 



(145) 



which is also consistent with Eqn. ()21|) for general ID power law potentials. Thus, for 
72 >> 1, one can easily arrange to have T^ev » T^. The superrevival and higher order 
terms in the expansion in Eqn. (jH)) all vanish, making this an ideal system to study exact 
revival behavior. 

The quantum revivals in this case are in a sense exact, since we have 



'ljj{x,t + Tr, 



E 

n=l 



-iEnt/h -iEnTrcv / h 



E 

n=l 



(146) 

for all times. The wave packet will also, however, reform itself at t = Trev/'^ (a so-called 
mirror revival or 'anti-revival') at a (possibly) different location in the well. To see this we 
note that 



ip{L -x,t + Trej2) = J2 anUniL - a;)e-^-^"*/?ie-iEnT,,„ 



/2h 



n=l 

oo 



t/h 



' -j^\n+lg— in TT 



(147) 



n=l 



-i){x,t) 



which implies that 

|V'(x,t + T,e„/2)|2= |^(L-a;,t)p (148) 

so that at half a revival time later, any initial wave packet will reform itself (same shape, 
width, etc.), but at a location mirrored (91) about the center of the well. The use of the anti- 
correlation function, A{t), defined in Eqn. ()143|1 . will then be useful not only for documenting 
the short-term, semi-classical periodicity of such packets, but especially in establishing this 
additional type of revival structure. This behavior is a special case of the more general 
Tret,/2 revival discussed in Eqn. (jHHj) . 

The behavior of the wave packet in momentum space for such mirror packets is also easily 
derived using the Fourier transform connection, in this case 



0(p,t) 

We simply note that 

(j){p,t + Trev/'2) = 



\j2Txh Jo 



^{x,t) e'P'^/^dx. 



(149) 



= [ i){x,t + Trej2)e'P''l^dx 

inh Jo 

ij{L-x,t) e'P^'^^dx 

1 



\f2jxh 
1 



v2vr^ Jo 

jpL/h 



_V2ttE 

-e'P^/^'(j){-p,t) 



ipiy,t) e 



-ipy/h 



dy 



(150) 
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so that 



|0(p,t + T,,,,/2)p = |0(-p,t)|2 



(151) 



and half a revival time later the initial momentum profile is also reproduced, except flipped 
in sign [p — > —p), so that the particle is also moving in the 'other direction'. 

It has been noted several times that the dynamical time dependence of wave packets in this 
system can be described by the free-space evolution of an infinite sequence of appropriately 
displaced i nitia l wave functions (j91|), (103), which can be derived using a 'method of images' 
technique ()104f ). 



2. Gaussian wave packets 



For definiteness, one typically considers Gaussian-like wave packets, corresponding ap- 
proximately to initial momentum- or position-space wave functions as in the t = limits 
of Eqns. ()74|1 and ()7(i|l . The expansion coefficients corresponding to such initial states can 
be well-approximated in an analytic form if we assume that the initial position-space wave- 
function, 

-^Cix, 0) = ^ e-(x-xo)V2b^ ^ipoix-xo)/h (^52) 



(where b = ah) is sufficiently contained within the well so that we make an exponentially 
small error by neglecting any overlap with the region outside the walls, and may thus also 
ignore any problems associated with possible discontinuities at the wall. In practice, this 
condition only requires the wave packet to be a few times Axq = b/ \p2 away from an infinite 
wall boundary. 

With these assumptions, we can then extend the integration region from the finite (0, V) 
interval to the entire ID space, giving the (exponentially good) approximation for the ex- 
pansion coefficients (j97| ) 



\Un{x)\ \i)G{x, 0)] dx 



— oo 



1 ] ^!^^gimTXo/L^-b^{po+mzh/Lf/2h^ _ ^-irnrxo/ L ^-b'^ {po-mrH/ L)^ /2h^-^ (^l^^^ 

2iJ Y La/tF 
because we can write 

sin = l(e^"--/^ - e--Wi)^ (I54) 

and we can perform Gaussian integrals such as 

-ax^-bx^^^ J^^byAa ^^55^, 
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in closed form. This explicit expression goes back to at least Born ()101l ) and in our context is 
very useful because it can speed up numerical calculation involving the expansion coefficients, 
such as the evaluation of A{t). It also accurately encodes the sometimes delicate interplay 
between the oscillatory pieces of the Gaussian (e"*^'"^'/'^) and the bound state (e='=*'^'^^/") 
wavef unctions, which can be difficult to reproduce in a purely numerical evaluation, and it 
does so in a way that is valid for arbitrarily large values of po, where the integrand would 
be highly oscillatory, as well as being valid even for small values of po (or no.) Finally, 
Eqn. (I153|) nicely illustrates how ipci^, 0) and the Un(x) must not only have an appropriate 
overlap in position space, but also must have an appropriate phase relationship between 
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their oscillatory terms. This phase connection leads to the exp(— 6^(po ± mrh/ aY /2fi^) 
terms, which can be understood from a complementary overlap in momentum space. 

With this approximation {xq,L — Xq » b), the normalization and energy expectation 
value conditions 

oo 

^|a„|2 = 1 (156a) 



n=l 



1 fc2 



n=l 



can be satisfied to arbitrary accuracy. For the illustrative numerical results presented here, 
we use the nominal values 

L = l, h=l, and 2m = 1 (157) 

and we will often denote po = rionh/L to define the central value of no used in the eigenstate 
expansion, which in the large n limit will also give 

1 . Anirh , .1/ L , ^ 

Apo = — p — so that An = = - — — . (158) 



a 



L {ah)7cV2 Axo27r' 



This is useful as it can provide an estimate of how many eigenstates must be included in 
a given expansion. In the limi t tha t no » An >> 1, we can use this expression (relating 
b oc Axo oc L/An) in Eqn. ()153|1 . which then reduces to the standardly used Gaussian 
expression coefficients in Eqn. ()24|) . 

For most cases we consider in this section, we will also use 

a = -^-^ so that Axq = 0.05 << L , Apo = 10 and An = — ^ 3 . (159) 
10^2 vr ^ ^ 

This choice of parameters also implies that the other relevant time in the problem, the 
spreading time, to, defined by 

and Trev /to >> 1 for any wave packet sufficiently localized to be contained in the well. 
While we will typically focus on wave packets characterized by rio >> 1, which is what is 
typically obtained in experimental realizations exhibiting wave packet revivals, we start for 
simplicity with po = packets. 



3. Zero-momentum wave packets 

We can make immediate use of Eq. p53|) by considering zero momentum {po = 0) wave 
packets; this case corresponds to placing a particle 'at rest' inside the infinite well potential. 
For such cases, the only natural periodicity in the problem is the revival time in Eq. ()145p . 
because there is no cor resp onding classical periodic motion. In this special case, the expres- 
sion for the ttn in Eq. ()153p simplifies even further to 



4677 _f,2„2^2/i2 . fUTlXo 



e-" " " sin — ^ (161) 
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which shows that for several special values of xq in the well, a number of the expansion 
coefficients will vanish for obvious symmetry reasons. 

For example, for x^/L = 1/2, all of the even (n = 2,4,6, . . .) coefficients are zero and 
the only non-vanishing terms in the expansion are the odd ones {n = 2k + 1) which have 
energies of the form 



E„ 



TT 



2mL2 



{2k + if = EoiAP + Ak + l) = Eo + 8Eo 



k{k + l) 



(162) 



The first term in Eq. (|162p contributes only to the same overall phase of the time-dependence 
of each term. The second term is of the form 8Eq times an integer a nd le ads to revival times 
that are 8 times shorter than the standard Trev = 2Tih/ Eq in Eq. (jl45p . We illustrate this 
behavior in Fig. d^a) where we also note that moving slightly away from xq/ L = 1/2 (as in 
Fig. mb)) removes these exact sub-revivals, while the required full revival is still present. 
(This was observed in Ref. (jOlf ) for any even parity eigenstate in the infinite well.) 

Similarly, for the cases of Xq/L = 1/3, 2/3, the a„ with n = 3k vanish, leading to special 
exact revivals at multiples of Trev/3 for these two initial locations and this behavior is also 
shown in Fig.[7I^c) where we find exact revivals at these shorter time intervals. For xq/L = 0.8 
(still far enough away from the walls to be considered reliably in this approximation), we 
also notice large partial (but not exact) revivals at 0.4Trev and 0.6Trev for similar reasons 
(because sin(4n7r/5) vanishes for n = 5k.) (We note that it is also possible to construct 
odd-parity wave packets that have different patterns of special revivals at other locations in 
the well.) In the same plots, we also indicate the values of the a„ expansion coefficients and 
note the patterns of vanishing a„ values for the cases of special symmetry. For these po = 
states, the distribution of n-values is far from Gaussian. 

We are more interested in semi-classical wave packets with po ^ and we next show, in 
Fig- El the effect of 'turning on' momentum values for an xq = L/2 wave packet. For pq = 
(top line, (a)), we have the special pattern of exact revivals at multiples of Trev/ 8 noted 
above, due to the vanishing of the even expansion coefficients (shown in the corresponding 
|a„p versus n plot in the right column). For a small, non-zero value (po = 37r ^Po, 
Fig.El^b)), only the exact revival at Trev remains, because the even expansion coefficients are 
no longer forced to vanish. The autocorrelation function decreases somewhat more rapidly 
from its initial value than in the po = case, because the particle is now slowly moving away 
from its initial position, in addition to spreading out. 

For still larger values of momentum, such as po = iOn (Fig. IH^c)), we see obvious evidence 
for the classical periodicity and th e first appearance of fractional revivals (jOll ) at rational 
fraction multiples of T^ev, as in Sec. III.DI The corresponding a„ now exhibit a more obvious 
Gaussian shape, with a spread. An, which closely approximates the result in Eqn. ()158p . 
For even larger momentum values (see the po = AOOn case in Fig. IHfd), for example), the 
classical period becomes much shorter than any obvious fractional revival time scale, and the 
shape of the expansion coefficient distribution is unchanged (same An), but simply shifted 
to higher values of n, as expected. 

In what follows, we will typically use this last case, uq = 400, in our visualizations. This 
implies that the appropriate relative time scales are 

— — = 2no = 800 and — = non 

J-cl J-cl 

so that much of the fractional revival structure is obvious, while a number of classical periods 
are present before significant spreading occurs. 
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4- Short-term, quasi- classical propagation 



We expect the short-term propagation of wave packets in the semi-classical limit we 
study to share many properties with bo th the classical system, as well as with the quantum- 
mechanical free-particle in Sec. Illl.AI For example, well before the initial wave packet 
nears one of the infinite wall barriers, we expect the wave packet to propagate with time- 
dependence as in Eqn. (fTUj). and with a decreasing auto-correlation function (Eqn. 
described by 



^/ree(^)|" 



VI + (t/2to)' 



exp 



(t/2to) 



l + (t/2to)2)J 



(164) 



The quasi-periodic wave packet in this confining potential must return to something like its 
original state at t = T^, so we expect also suspect that the behavior near that point will be 
approximately described by Afree{Tci — t). For later convenience, we will also find it useful 
to single out the dispersive pre-factor above as the 'envelope' and define 



I ^enn (^) I 



Vl + (t/2to)2 



(165) 



To see to what extent the free-particle autocorrelation function is relevant, we plot in Fig. IHl 
|y4(t)p over the first classical period, using Eqn. ()76|) for the Gaussian wave packet, and 
increasing numbers {N) of eigenstates used in the expansion (solid curves.) For comparison, 
we also plot \Afreeit)\'^ (dashed) and \Afr-eeiTci — t)P (dot-dashed) and find that they are 
indeed good representations of the numerically evaluated result for the Gaussian packet over 
much of the early and late parts of the classical period. The agreement improves as increasing 
numbers of states used in the expansion can more approximately reproduce the exponential 
suppression predicted by Eqn. fjl64|) . We also plot (as diamonds) values of |A(t)P from 
the generic expression given by Nauenberg (j4Qi) for the short-term quasi-classical behavior, 
discussed in Appendix [Bl and note that especially for this case where there are no higher 
time scales, it is an excellent approximation. 

Over longer times, we know that the wave packet will spread significantly and that 
A{t),A{t) will not be exactly periodic, but will decrease due to the dispersion, presumably 
described by Eqn fll65p . To confirm this, we plot both |A(t)P and |^(t)P over the first 
15 classical periods (corresponding to almost 15TcL/to = 15/7r k, 5 spreading times) in 
Fig. along with the 'envelope' function in Eqn. ()165j) and note that the initial decrease 
in the maximum values of the autocorrelation and anti-correlation functions is indeed well 
described by this prescription. 

Continuing in this vein, we next examine, in Fig. ^2 both |v4(t)p and |74(t)p, now 
evaluated at integral multiples of T^/ ove r the first 150 classical periods. In this case we note 
that the dispersive prediction of Eqn. (|165p is a good approximation for over 40 periods, 
while the anti-correlation function (evaluated at t = nTc/) slowly grows from its exponentially 
suppressed value at t = and eventually crosses the \A{t)\ curve. For comparison, we have 
plotted (horizontal dotted line) the value of |A(t)P which might be expected if the wave 
packet were an incoherent sum of eigenstates, with no correlations, given by 



\A{t)f 



0.089 



(for An = IO/tt) 



(166) 



The relationship between A{t) and A{t) at some of the first few fractional revival times is 
also obvious 
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5. Revivals, fractional revivals, and mirror revivals 

The longer term structure of both A{t) and A{t), showing evidence for revivals, fractional 
revivals, and mirror revivals, is illustrated in Fig. only the interval (0,Tret;/2) is shown 
as the plot is symmetric about Trevf^ for this system. Many of the locations of possible 
low-lying fractional revivals are shown as vertical dashed lines, as well as the incoherent 
value, l^jncl, from Eqn. ()166|) . We can then examine the detailed behavior of A{t), A(t) near 
various specific fractional revival times. 

We begin, in Fig. with t = T^ev/'^i the p/q = 1/ 2 (or mirror or anti-revival) time, 
where we know that the wave packet will reform (Eqn. (|148|) ) on the opposite side of the 
well, with momentum in the opposite direction (Eqn. The upper plots in Fig. IT!^ 

show the behavior of \A(nTci)\'^ over the entire revival time (upper left) evaluated at integral 
values of T^, with the mirror revival indicated by the arrow, as well as the typical behavior 

of both |A(t)P (solid) and |A(t)P (dashed) near T^ev (or t = as well) in the upper right 
corner. The lower right shows the autocorrelation and anti-correlation functions for one 
classical period on either side of Tre^,/2, showing how the anti-correlation function is unity 
at Tre„/2, with A(t) being Tci/2 out of phase with it, as described very generally in Eqn. (jSHI)- 
Finally, the position-space probability density is plotted (lower left) where the initial wave 
packet was located at xq/L = 2/3 (solid), and was given positive momentum (illustrated by 
the arrow), while the value at t = Tj-ev/^ (dashed) is indeed centered at the mirror location 
(L — xq) and a short time late r (A t = Tci/20, dotted curve) the packet is clearly moving to 
the left, consistent with Eqn. ()15ip . 

We next turn our attention, in Fig. El to an example of a fractional revival at T^ev/^- 
In this case, the plot of \ip{x, t)p in the lower left shows the pattern of three 'mini' versions 
(dashed) of the original wave packet (solid), moving in a highly correlated manner (note the 
time-development at a time Tcz/20 later, dotted curve, indicated by the arrows.) This is also 
reflected in the plots of A{t) in the region around t = T^ev/'^ where there is a periodicity of 
Tc//3 due to the 'mini' packets (labeled a,b,c) being correlated with the initial wave form, 

and the same shorter periodicity is also apparent in the behavior of A{t), as predicted in 
Eqn. 

The case of the T^ev/^ revival, shown in Fig. is interesting since in this case the 

maximal values of |A(t)|, \A[t)\ near this fractional revival are identical (1/2) and they are 
both in phase (see the lower right of Fig. ^1) Once again, the two 'mini' packets are highly 
correlated in position-space (one packet reforming at xq, the other aX L — xq as befits a 
mirror revival) and in momentum-space (with their momenta in opposite directions) just as 
in Eqn. (gSj). 

Finally, if we look at many other times at random which are not near a fractional revival, 
as in Fig. UHl we see much smaller values of |y4(t)|, |A(t)|, typically consistent with the 
incoherent limit of Eqn. ()lfifi|l . and position-space probability densities (as in the lower left 
of Fig. IT6|l which exhibit far less obvious structure, being much more consistent with rapid 
oscillations about a uniform or 'fiat' classical probability distribution, namely P{x) = 1/L 
(denoted here by the two horizontal arrows.) This is typical of a collapsed phase where 
the wave packet is described more by an incoherent sum of energy eigenstates, with the 
probability density averaging about a semi-classical value, as in Eqn. (p?T|) . which in this 
case is the trivial constant value. 

Since the phenomena of wave packet revivals (and fractional revivals) is first and foremost 
about the important phase relationships between the energy eigenstates, as well as between 
the individual 'mini'-packets, it is instructive to examine both the real and imaginary parts 
of A{t), both over the long term revival time scales, as well as near fractional revivals or 
during the collapsed phase. 

We begin by examining in Fig. El both the real and imaginary parts of A{t) as well as 
a parametric plot (in an Argand diagram) over half a revival time. The plots show A{nTci) 
sampled at multiples of the classical period (just as in Figs.fTSl-ITBj). illustrating the intricate 
phase relationship present. The dashed circle corresponds to the value of \A\'^ = 1/2 and we 
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indicate the location of a 1/4 fractional revival showing how the autocorrelation function 
'lingers' tangentially during the fractional revival period. Close-up versions of the fractional 
revivals at Trev/4: and T^ev/'i are shown (in the same format) in Figs. 1181 and 1191 respectivelv. 
where the obvious phase correlations between the 'mini'-packets are clearly present. For the 
case of the Trev/4: fractional revival, the approximate wavefunction in Eqn. PHj) implies that 
the autocorrelation function near that time can be generally written in the form 

Ait ^ Trev/A) = 4= [e*"/'A/(t) + e-'^'^Aei{t + T,i/2)\ (167) 
v2 

where we define 

/ + 00 
rci{x^)Mx,^)dx (168) 
-oo 

and we see that this form is consistent with the visualization in Fig. ^1 

For the T^ev/S case in Fig. [121 the wavefunction in Eqn. can be used to write 

A{t ^ TrejS) = {iA,i{t) + e-2-'/3 [A^t + T,i/3) + A,i{t + 2T,z/3)]) /Vs (169) 

which describes the magnitude/phase relationships seen in Fig. IT9l 

Finally, the behavior of these quantities during a typical time, T*, in the collapsed phase 
is shown in Fig. 1201 where in this case the dotted circle corresponds to the incoherent value 

\A. 12 _ I 14 
l^incl / J I "'Til • 



6. Expectation value analysis 

While an analysis of the short-term quasi-classical and long-term revival structure of 
wave packets using autocorrelation function methods is often the most directly comparable 
to important experimental observables, the visualization of the time-dependence of quantum 
wave packets through their expectation values, both in position- and momentum-space, can 
also be a valuable tool ( lQ5j for understanding many of the effects arising in the quantum 
mechanical time evolution of wave functions. For example, the time-evolution of the uncer- 
tainty principle product, Ar ■ Apr, has been used (jS^), (jsil) to distinguish between regimes 
of classical versus non-classical behavior in Rydberg atoms. 

In this section, using the same parameter values as in earlier analyses, we examine the 
expectation values, (x)t and {p)t, and uncertainties, Axt and Apt, for Gaussian wave packets 
in the infinite well (94), both for short times (several classical periods) and over the entire 
revival time. 

Once again, com parison s to both the classical limits and the results of the free-particle 
wave packet in Sec. IIII.AI can be useful. For example, the time-dependence of a classical 
particle in the infinite well potential, starting with velocity +vo, at the center of the well, is 
shown in Fig.|^ with cusps in x{t) at the 'bounces' at the infinite walls, and a discontinuous 
change in momentum (from ±mf o to ^mvo) at the same points. We might expect somewhat 
similar behavior, softened by quantum effects, at least at early times for the quantum prob- 
lem. For comparison, we recall that the free-particle Gaussian wave packet in Eqns. (f7^ 
and (f7B)) is characterized by 

{x)t=pot/m + xo , Axt = Axo^/lTJtJtof (170a) 
{p)t = Po , Apt = Apo. (170b) 

The short-term, quasi-classical time-development can be seen in Fig. |221for the position- 
space variables and in Fig. 123 for the momentum-space quantities, over the first 10 classical 
periods (approximately 3 spreading times.) The classical looking behavior of {x)t is apparent, 
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but with the maximal values near odd integral multiples of Tci/2 (corresponding to the 'hits' 
on the walls) becoming smaller and smaller, as the wave packet spreads, becomes wider and 
finds it increasingly harder to get near the wall. The results for Ax^ are also consistent with 
expectations, namely there are small 'dips' in Axt at t = {2n + l)Tc//2 arising from the 
'compression' of the packet ( jlOGl ) as it strikes the wall, superimpose d on a uniform increase 
(the dotted line) consistent with the free-particle spreading in Eqn. ()170a|) . 

In our exemplary cases, the initial wave packet is characterized by positive momentum 
Pq = 4007r ^ 1260 and we see distinct evidence of 'fiips' between ±po, at times consistent 
with the classical 'bounces' in the {p)t plot at the top of Fig. |2S1 The initial momentum 
spread from Eqn. ()159|) is Apo = 10 and we expect this to stay roughly constant between 
successive collisions with the wall. During the collision times, however, the spread arising 
from the momentum-space probability density is dominated not by the intrinsic spread of 
a single peak, but rather by the 'distance' between the two peaks at p = ±po- In the limit 
where P{p) can be approximated by two highly peaked features at ±po, namely 

Pip) = l[Sip-Po)+5ip + Po)] (171) 

the corresponding spread in momentum is actually closer to Apd = po, and this type of 
behavior is also evident from Fig. |2S1 

For the autocorrelation function analysis, we found that the approach to the collapsed 
state could be discussed in terms of the decrease of \A(t)\ from its highly correlated value 
of uni ty ne ar classical periods to something closer to the incoherent value of \Ainc\ given by 
Eqn. ()l(i(i|l . In more physical terms, we expect the increasingly wide quantum wave packet 
to approach a collapsed state which is uniformly spread over the entire well, consistent with 
the classical probability density Pd{x) = 1/L. For that distribution, we have position spread 
is given by 

{x)ci = ^ , {x^)d = 1 SO that Ax,i = ~ 0.288L. (172) 



To visualize this approach to the collapsed state, we plot in Fig. |^ the time- dependent 
position spread, Axt corresponding to three different values of Axq and note that they all 
'saturate' at Ax fiat = L/\/l2, at different times, due to the difference in spreading times 
(to = 2m(Aa:o)^/^) . In each case, the bound state values initially follow the free-particle 
prediction of Eqn. p7()aj) for short times. The time scale for this collapse can be estimated 
by equating 

-|= = Axd = Ax(t=T,„„) = Axov^l + (T,oji/to)2 (173) 

giving 



which defines another time scale, intermediate between to and Trev 

The same typ e of approach to the values appropriate for a collapsed state can also be 
seen (in Fig. 1^ for {x)t (approaching {x)ci = L/2), for {p)ci = (equal admixtures of left 
and right moving states), and for Apd = +po (same reason) over the same time scales, and 
with the same dependence on the spreading time (via its dependence on the initial width.) 

Finally, the behavior of selected average values over the entire revival time, illustrating 
both the revival and mirror revival times, as well as 'close up' views at intermediate times, are 
shown in Figs. OBI 1271 andl28l The plots are shown for times given by t = {n + l/8)Tci (solid) 
and (n-|-5/8)Tc/ (dashed) to help distinguish standard versus mirror revivals. The magnitude 
of excursions from the classical value of {x)ci = L/2 for the collapsed state are shown on an 
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expanded scale in Fig. 123 which includes a fractional revival at t = T^ev/^ = 133Tcl and 

t = Trev/5 = IQOTa. 

For the time-dependent spread in position, Axt, in Fig. 123 we see that Axt returns to the 
initial value, Axq, at Trev (full revival) and Trevf^ {p/q = 1/2 or mirror revival) as expected, 
but also at Trev I ^ and 3Tj.e„/4. These last two instances are special cases and are due to the 
choice oi Xo = L/2 for the initial wave packet (compared to = 2L/3 as shown in Fig. IT^ 
where the two 'mini' packets are superimposed at the same location (the midpoint of the 
well), but moving with opposite momenta. This is consistent with the {p)t and Apt values 
at the same times in Fig. |2H1 

We note that the observed agreement of a number of expectation valu es an d spreads with 
results obtained from the classical probability distributions (as in Eqn. ()172j) ) during times 
much longer than the spreading time, while consistent with the approach to a purely 'flat' 
distribution, can also be thought of, during periods of fractional revivals, as being due to 
the correlated behavior of a number of 'mini-packets' or 'clone' wave packets, which can give 
similar results ^10(1) . 



7. Phase-space picture of fractional revivals using the Wigner function 

Another useful way to visualize the correlated position- and momentum-space structure 
of wave packets in a ph ase-space type picture is by using the Wigner quasi- or pseudo 
probability density (jlOTf ) defined by 

1 j' + OO 

W{x,p;t) = — ij*{x + y,t)^{x-y,t)e^''Py'^dy (175) 



Txh 



oo 
oo 



1 r+oo 

= ^ (l)*{p + q,t)(j){p-q,t)e~^"^'/''dq. (176) 

TT/i J-oo 

(The Wigner distribution f or G aussian wave packet solutions in the infinite square well 
has been discussed in Ref. fllOSf ) to which we refer the reader for many calculational and 
visualization details.) The genera l pro per ties of the Wigner distribution have been discussed 
in a number of accessible reviews ()109[ ) - ()117f ). where, among other features, it is noted that 
integration of W{x,p]t) over one variable or the other is seen to give the correct marginal 
quantum mechanical probability distributions for x and p separately, since 

/ + CO 
W{x,p;t)dp = mx,t)\' = PQM{x,t) (177) 
-oo 

/ + 00 
W{x,p;t)dx = \(j){p,t)\' = PQM{p,t) . (178) 
-oo 

The Wigner function is also easily shown to be real, but need not, however, be positive 
definite, hence the name quasi- or pseudo-probability density. 

The Wigner distribution for a standard time-dependent Gaussian free-particle solution 
of the form in Eqns. (f7^ or (fTBj) is easily obtained, using standard Gaussian integrals, and 
is given by 

W^^\x,p; t) = — e-'^'(P-Po)' Q-ix-xo-pt/mf/f^^ ^^7g) 

irh 

where P = ha; in this case, the ultra-srnooth Gaussian form does give a positive W{x,p; t), 
as has been discussed in the literature (llir 
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Another result which will prove useful in what follows is the expression for the Wigner 
function for the case of a linear combination of two ID Gaussians, characterized by different 
values of Xq and Pq. For example, if we assume that at some instant of time we have 



= 'yip(^G){x;XA,PA) + 5ip(G)ix;XB,PB) 

1 



;i80) 



7 



IT 



+ 6 



^-(x-XBf/2l3^ ^ipB{x-XB)/fi 



the corresponding Wigner function is given by 
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where 
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X 



Xa + Xb 



and 



p 



Pa+Pb 



;i82) 



2 "2 

In this case, the Wigner function is characterized by two smooth 'lumps' in phase space, 
corresponding to the values of {xa,Pa) and {xb,Pb) of the individual Gaussians, but also 
by an oscillatory term, centered at a point in phase space defined by the average of these 
values; the oscillations are either in x (if pa 7^ Ps), p (if xa 7^ x p) or both. (See Ref. ()119f ) 
for a similar result for the case when Pa = Pb = 0] see also Ref. (|63h for a related expression 
involving coherent states.) 

For a general, time-dependent wave packet constructed from energy eigenstates of the 
form in Eqn. (jH}, the Wigner distribution is given by 
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where, in general, we must calculate both diagonal (m = n) and off-diagonal (m 7^ n) terms 
of the form 



W^"''''\x,p) 



1 



+00 



Um{x + y) Un{x - y) e^'Py/^dy 



(184) 



where we have assumed that the individual position-space bound state eigenf unctions, Un{x), 
can be made purely real. The off-diagonal Wigner terms are, however, not real, but do satisfy 



[iy('"'")(x,p)]* = (a;, p)- 

For the infinite square well, the diagonal (m = n) terms are given by 



Un{x + y) Un{x - y) e 



2ipy/h 



dy 



(185) 



;i86) 



and the limits of integration are determined by the restriction that the m„(x±?/) in Eqn. (jl4ip 
are non-vanishing only in the range (0, L) and they must therefore simultaneously satisfy 
the requirements 

0<x + 2/<L and Q<x-y<L. (187) 
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This leads to upper and lower bounds for the integral over y in Eqn. ()186|) which depend on 
X via 



— X < y < +x 
-{L - x) <y < +{L - x) 



for < a; < L/2 
for L/2<x<L^ 



(188) 
(189) 



Thus, over the left-half of the allowed x interval, (0, L/2), we have 
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nn{x + y) 



sm 



mT{x — y) 
L 
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A{p/h - mx/L) 



4:{p/h + mT/L) 



(190) 



f 2mix\ sm(2px/h)^ 



while over the right-half of the interval, (L/2 , L ), on e simply makes the replacement x 
L — x. This form has been derived before ()115[ ) - (jllTf ). although in at least one reference it is 
written in terms of Bessel functions {jo{z)) which somewhat obscures its simple derivation. 

For a general wave p acke t solution for the infinite square well, we also r equi re the off- 
diagonal terms in Eqn. ()184|) . Using the position-space eigenstates in Eqn. ()141|) . over the 
interval (0,L/2) we find that 
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2pL/h + [m + n)n 

-i{ra-n)nx/L ?>V^[{2p / k - {m + n)TX / L)x\ 

2pL/h — (m + n)7r 
_^+^i^+n)nx/L sm[{2p/h+ {m-n)7i/L)x] 

2pL/h+ (m — 71)71 
_^-i{m+n)^x/L sin[(2p/;i- (m-n)7r/L)x] 
2pL/fi — (m — 71)71 

and it is easy to check that this result reduces to the expression in Eqn. ()19U|) when m = n. 
In order to extend this to the interval [L/2, L), it is important to note that the substitution 
X — *• L — X should be made only in those terms arising from the integration over dy, namely, 

the sm[{2p/h± [m ± n)/L)x] terms. 

Using these e xpres sions, the expansion in Eqn. ()183p . and the Gaussian expansion coef- 
ficients in Eqn. ()153|) for the an, we can evaluate the time- dependent Wigner function for 
the infinite well. As examples, we show in Figs. EUl and IHUl W{x,p;t) versus {x,p) at two 
fractional revival times, namely Trev / 3 and Trev/4:, for direct comparison to the position- 
space probability densities in Figs. and ^1 for the same values. (The same initial values 
of xq are used in each case, but we use no = 40 here for ease of visualization.) The results 
are consistent with the expression in Eqn. with smooth isolated Gaussian 'lumps' and 

oscillatory cross-terms; for example, in Fig. 1291 the cross-term between the two Gaussians 
which have the same momentum {—A07rh/L ^ —126 in this case, so pa = Pb), but different 
central locations [xa Xb) is seen to be oscillatory only in the p variable. These views make 
clearer the Schrodinger-cat like behavior of the split wave packet at fractional revival times. 
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8. Quantum carpets 



A semi-classical visualization of the initi al spre ading of bound state wave packets in 
the infinite well was provided by Born ()101f l. ()102f l who examined the diverging classical 
trajectories of particles with differing initial velocities (momenta) in (x, t) plots where such 
classical paths can be identified with individual 'world- lines'. An example of such a plot 
is shown in Fig. |^ which illustrates how the spreading can lead to an almost uniform 
probability density, giving rise to calculations such as in Eqn. ()172j) . 

The corresponding quantum mechanical picture of such behavior can be obtained by 
examination of the quantum mechanical position-space probability density as a function of 



time, given by plots of t) 



[X. 



t)P versus (x, t). As examples, we show in Fig. IH^ 



the probability density over the first classical period, comparing it to the more interesting 
structures present over the same time interval starting at the T^ev/^ fractional revival, as 
seen in Fig. 1221 

To visualize these over even longer time scales, we note that because of the periodicity 
of the system given by the revival time, it suffices to plot P{x,t) = \4'{x, t)\'^ over the two- 
dimensional {x,t) space given by (0, L) and (0,Trei;/2), such as in Fig. IMl which were first 
produced by Kinzel (120); such plots have come to be known as quantum carpets. The rich 
structure of ridges and canals apparent in such plots, which are clearly correlated with the 
spatio-temporal stru cture of reviv als and fractional revivals, have been discussed in a number 
of approaches (|9^. ()120l ) - ()128l ). including making use of a traveling wave decomposition 
of the wavefunction or u sing Wigne r quasi-probability densities. 

One useful approach (J123'), (Jl23) to the visualization of this pattern formation, based on 
an expansion of the wavepacket in traveling waves, begins by writing the stationary state 
solutions in the form 



2 /n'KX\ 
sin I 1 
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(192) 



so that 



^{x,t) = ^a„^„(a;,t) 



(193) 



n=l 



for any more general time-dependence state, including localized wave packets. The proba- 
bility density can then be written as 

P{x,t) = ij*{x,t)^{x,t) = P,i{x,t) + Pgc{x,t) (194) 
where one defines classical (Pd) and quantum carpet (Pqc) contributions by 



|g«(n-m)7r[x/L+2{n+m)i/Treu] _j_ ^—i{n—m)TT[x/L—2(n+m)t/Trev]^^-^g^'^ 



n,m=l 
oo 



^ n,m=l 



^i{n+in)n[x/L+2{n~m)t/Trev] _^ ^-i(n+m)TT[x/L—2{n—m)t/Trev]^ (196) 



The contributions to the classical component (Pd) are dominated by {x,t) locations at 
integral multiples of [x/L ± 2(n + ■m)t/Trev]- For localized wave packets, the n,m values 
will be peaked near no and using the fact that Trev = '^n^Td, the factors in the exponents 
can be approximated [x/L ± 2t/Td] which are therefore similar to those seen in Fig. I^for 
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classical trajectories. Since Td << Trev, or equivalently, since n,m » 1, these world-lines, 
when plotted over the ranges (0, L) and (0, Trev), have very 'flat' slopes and, also taking into 
account their diverging paths, soon cover the {x, t) plane almost uniformly, as illustrated in 
Fig. 123 providing a 'background'. 

The quantum carpet [Pqc) terms, on the other hand, have dominant contributions from 
terms in the exponentials of the form [x/L ± 2(n — m)t/Trev]- The (n — m) terms will only 
contribute significantly for n,m values for which \n — m\ < 0{An). This implies that the 
world-line slopes for these factors will cover a significant fraction of the entire {x,t) 'plane', 
in a few 'bounces', with two typical cases shown in Figs. EHl (b) and (c); the interference 
of such terms with the classical background can then give rise to the observed patterns of 
ridges and canals. 

For initial wave packets which aremore highly localized (smaller Axq, larger An), more 
obvious features can be resolved ()129f ). similarly to the dependence seen in Fig. 121 The 
relative 'sharpness' of the features can also be seen to arise from the differing pre-factors in 
the exponentials. For example, the {n — m) pre-factors in the Pd terms will be dominated 
by small values of \n — m\, so that relatively large changes in the corresponding [x ± 2{n + 
m)t/Trev] terms will still contribute significantly, enhancing the uniformity of the classical 
contributions; on the other hand, for the Pqc terms, the (ra+m) pre-factors are large ((9(2?7,o)), 
so that small changes in the accompanying [x ± 2|n — m\t/Trev\ factors will rapidly lead to 
cancellations, helping to define the 'contrast' seen in the quantum carpet images. 

It has also been observed that the 2D (in (x, t)) probability density patterns observed 
here are similar to what have been described as 'scars' ()13flf ) in spatially two-dimensional 
{{x,y)) quantum mechanical systems. We also note that si milar techniques have been used 
to examine other systems which also exhibit exact revivals ()125[ ). 



E. Variations on the infinite well 



Given the pedagogical familiarity and simplicity of the infinite well potential, a number 
of variations on this system have been studied and we briefly mention three of them. 

We have noted that the important time scales for wave packet propagation are determined 
by the energy eigenvalues through Eqns. ©, (fT^- and (j^, so it is clear that systems with 
identical (isospectral) or closely related energy spectra will have similar patterns of classical 
periodicity, wave packet collapse, and revival. For example, pairs of systems defined by 
superpartner potentials which are related by one-dimensional supersymmetry (131) - (13^1) 
have energy spectra which are identical except for the ground state energy and would be 
expected to exhibit almost the same pattern of revival behavior, at least for n^ >> 1. The 
construction of superpartner potentials (134} starts with the assumption of a ID Hamiltonian 
with a zero-energy ground state, namely one such that 

Hi-)M^)=[-^^,+V^-)ix)^M^) = 0- (197) 

Since ipo{x) is assumed known, one can define a superpotential for the problem by the 
identification 

Wix) ^ ( . (198) 



^2m y-ipoix 

The pair of systems defined by the Hamiltonians 
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can then be shown to have the same energy level spectrum, with E^^^ = except 

that the zero-energy ground state of {Eq~^ = 0) has no counterpart in H{+)- The 
wavefunctions in the two systems are related by generalized raising and lowering operators. 

The standard infinite well potential of Eqn. ()139|) can be put into the form of -f^(-) by 
subtracting the appropriate zero-point energy to write 

y(_)(x) = na^)-|^ with 4-)(x) = y|sin(^) (200) 
with quantized energies given by 

The superpotential is then given by 

iy(x) = --^^cosf^) (202) 



/2mL V L 

and the superpartner potential of the infinite well potential is given by 



2mL2 



, TTX\ 

2 csc^ ( — j 



(203) 



so that large n wave packets in this potential will have the same pattern of semi-classical 
periodicity, and exact revivals, as for the infinite well. We note that this pot entia l is a 
special case of a larger class of Poschl- Teller potentials and the authors of Ref. ()135f l have 
constructed coherent states (see also Refs. ([TU l - (fl3h) for such systems, making using of 
their SU{1, 1) symmetries, while the r evival structure in Poschl- Teller and Rosen-Morse 
()125l ) and other related potentials (13&} have also been considered. Revivals and fractional 
revivals in other SU{1, 1) symmetric systems have also been discussed f(137.) . 

We note in passing that the same approach can be used with the harmonic oscillator 
potential, defining 

= - — (204) 

but trivially returns the supersymmetric partner potential 

Vi+){x) = ^x' + — (205) 

and the raising/lowering operator formalism is the standard one seen in textbooks. The use 
of a phenomenological supersymmetry involving atomic energy levels, has been discussed 
()138l ) as it relates to the revivals and fractional revivals of Rydberg atoms and radial squeezed 
states. 

The finite well has been studied by several groups ()139l ). ()140[ ) with an eye towards pro- 
viding a more realistic model system, focusing on the appearance of superrevival time scales 
since the energy eigenvalue spectrum is no longer exactly quadratic. Aronstein and Stroud 
(fT41i) have provided a useful description of the bound state energy eigenvalue spectrum 
which is helpful for the determination of both the global energy spectrum as well as local 
approximations, useful for wave packet construction. Finally, the addition of a 5-function 
to the infinite well problem is a staple of the pedagogical literature ffl4Z) . an d its effects on 
the revival structure of wave packet development has been discussed in Ref. ()143l ). 
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F. The quantum bouncer 



1. Energy eigenf unctions and eigenvalues 

One of the more famihar systems of classical mechanics, the 'bouncer', has been con- 
sidered in the context of quantum mechanical wave packet propagation and revivals. This 
system, defined by the potential 

, l-oo for 2 < 
V{z) = { - , (206) 

Fz for z > 

corresponds to a particle under the influence of a constant force F (such as gravity), but with 
an infinite wall at ^ = from which the particle bounces. While this system has long been 
a popular non-trivial example in the pedagogical literature and in collections of problems in 
quantum mechanics, advances in experimental techniques ()144( ). using reflection from laser - 
i ndu ced evanescent waves, have allowed the observation of such 'bouncing' atoms ()14r 



f l47f 



( 146 ) and even the quantum mechanical interference between different atomic trajectories 



The classical periodicity for this system is simply twice the time a particle takes to fall 
from its maximum height, Zo, and is given by 



(207) 

and several groups ()148f l - ()151f l have considered different aspects of the time-development 
of (Gaussian) quantum mechanical wave packets, initially centered at zq, including their 
subsequent revival behavior. (The effect on the bouncer of the application of a pe riodi c 
external driving force, especially on the structure of revivals, has also been discussed (jl52f ). 

diss).) 

The Schrodinger equation for this system can be written as 

' +FziP{z) =Ei){z) (208) 



2m dz^ 



and the change of variables 



py + a, where p = i ^— ^ j and ^ = (^09) 



reduces it to the familiar Airy equation 



yi,{y) . (210) 



The condition that the wavefunction be well-behaved as z/y +oo excludes the divergent 
Bi{y) solution, while the boundary condition at the infinite wall, il){z = 0) = 0, implies that 

Ai (^--^ = = Ai{-yn) (211) 

and the quantized energy eigenvalues are determined by the zeros (?/„) of the well-behaved 
Ai{y) solution. (This problem is discussed in a variety of textboo ks, b ut the first analysis 
using the Schrodinger equation seems to have been done by Breit ()154f ).) 
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The energy spectrum then has the form 



f^2p2\ 1/3 



E^ = yJ^] . (212) 
\ 2m J 

The Airy function zeros (|/n) can be derived from standard handbook results, but the energies 
in the large n limit relevant for wave packet construction can also be easily obtained from a 
WKB approximation. The WKB quantization condition is 

J ^2m[En -Fz)dz=\n^-^ - J fnx (213) 

where we use the appropriate matching coefficients, Cl = 1/2 for the infinite wall, and 
= 1/4 for the linear barrier. The resulting energies are given by 

/Stt \ fh^F'^\ 

En= [—{n + ^/A)\ f^j where n = 0,1, 2,..., (214) 



which agrees with the Airy function analysis ()149f ) to this order in n. (We note that ex- 
perimental evidence for the q uanti zed nature of the bound state spectrum of neutrons in 
the Earth's gravitational field (155!) has recently appeared, motivating additional studies of 
wave packet dynamics (156) of such systems.) This expression can also be used to associate 
the central value of no in the eigenstate expansion with the initial position, zq, via 



Fzo = Einuiai = ( y('^o + 3/4) j ( ] (215) 



which will prove useful. 



2. Classical period and revival time 



The classical periodicity of the system, as encoded in the energy eigenvalue spectrum, is 
given by 

37r%o + 3/4) 



27ih 
\E'{no)\ 



E{no) 



2mzQ 



(216) 



(where we have used the Zq/uq identification in Eqn. (j215j) ) and is consistent with the classical 
result. As an example of the semi-classical, short-term time-development of a Gaussian wave 
packet in this potential, we show in Fig. EHl the position- and momentum-space probability 
densities for such a packet with 



using the parameters 



zo = 25 



h 



and 



Azo 



2m 



(217) 



(218) 



so that Tci = 10. On the left of Fig. |2ol the position-space probability is initially seen to 
spread in a manner which is numerically consistent with Eqn. (|100|) . while the calculated 
position value {z)t (solid curve) agrees well with the classical expectation for the trajectory 
(dashed) except, of course, for the cusp at the 'bounce'. The packe t exhi bits the standard 
'interference' pattern during the colhsion with the wall ffl06[ ). ()157l ). ()158l ) at the 'bounce', 
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and then reforms into something hke the initial packet (compare to the dotted initial packet 
superimposed on the t = 10 case), only wider. 

For the momentum-space distributions (shown on the right of Fig. ()36|) ). we also see 
featur es of both the classical motion and the uniformly accelerated wave packet, as in 
Eqn. (|95ap . The expectation value of momentum {p)t, calculated from |0(p, t)p and plot- 
ted as the solid curve, is once again consistent with the classical trajectory (dashed curve), 
except near the discontinuous, impulsive change in momentum values at the 'bounce'. The 
shape of the momentum-distribution follows the form expected from Eqn. ()94j). namely 
uniform translation with no change in shape, from t = 1 ^ t = 3 and then again from 
t = 8 — t = 10, that is, during the time when it is not in collision with the wall, but with a 
definite final change in shape, compared to the initial |(/)o(p)P superimposed on the t = 10 
result, resulting from the collision with the wall. The dotted vertical lines indicate the values 
of p = 0, and also the classically expected minimum and maximum values of momentum 

given by ±Pm = ^V'^fnE = ±\/2mF zn. T h e cha nge in shape can be understood, in great 
part, from purely classical arguments ()149l ). ()151f ). 
Using Eqn. (fl^ . the revival time is then given by 

_ Arch _ IQmzl 

\E"{n,)\~ nh ^ 

which can be compared to the similar result for the infinite well in Eqn. ()145j) which obviously 
has t he same dimensions, but a differing dependence on the initial energy. Gea-Banacloche 
()149f ) examined the explicit time-dependence of Gaussian wave packets for the 'quantum 
bouncer' and noted that packets returned to their classical periodicity, but half a period out 
of phase, at a time half that of Eqn. ()219|) . which he described as the revival time and this 
would correspond to & p/q = 1/2 revival in the la ngua ge of Averbukh and Perelman (j3^: 
discussions of the structure of fractional revivals ()148f ) in this system have also appeared. 
The pro bability density at t = and near the first two (complete) revivals are shown in 
Fig. IHZI showing the revival structure becoming increasingly approximate as the effect of 
higher order {{n — n^xfj"^ super) terms becomes important. 



As noted in Sec. III. C| during the collapsed phase, the probability density will be consistent 
with an incoherent sum of the individual probability densities, as in Eqn. ()31|) . We observe 
this behavior at the bottom of Fig. l37l where we plot |'?/'(x, T*)p at a time T* which is not 
close to any obvious fractional revival time. The anti-correlation between the 'wiggliness' of 
the wavefunction and its magnitude, familiar from semi-classical discussions of stationary 
states or from the WKB approximation, is apparent. For additional comparison, we plot a 
purely classical probability density, given by 

PAz) = - , ^ , (220) 



based on simple ^^How much time, dt, does the particle spend in a given position bin, dx f 
arguments or, equivalently, given by the pre-factor of a WKB wavefunction. In this case, 
the classical distribution is given by 

Pc/(^;^o) = -^^= (221) 

and is shown in Fig. IHTT b) as the dotted curve. The visualization of the time-dependent 
expectation values in both t he s hort-term semiclassical and the long-term revival phases 
have been discusse d in Ref. (115 If ), where, for example, the classical position-space average 
values, using Eqn. ()221|1 . are given by 

(^) = ^ , (^2^ = M and A;2 = ^ ^0.3^0 (222) 
3 15 V45 



and quantum wave packets collapse to near these values during much of the time between 
revivals. 
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G. 2D rotor and related systems 



1. Two dimensional free quantum rotor 

The problem of a particle confined to a circle of radius L, but otherwise free to rotate, 
defines the free quantum rotor, and is described by the Schrodinger equation 

(223) 



2/ de'^ 



where I = /iL^ is the moment of inertia. (We will occasionally denote the particle mass by /i 
to avoid confusion with standard notation for angular quantum numbers when appropriate.) 
The quantized energies and normalized wavefunctions are then given by 

f-2 2 1 

Em = ^ and ^„(^) = -=6^"^^ (224) 
2i v2vr 

for m = 0,±1,±2,.... The solutions with +m, — m for |m| ^ are doubly degenerate, 
corresponding to the equivalence of clockwise and counter clockwise motion. The angular 
wavefunctions can also be written in the form 



1 / V 27r for m = 

cos(m0)/v/? for m > (225) 
sin(m6')/y7r for m > 

which displays the pattern of degeneracies in a different way. 

The classical system would have an energy, classical frequency, and period associated by 

E = \luj^ and Td = — = 27iJ^ . (226) 
2 U \ 2hi 

The corresponding quantum mechanical period requires 



so that 



, ^, , fv-m ^ / 2Em 

\E'J = -^ = N^ (227) 



Tci = ?^^= 2n/^/2EjI (228) 



as expected. 

For this system, with purely quadratic dependence on the single quantum number, the 
revival time is given by 

_ 4vr/ _ 47r/iL2 

-'-rev Z Z • yAZt)J 

h n 

This system can, of course, be thought of as a free particle with periodic boundary condi- 
tions, to be compared to ID infinite well (with reflecting boundary conditions) which has 
a very similar energy spectrum. Localized (angular) wave packets can then be constructed 
with many of the same properties as those in the ID infinite well, including zero (angular) 
momentum states with shorter revival times, and fractional revivals at p/q multiples of T^ev 
The related problem of quantum diffusion (159) on a circular one- dimensional lattice has 
been found to have similar revival-like behavior, while a propagator approach (jOSh . (99) for 
the ID problem has also shown explicit evidence for fractional revivals. The time evolution 
of the quantum rotor has also been examined in an interesting way in Ref. (,64) . 
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2. Quantum pendulum 



A related problem, with a far richer structure of classical time-dependence and quantum 
energy eigenvalue spectra, is the quantum pendulum, defined by the potential 

\/(^) = -^cos(e) (230) 

where Vq = mgl for a pendulum under the influence of gravity, or Vq = qEL for a point c harge 
q acted on by a constant electric field, E. This problem was first studied by Condon (fl6Q) 
in the early days of quantum theory, and is a staple of the pedagogical literature, appearing 
as an example of perturbation theory in many textbooks and collections of problems. The 
corresponding Schrodinger equation can be written in the form 

^2 d^^ie) ^ ^^^(^) _ (231) 



21 



This problem is of interest since it reduces to the free-rotor case in the high energy limit 
{En » +Vo), while it has the harmonic oscillator case in the limit of small oscillations, with 
predictable anharmonic corrections derivable from the cos(^) term. The classical problem is 
also more interesting as the periodicity can be evaluated in both limits, as well as in general 
case (using elliptic integrals), with the limiting case of E ^ +Vo defining the separatrix 
where the classical period diverges, namely the 'stuck on top' point. 

The quantum mechanical problem can be written as one of the familiar equations of 
mathematical physics, namely Mathieu's equation. Recently, this problem has been studied 
both in the context of understanding the connections between the classical periodicity and 
the qua ntum energy eigenstates, but also in terms of the revival behavior. The authors 
of Ref. ()161l ) use fourth-order perturbation theory and numerical evaluation of the energy 
eigenvalues to discuss the classical periodicity in the low-energy, high-energy, and separatrix 
limits, where the classical divergence in is 'softened' by quantum effects, while also dis- 
cussing the energy dependence of the quantum revival and super-revival times. For example, 
the lowest-order anharmonic correction to the energies in the low-energy limit is given by 

E^^^ = ^i2n' + 2n+l) (232) 

which gives a (non-infinite) revival time 

_ '27Th _ 327r/ 

^^^^ " jE^ ~ ~r 

which is independent of Vp a rid 8 times larger than the revival time for the high-energy, 
pure-rotor limit in Eqn. ()229|) . 

Revivals in 3D rotational systems (with energy eigenvalues proportional to Z(/ + 1)) have 
also been studied (j38|), (flfil . (IT6l . 



IV. TWO- AND THREE-DIMENSIONAL QUANTUM SYSTEMS 

A number of integrable two-dimensional infinite well or quantum billiard geometries lend 
themselves to the study of quantum revival behavior in systems with several quantum num- 
bers, and we focus here on three polygonal billiard footprints, namely the square (A^ = 4), 
equilateral triangle {N = 3), and circular (A^ oo) infinite wells. Discussions of time- 
dependent wave packet solutions of the first and third cases go back to at least de Broglie 
(0) and also provide useful examples of the connections of the short-term time-development 
of wavepackets to classical periodic orbits. 

We note that evidence of wave packet revivals has been pr esent ed from numerical simu- 
lations of quantum systems with classically chaotic behavior ()165D . with specific examples 
being illustrated for the stadium billiard. 
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A. Two dimensional infinite well and variations 

For the two-dimensional infinite square well (with dimensions L), the 

problem simplifies to two copies of a single 1-D infinite well because of the separability of the 
potential. For example, the energy eigenvalues, E{nx,ny), and position-space eigenstates, 

W{n,,ny){x,y), are given by 

/i^7r^(n^ + n^) 

E{n^,ny) = ^ and W(„^,„^)(x, = M(„^)(x)u(„^)(?/) (234) 

where the appropriate quantum numbers and the u„(x) are given by 

Eqn. (11411) . The two non- vanishing revival times are given by Eqn. (j63p and are simply 
related to each other via 

Am T ^ 

rev rev y^^^ ) 

with no cross-term (Trev'^^^) present. Therefore, the quantum revival structure is very simply 
related to that of the 1-D infinite well, including the possibilities of special 'symmetric' 
revivals for zero-momentum (po = (0,0)) wave packets at particular locations, such as for 
initial values of (xo,?/o) = {L/2^L/2) or (L/3,2L/3), as well as the same rich structure of 
fractional revivals. For rectangular infinite wells with commensurate [L^ x Ly, L^/ Ly = p/q) 
or incommensurate [L^/Ly ^ p/q) sides, the structure of the revival times may be more 
complex and interesting examples have been given in Refs. (jlll) and ("SO). 

The pattern of closed or periodic orbits is also interesting as the path lengths for closed 
orbits in the 2-D square billiard can be readily deduced from simple geometric arguments 
and are given by 

L{p,q) = 2L^/p^ + q^ (236) 

where 2p, 2q count the number of 'hits' on the horizontal and vertical walls respectively, 
before returning to the starting point in phase space. The corresponding classical periods 
for such closed trajectories are given by 

y(po) ^ ^237) 

where Vq is the classical speed. Such orbits can be produced by point particles in the 2-D 
billiard, starting from any initial location^ {xq, yo), inside the box, provided they are 'pointed' 
appropriately, namely in the tan(^) = q/p direction. The values of 



9 = tan-\q/p) and Tj°Vr = ^/p^ + q^ (238) 

(where r = 2L/vq is the period for the simplest, 'back-and-forth' closed trajectory) for many 
of the low-lying cases are tabulated in Ref. (^97|) . 

The condition for periodic orbits in Eqn. (j65|) can be implemented very easily in this case, 
and we will examine this example in detail. For such closed orbits to occur we require that 

or Uy = n^lq/p). Substituting this back into Eqn. ()234|) . as well as equating the quantized 
total energy, E{nx,ny), with the classical kinetic energy, gives 

1 



-mvQ < > E{nx,ny) 




(240) 
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or 

n = ( P and n = ( ^ (241) 

so that 

ijr'=pijr'=pf|^) = ^^^^^ (242) 

\h7rn^J vo 

which is consistent with the purely classical, and geometrical, result from Eqn. (j237p . 

We illustrate the short-term time- dependence of such wave packets, through plots of 
|A(t)P = \Ax(t)Ay(t)\'^ versus t, in Fig. EHl The results shown there are for wave pack- 
ets characterized by initial positions (xc^/o) = (-^/2,-^/2) and initial momenta given by 
{Pox,Poy) = {Pocos{9),posm{9)) where we use po = AOOtc and vary 9 = tan~^(poi;/PoxO! '^^ 
also use Axq = Ayo = 0.05 and the same physical parameters as in Eqn. p57|) . With 
these values, the classical period (for the simplest back-and-forth motion), r = 2L/vq, the 
spreading time, to (from Eqn. (f77j)). and the revival time, Trev (from Eqn. ()235p ). are given 
numerically by 

(2u)L 1 . , , 

r = = ^ 0.8 X 10■^ 243 

Po 4007r ' ^ ' 

to = (^^^ Axl = (0.05)2 ^ 2.5 X 10-3 (244) 

Tre. = ^ = -^ 0.64. (245) 

The wave packet is seen to exhibit a reasonable number of classical periods before significant 
spreading occurs, with the revival time scale being much larger than both; the locations 
of classical closed or periodic orbits are denoted by stars, with the closed orbit patterns 
indicated for some of the simpler 2D trajectories. 



B. Isosceles (45° - 45° - 90°) triangle billiard 

The energy eigenvalues and wavefunctions for a special 2-D triangular billiard 'footprint' 
can be easily deri ved (166) - (168) from those of the two-dimensional square infinite well 
solutions in Eqn. ()234|) . The standard results for the 2-D square well 



E{n^,ny) = ^—^ {nl + nj) and W(^n^^ny){x,y) = n(„^)(x)n(„^)(?/) (246) 

hold for any integral n^jUy > 1 and for rix = Uy there is a single state, while for rix 7^ n^, 
there is a two-fold degeneracy. Linear combinations of these solutions can be written in form 

= ^ -"m(^)«w(?/)] {m^n) (247) 

^(n,L)(3^'2/) = ^ + ^M(a;)^i(n)(2/)] {m^n) (248) 

^{n!n)(^'l/) = U(n){x)U(n){y) (249) 

which have the same energy degeneracy, but exhibit different patterns of nodal lines. These 
alternative forms are useful since they allow one to discuss the energy eigenvalues and 
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eigenfunctions of the 45° — 45° — 90° triangle billiard formed by 'folding' the square along a 

(n,m) ^ 

{n,m) * 

The allowed eigenvalues for this case are still given by 



diagonal, since the ^) ^'^^ satisfy the appropriate boundary condition along the new 

hypotenuse as one can easily see from Eqn. ()247|) that L,\(a;, y = x) = 0, by construction. 



E(^-^^y) = ^TZrF^(^l + ^l) E{n,n) = —-{n^ + n'') (250) 



but now with only a single (n^,, Uy) {n, n) state allowed, with corresponding wavefunctions 

given by Eqn. ()247|) . but multiplied by to account for the different normalization needed 
in the smaller area billiard. 

Th e sin gle common revival time in the 45° — 45° — 90° triangle billiard is then still given by 
Eqn. ()2H5j) and localized wave packets can also be constructed, now using the appropriately 
normalized analogs of the wavefunctions in Eqn. ()247|1 . once again, provided they are kept 
away from any of the infinite wall boundaries. 

The structure of the classical closed or periodic orbits in this case is the same as for 
the square billiard, since all of the standard {p, q) orbits in the 2-D square well are simply 
reflected off the ne w di agonal wall (along the hypotenuse), giving rise to the same allowed 
orbits as in Eqn. (j238p . The only new feature in the semi-classical propagation of such 
wave pac ket solutions in this 'folded' geometry is the existence of a special isolated closed 
orbit 13), (EM) at 135° (one which bisects the 90° right angle, bouncing normally off the 
hypotenuse, and returning to the corner) which has path lengths which are multiples of 

(Vl^ + 1^)^/2 = \/2L/2, namely half that of the standard (p, g) = (1,1) features. When 
we construct wave packets using parameters appropriate to this geometry, we see twice as 
many features in the A{t) plot for the (1, 1) case because of this special classical closed orbit, 
but otherwise reproduce the results shown in Fig. EH Extensions to 3D cubic or rectangular 
parallelepiped billiard systems are also possible. 



C. Equilateral triangle billiard 



1. Energy eigenvalues and eigenfunctions 

It is perhaps an under-appreciated fact that the energy eigenvalues and position-space 
wave functions for a particle in an equilateral (60° — 60° — 60°) triangular infinite well 
(or billiard) of side L are available in closed form. They have been discussed by many 
au thor s, in a variety of different contexts, using complementary met hods of derivation (163) 
- ()172l ) and more recently in the context of wave packet revivals ()173[ ). For definiteness 
in what follows, we will assume such a triangular billiard with vertices located at (0,0), 

(L/2, v^L/2), and (— L/2, v^L/2) and we will denote the particle mass by /x to avoid 
confusion with standardly used notation for various quantum numbers. With this notation, 
the resulting energy spectrum is given by 



E{m, n) = 




(251) 



for integral values o f m, n, with the restriction that m > 2n. (In what follows, we will use 
the notation of Ref. (11711 ) except for a minor relabeling, for the energies and wavefunctions.) 
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For the case of m > 2n, there are two degenerate states with different symmetry properties 
which can be written in the forms 
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One can confirm by direct differenti ation that they satisfy the 2-D Schrodinger equation, 
with the energy eigenvalues in Eqn. ()25H1 . as well as the appropriate boundary conditions. 
Extending earlier results, we have here also included the correct normalizations, since we 
are, of course, interesting in expanding Gaussian wave packets in such eigenstates. 

For the special case of m = 2n, there is a single non-degenerate state for each value of n, 
given by 
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Clearly these states satisfy 
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and the w^m^2n n)(^ ^ y) states also satisfy 
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0. 



The pattern of energy level degeneracies, and wavefunction symmetries is thus very similar 
to that for the 2-D square b illiar d, esp eciall y when the solut ions for that problem are written 

and (12491) . (See Ref. ()174| ) for a discussion of such pairwise 



in the form in Eqns. 



degeneracies in 2D polygonal billiards.) 



2. Classical periodicity and revival times 

Turning now to the time-dependence of wave packets in this geometry, the lon g-te rm, 
revival time scales in this two quantum number system are given by Eqns. (jp^ and ()64|) as 

\d^E/dm^\/2' '■^^ \d^E/dn^\/2' ""'^ \d^E/dmdn\ ^ ' 
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which all give the single common revival time 



rpirn) T^(™) rplrn,n) _ rp 



9/iL2 



(258) 



and exact revivals (just as for the 2-D square well) and fractional revivals are present in this 
system, with a single revival time guaranteed for any and all possible initial wave packets. 
Thus, both the = 3 and N = 4 polygonal billiards (the equilateral triangle and square) 
exhibit similar and simple energy eigenvalues and exact quantum revivals. Just as with the 
ID and 2D infinite wells, we note that special, shorter-time scale revivals are also possible 
in the equilateral triangle case for zero-momentum {{pox,Poy) = (0, 0)) wave packets initially 
centered at 'symmetric' locations wit hin the triangular billiard, such as at the geometric 
center, and half-way down a bisector ()173f l. 

The short-term, classical periodicity of quantum wave packets in this geometry can also 
be determined from calculations of 
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The condition leading to closed orbits can then be written as 
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If we substitute this condition into the energy spectrum in Eqn. ()25H1 . as well as equating the 
quantum energies with the classical kinetic energy, [j,Vq/2 (where Vq is the classical speed), 
we are led to the association 
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This implies that 
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The period for classical, closed/periodic orbits is then given by 
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where 



L{p, q) = d{p, q) = LV3\/p'^ + pq + q 



(265) 



(266) 



are the corresponding path lengths for periodic orbits. The possible classical closed or 
periodic orbits can be derived from a geometric construction (involving tiling of the 2-D 
plane, just as for the square case) giving this same result. (See Refs. f)l74) and f(175| ) 
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for discussions of the derivation of energy eigenvalues and eigenfunctions for convex plane 
polygonal billiards using tiling methods.) 

We note that the special cases of 2m = n and 2n = m correspond to g = and p = 
respectively ^ bot h of which give L{p, q) = v^L. In these cases, one of the classical periods 
from Eqn. (j260|l formally goes to infinity and this can be understood classically from the 
corresponding path length, which corresponds to periodic, 'back and forth' motion from one 
corner, along a bisector, t o the opposite side, but with no repetition in the complementary 
direction. (See also Ref. ()171f ) for a derivation of the quantized energies from which this 
effect can be also inferred.) 

Be cause of the (relatively) simple form of the allowed wavefunctions in Eqns. ()252j) and 
()253j) . involving trigonometric functions, just as for the 1-D infinite well, one can evaluate 
the expansion coefficients for any 2-D Gaussian wave packet by extending the region of 
integration from the (finite) triangular region to the entire 2-D space, so long as the initial 
wave packet is well away from any of the walls. The required Gaussian-type integrals are of 
the forms 



g^Po(x-xo)/?ie-(x--o)V26^ cos ( — 1 dx 
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with similar expressions for y-integrations. The resulting closed form expressions for the 
expansion coefficients can be used in calculations of the auto-correlation function to illustrate 
the long-term revival structure of wave packets, as well as the s hort-term, semi-classical 
propagation giving rise to closed orbits of the form in in Eqn. ()2fj5|l : for ex ampl e, the analog 
of Fig. EHl for the equilateral triangle billiard have been presented in Ref. ()173l ) . 

A 'folding' of the equilateral (60° — 60° — 60°) triangle along a bisector yields another 
special triangular geometry, namely a 30° — 60° — 90° right triangle. The energy eig enval ues 
and eigenfunctions for this case can also be trivially obtained from those of Eqn. ()252|) as 
they satisfy the new boundary condition along the 'fold'. Such a system will have the same 
energy eigenvalues as in Eqn. 1)2511) (but with only one (m, n) com bination, with m > 2n, 



possible) and the same common revival time, Trev, as in Eqn. ()258)1 . Wave packet solutions 
can also be constructed (remembering to include an additional factor of -\/2 to account for 
the normalization difference) from the {x, y) in Eqn. 



Finally, we note that Liboff (jl 
eigenstates for the hexagonal {N 
revival behavior. 



has displayed a subset of the wavefunctions and energy 
6 regular polygon) which can be seen to also give exact 



D. Circular billiard and variations 



1. Energy eigenvalues and eigenfunctions 

We finally turn our attention to the iV — oo limit of the A /"-sid ed regular polygonal 
billiard, namely the circular infinite well, which has been explored ()177[ ) in terms of quantum 
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revivals. The potential for this problem can be defined by 



, for r < R 

Vc{r) = { (269) 
oo for r > i? 

and we once again use the notation /i for the particle mass. The (unnormalized) solutions 
of the corresponding 2-D Schrodinger equation are given by 

U(^m){r,9) = J\^\{kr)e'^' (270) 

where the quantized angular momentum values are given by = mh for m = 0, ±1, ±2, ... 
and the J\m\{kr) are the (regular) Bessel functions of order |m|. The angular eigenstates can 
also be written in the form in Eqn. 



The wavenumber, k, is related to the energy via k = ^j2^E/h? and the energy eigenvalues 
are quantized by application of the boundary conditions at the infinite wall aXr = R, namely 
J\m\{z = kR) = 0. The quantized energies are then given by 

TTi ^ [^(m,nr)] ^r,^-, \ 

Em = 2/.i?2 (271) 

where ^(m,nr) denotes the zeros of the Bessel function of order |m|, and counts the number 
of radial nodes. A general time-dependent state in this system can be written in the form 

oo +00 . i.^g X 

^{r,e,t) = Y, Yl <^(m,nr)[Nim,n.)J\„,\{k^m,nr)r)] ( -= j e-'^t™-)*/'^ (272) 
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where the radial normalization factors are determined by insisting that 

R ^ 

[N(m,nr)J\m\{k(rn,nr)r)] T dv = I . (273) 



The energy spectrum is doubly degenerate for \m\ ^0 corresponding to the equivalence 
of clockwise and counter-clockwise (m > and m < O) motion. We therefore see a pattern 
of degeneracies very similar to that of both the square and equilateral triangle wells, with 
two equal energy states for each (|m| > 0,n.,.) value, and a single one for each (m = 0,71,^). 
Because the quantum number dependence of the energy eigenvalues is the determining factor 
in the structure of wave packet revivals, we need to examine the m, dependence of the 

As a first approximation, one can look at the large z behavior of the Bessel function 
solutions for fixed values of \m\, namely 

J\m\i^) ^ yjcos (^z - ^ - ^) + ■ • • . (274) 

With this approximation, the zeros are given by 

\m\ Tc TT ( Iml 3\ 

2(m,n.) - — - ^ ~ (2^r + 1)- Or ^^(m,™.) ~ ( + — + - 1 TT . (275) 

If this result were exact, the quantized energies would depend on two quantum numbers, 
in at most a quadratic manner, and there would be exact wave packet revivals, just as for 
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the 2-D square or equilateral triangle billiards. However, there are important corrections to 
this first-order formula which means that the Bessel function zeros are not given by exact 
integral values: however, a useful approximation in the large quantum number limit relevant 
for wave packet expansions can be obtained in a straightforward and accessible way by use 
of the WKB approximation. 

The two-dimensional Schrodinger equation for the radial wavefunction can be written in 
the form 

h"^ [(pR{r) ldR{r)\ h^m'^ 



and to recast this into a one-dimensional equation suitable for the WKB approximation we 
write u{r) = R{r)/^ to obtain 

^ ^ H ^ :-^u{r) = Eu{r) . (277) 



2yU dr'^ 2yur^ 

In order to obtain WKB wavefunctions with the correct behavior (same phase for la rg e r as 
the exact solutions), one makes use of the Langer transformation (|5^. ((52), ()178f l - ()181[ ) 
which effectively replaces — 1/ 4 wit h m"^- This substitution is valid for all but s-states, 
which must be treated differently ()182f ). ()185[ ). 

In this approach, we note that in the radial direction the particle moves freely up to 
the infinite wall at r = i?, but is subject to an effective centrifugal potential given by 
Veff{r) = L^/2/ir^ = (mh) /2fir'^. The classical particle cannot penetrate this centrifugal 
barrier and has an associated inner radius or distance of closest approach, Rmin, given by 

m'^ff' \m\h 

^eff{Rmin) = 7, — ^2 ~ ^ Rmin = „ ■ (278) 

We can also write this in the useful form 

Rmin = where E = — — (279) 

z 2fiR^ 

more directly in terms of the desired dimensionless variable z, which is equivalent to the 
energy eigenvalue. 

The WKB quantization condition on the radial variable, r, is then given by 



r K{r)dr={nr + CL + Cii)7T where K{r) = J^^l-^. (280) 

The matching coefficients are given by = 1/4 and = 1/2 which are appropriate for 
'linear' boundaries (at the inner centrifugal barrier) and 'hard or 'infinite wall' boundaries 
(such as at r = i?), respectively. The WKB energy quantization condition for the quantized 
energies, in terms of rir explicitly and |m| implicitly, via both the E and Rmin terms, can 
then be written in the form 



^l-^rfr = (n, + 3/4)7r. (281) 

(This result was obtained earlier in Ref. ()183f ). We also note that an improved approxi matio n 
for m 7^ eigenvalues has also been obtained through the use of periodic orbit theory ()184[ ).) 
The integral on the left can be evaluated in the form 

~ dr = Jr' - RLr, - Rm^n sec-' (-^) (282) 
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where we define x = Rr, 



IR 



\m\ z. 



(i.e., Rmin/R « 1 or \m\/z « 1) to obtain 
x"^ — X sec~^(l/x) = 1 



This result can be expanded for small values of x 
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The WKB quantization condition in Eqn. 
form 
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can then be written, in terms of z, in the 
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If we keep only the first two terms on the left-hand side, we find that 

z ^ {rij. + \m\/2 + 3/4)7r = ZQ{m, rir) 



(284) 



(285) 



which is the lowest-order result obtained directly from the limiting form of the wavefunction. 

To improve on this result, we simply keep successively higher order terms, solving itera- 
tively at each level of approximation using a lower-order result for z, and we find the much 
improved approximation 
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which we have confirmed numerically is an increasingly good approximation, especially for 
fir » 1. For the study of wave packet revivals, we only require the energy eigenvalue 
dependence on m, Ur to second order, but higher order terms such as those above might be 
useful for super-revivals and even longer-term time- dependence studies. 

For the special case of m = 0, this W KB t echnique only returns the zero-th order result, 
z = ziOjUj.), and the authors of Ref. ()177f l. motivated by the form of the expansion in 
Eqn. ()286j) . numerically fit the first 50 lowest-lying m = zeros to a similar form and find 
the result 

'io,nr) = Zo{0,n^r) + ^^^ + --- . (287) 

which gives a much improved approximation. Friedrich and Trost ^185f) have obtained an 
improved WKB approximation for the 2D circular well for the m = case which, when 
expanded in terms of -^(o,nr), gives the first two terms in Eqn. ()287|) . so we will use that 
next-to-leading result in w hat f ollows. 

Using Eqns. ()286|1 and ()287|1 . we can evaluate the energy eigenvalues to quadratic order 



in nr,m in order to probe the revival structure of wave packets. For the special case of 



m = 0, we find that 
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while for the more general case with m 7^ 0, we find 
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The fact that these energies depend on non-integral values of the effective quantum 
numbers is remin iscen t of th e case of Rydberg wave packets in alkali-metal atoms due to 
quantum defects ()186f l. ()187l ) and methods similar to those used there might prove useful. 
In what follows, however, we simply examine the time- dependence of typical m = and 
m 7^ wave packets directly. 
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2. Wave packets and time scales 



We begin by focusing on the special case of zero-momentum wave packets, centered at 
the origin, namely with vanishing values of {pox,Poy) and (xo,?/o)- For such states, where 
only the m = eigenstates contribute, we can write the energy eigenvalues from Eqn. 
in the form 
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where = nr(2nr + 3) is an integer (neither even nor odd in general). The last term 

in the square brackets is independent of Ur and will make the same, constant, overall phase 
contribution to the autocorrelation function, so we focus on the /(n^) term. Since this integer 

has no special evenness/oddness properties, its contribution to the phase of each |a(n,nr)P 
term in Eqn. (j272|) will be identically unity at a revival time given by 
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Thus, at integral multiples of 4To, we expect nearly perfect revivals because of the almost 
regularly spaced structure of the m = Bessel function zeros. At the se rec urrences, we can 
also predict the overall phase corresponding to the last term in Eqn. ()290|1 . namely 



-j?i27rV4Mfl^ (4To)(9/8+l/27r2 ) 
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where F = 1/4 + l/vr^ ^ 0.351. 

To investigate these predictions numerically, we ha ve used an initial Gaussian wave packet 
consisting of a product of two forms as in Eqn. p()7|l with the specific values 



2m, h, R=l 



and 



10^2 



so that Axq = Ayo 
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(293) 



The expansion coefficients a(^m,nr) iu Eqn. ()272|1 are calculated by numerical evaluation of 
the required 'overlap' integrals of the initial state with the eigenstates. 

Using the expansion coefficients for this state, we plot the modulus squared of the autocor- 
relation function, |74(t)|2, in the bottom plots of both Figs. IH^ and EUl with time 'measured' 
in units of Tq. The almost exact revival structure at integral multiples of 4To is evident. 
As a further check, we can evaluate the phase of Ajt) a t eac h revival and find that to an 
excellent approximation it is given by —nFn, as in Eqn. ()292j) . 

We next move away from the special case of the zero-momentum, central wave packet by 
considering individually the case of Xq 7^ and poy 7^ (but not both). In each case, the 
average angular momentum of the state is still vanishing, but m 7^ values of the expansion 
coefficients are now required. We must now use t he m ore general case for the energies, which 
to second order in m ^ 0, rir, are given by Eqn. ()289j) as 
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where 



l[nr) = n^.{2n,,. + 3 + 2\m\) 



and 



l{\m\) = \m\{\m\ +3)/2 



(295) 



are both integers, again, with no special even or oddness properties. We can then write 
these energies in the form 
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At integral multiples of the m = revival time, tj^ = N{ATq), the first two terms give 
g-Af(27ri) _ -[^ phases to each (m, rir) term in the autocorrelation function, while the last term 
gives an overall, (m nj.)-independent phase, just as in the m = case. The other term, 
however, gives a contribution 

g-(27ri)(m2Ar)(2/7r2) ^^297) 

which depends on m explicitly and which therefore eliminates the revivals, increasingly so, 
as the wave packet is dominated by m 7^ terms. However, because of a seeming numerical 

accident, at integral multiples of 5TrS~°'' = 20To, we recover approximate revivals due to 
the fact that 5 x (2/7r^) = 1.013. We thus find approximate revivals for the more general 

m ^ case given by Tje'S^") = (7rV2)T4™=°^ ^ 5T4^=^°^^. 

This effect is illustrated in more detail in Figs. IH^ and EUl where we plot |A(t)P versus t 
as we move from the central, zero-momentum wave packet by first moving away from the 
origin {xq 7^ 0, in Fig. IH^ or having non-zero momentum values [poy 7^ 0, in Fig. EUl) In 
each case, as we increase the parameter [xq or poy), we necessarily include more and more 
|m| 7^ eigenstates. For even a small mix of such states, the Trev~'^^ revival periods at 

most integral multiples of ATq disappear, while evidence for the more general Tr™"^°'' = 20ro 
revivals remains evident. 

We note that this 'lifting' of a seemingly 'accidental' degeneracy in the pattern of revival 
times is somewhat similar to the special case of a zero-momentum Gaussian wave packet in 
a 2D square or equilateral triangular billiard, initially placed at the center. 

This pattern of revival times depending on two distinct quantum numbers is also some- 
what reminiscent of that encountered in a rectangular billiard with differing sides of length 
Lx, Ly where if the sides are incommensurate one would expect a less elaborate revival 
structure. Since the revival times typically scale L , the appearance of a vr scale 

factor which can give rise to very close to an integer ratio lO/vr^ ^ 1 (to within 1.3%) 
is appropriate; in this case, the relevant length scales for the radial quantum number and 
azimuthal quantum numbers are most likely multiples of R and 27iR respectively, so that 
relative factors of vr^ in the revival times can appear naturally. 

The presence of the Am 7^ revivals becomes increasingly less obvious as the average 
angular momentum is increased away from zero (with both Xq and pQy now non- vanishing) , 
since the required energy eigenvalues are in a region of large \m,\/z where the lowest-order 
approximation (from Eqn. f)285|) ) of evenly spaced z values becomes worse. We also note 

that even with (L) = 0, as we increase poy, the spread in m values required also increases, so 
that the overall number of states required to reproduce the initial Gaussian, and which have 
to 'beat' against each other appropriately, increases as well, making revivals more difficult 
to produce. 



3. Variations on the circular billiard 



The discussion of the circular billiard can be easily extended to the case of the half- 
circular footprint, with the addition of an infinite wall along any diameter. The resulting 
energy eigenvalue spectrum is obtained in a similar way as in the cases of the square and 
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equilateral triangle billiards when 'cut' along a diagonal. The angu lar wavefunctions for 
the half-circle problem can be obtained from the form in Eqn. ()225j) by choosing only the 
sin(m^) (m > 1) forms which vanish on the new boundary. Thus, the energy eigenvalue 
spectrum consists o f one c opy of the m ^ ener gies o f the full well, very similar to the cases 
encountered in Sec. llV.li^ and at the end of Sec. llV.Cl However, because purely m = wave 



packets are no longer possible, the existence of recognizable quantum revivals is less obvious 
While we have focused on the long-term revival structure of quantum wave packets in 
the circular well, the short-term, semi-classical propagation leadin g to cl osed orbits in this 
geometry can also be studied using the same methods as in Sees. llV.Xl and lIV.Cvl and we 
present such an analysis in Appendix [Ej The 2D annular billiard, where an additional 
infinite wall at r = i?j„ < i? is added, can be studied with the same WKB methods used 
here, and we discuss that case in Appendix |El as well. 

The problem of the spherical billiard, with an otherwise free particle confined to a circular 
region of radius R, is an obvious extension, with the (unnormalized) solutions in spherical 
coordinates given by 

^(r,^,0) = jKA;(z,„0^)^m(^''/') (298) 

and the quantized energies determined by the zeros of the regular spherical Bessel functions 
via 

h"^ 2 

E{i,nr) = h/,n.))] where jK^(/,n,.)) = 0- (299) 

For such a central potential, the energy eigenvalues do not depend on the m quantum 
number, so the corresponding classical period (T^J^^) and revival times {Tr'l!^ ,Trev'^\Tre^^) 
decouple from the problem. For the special case of spherically symmetric (Z = 0) solutions, 
one has 

Uz) oc ^ (300) 

z 

and the / = eigenvalues are exactly quadratic in the quantum number; this implies that 
initially central (fq = (0,0,0)) Gaussian wave packets with vanishing initial momentum 
(Po = (0, 0, 0)) will exhibit exact revival behavior, with no higher order time scales present, 
while all other packets, which necessarily include / 7^ components, will have less obvious 
revival behavior. 

Finally, the addition of a single infinite wall along any r adius (not diameter) of the 2D 
circular billiard (adding a 'baffle^ can be analyzed in detail (11881 ). with the result that half- 
integral values of the 2D angular momentum (m = 1/2,3/2, ...) are allowed. Then, noting 
that the 2D {Jm{.z)) and 3D {im{z)) Bessel functions are related by 

jm{z) = ^j^Jm+l/2iz) , (301) 

we see that wave packets constructed from only m = 1/2 eigenstates in this geometry will 
also exhibit purely quadratic dependence on the radial quantum number and also have exact 
revival behavior. 



V. EXPERIMENTAL REALIZATIONS OF WAVEPACKET REVIVALS 

The existence of revival and fractional revival behavior in quantum bound states, first 
found numerically in simulations of Rydberg atoms (|28h. has led to a number of experimental 
tests in atomic, molecular, and other systems. We briefly review some of the experimental 
evidence for quantum wave packet revivals, while noting that excellent reviews of wave 
packet physics ^W) - (|23h . (|37h have appeared elsewhere. In addition, we discuss other 
experimental realizations of quantum revivals found in the occupation probability of a two- 
state atom system in a quantized electromagnetic field, described by the Jaynes-Cummings 
model, and in the behavior of the macroscopic wave function of Bose-Einstein condensates. 
We also discuss revival-like behaviors which arise in various optical phenomena. 
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A. Atomic systems 



As mentioned above, early formal studies (P) - ([l6l) of t he cons truction of Coulomb wave 
packets were expanded upon by a number of authors ()189l ) - ()201l ). including suggestions for 
experimental production of localized wave packets. Since excellent reviews flToh - (|23h. (|37f l 
exist on this subject, we limit ourselves to a few comments. 

The production of localized Rydberg wave packets makes use of the fact that the large 
n energy levels in such atoms are very closely spaced; for example, the energy difference 
between successive Rydberg levels for ng = 50 (An = ±1) is approximately Ai?„ = 2 x 
10^^ el^. Short-duration (At ~ Tpump) laser pulses necessarily have a range of energies 
{AE = h/At) and can therefore simultaneously excite a number of states (experiments 
have been done where wavepackets containing of order 3 — 10 states have been produced, 
corresponding to roughly An = 1.5 — 5 in our notation.) One-electron-like atoms such as 
potassium and rubidium have ionization potentials of order 4 eV, requiring laser wavelengths 
of roughly this order to excite states to E^, = —Eq/uI < for no >> 1. As examples, several 
classic experiments jHI ) f (j3^ . (|33[)) used potassium (rubidium) atoms which have ionization 
potentials given by 4.341 el^ (4.177eV^) and produced excited Rydberg states with no = 65 
(no = 62, no = 46.5) using single photon excitation corresponding to laser wavelengths of 
2858 A (297 nm); an earlier experiment (j26h also used rubidium, but utilized two-photon 
excitation with A = 594 nm, and contained only 2 — 3 component states. 

The time-development of the wave packet can be subsequently probed with a second laser 
pulse, as a function of the time delay from the first pulse. A classical description (|30t l of the 
process involves noting that the rate of energy absorption from the laser pulse by the wave 
packet is given by R{t) = 3{t) ■ E(t) where J(t) = ev()f:) is the classical electron current and 
E(t) is the time-dependent electric field. This rate is large (small) near the inner (outer) 
Keplerian turning point, where the electron speed is big (small), so that when the electron 
is near the nucleus, energy from the 'probe' laser can be efficiently absorbed, resulting in an 
increased probability of ionization (which is the observable signal), a method first proposed 
bv Alber, Ritsch, and Zoller (flTl ). Other experiments ()33h have used different techniques 
(1181 ) to monitor the dynamic behavior of the packets. 

For systems described by the Coulomb-like spectrum in Eqn. ()71|) . such as large-n Ryd- 
berg states where Z is effectively unity, the classical period will be given by 

Tf °"')(no) = ^Jj^ = (1.52 X 10~''s){no)' ^ 20ps{no/50f. (302) 

and the corresponding revival time is then given by 

T,(£°"') = (2no/3)Tf ^ 670 (no/50)^ . (303) 

For example, the original numerical predictions of Parker and Stroud (j28h used an effective 
no = 85 which gives 

Tci = 93.5ps and Trev = 5.3ns (304) 

so that the p/q = 1/2 revival shown in Fig. |3 is indeed observed to occur at roughly 
Trevf^ = 2.7ns as seen in their simulations (where Averbukh and Perelman (34) added the 
appropriate labeling to the original data, as in Fig. HI) 

In one of the early experimental observations of fractional revivals by Yeazell and Stroud 
(l3ll ). the initial wave packet was excited with no = 72 (and including roughly 5 states, or 
An ~ 2.5) which gives 

Tci = 57 ps and Trev = 2.7 ns. (305) 

In the data from Ref. fi31i ). reproduced here in Fig. |^ that initial periodicity is clearly 
visible, and a fractional revival, with local periodicity roughly half that of T^, corresponding 
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to a p/g = 1/4 fractional revival at Trev/^ ~ 680 is also apparent (note the two closely 
spaced vertical dashed lines); similar results were obtained by Meacher et al. (jS^)- 

Experimental data on rubidium from Ref. (33), shown in Fig. El corresponding to = 
46.5 and An ^ 3 (6.5 states), gives = 15ps, and a p/q = 1/2 revival is evident, with 
the initial local classical periodicity, at = 237 ps (which is labeled there as Trev) The 

inset shows similar data from the same experiment, but obtained with no = 53.3 and almost 
10 states excited in the wavepacket expansion, where fractional revivals up to order Trev/7 
are now visible, due in part to the larger value of An used, and connections such as those in 
Eqn. (j59p . Observations of fractional revivals in two-electron atoms such as calcium (20^ 
have also been reported. 

In contrast to most studies where localized electron wave packets are excited by short 
optical pulses, it has also be en sh own possible to generate Rydberg wave packets using THz 
frequency half cycle pulses f)203l ). While the short-term Kepler orbit motion and longer- 
term revival behavior of packets produced in this way are comparable, the THz wavepacket 
is initially delocalized and only becomes localized after half a revival time. 

States where the effects of the inner atomic core are important are often described by 
quantum defects, where the effective principal quantum number is given by n* = n — 6{n, I), 
with an angular-momentum quantum number dependent correction characterizing the effect 
of the inner core electrons. The eff ects of qu antum defects on the structure of wave packet 
revivals has been examined in Re fs. |186f l and (jlSTf l. while its effect on the detailed recurrence 
spectrum has been documented ()33l) experimentallv. 

The ti me-e volution of localized electron packet s in an exter nal electric field (Stark wave 
packets) ()204f l has been studied experimentally ()205f l - ()208l ). The effect of the field on 
the energy spectrum is most easily determined by solving the Schrodinger equation using 
parabolic coordinates (jH3) where the principal quantum number, n, can be written in terms 
of two parabolic quantum numbers, ni,n2, and the azimuthal (magnetic) quantum number, 
m in the form 

n = ni + n,2 + \m\ + 1 . (306) 

The linear (in applied field) correction to the energy spectrum due to an external electric 
field of the form V{z) = —eFz can be written ^84i) as 

E^^\k,n) = ^nk[eFao] (307) 

where k = rii — n2 and = h'^/me'^ is the Bohr radius. For fixed n, the spacing between 
adjacent energy levels is given by AE^^'^ = 3n[eFao], due to the form in Eqn. ()306|) which 
implies Ak = 2 jumps, which is therefore similar to the even spaced levels of an oscillator 
spectrum. The classical periodicity associated with the k quantum number will then be 
given by Eqn. (jHS)), suitably modified to read, 

(fc) 27r^ 2-nh 2.6ps 

^ 2\dE/dk\ ^ 3n[eFao] " n[F / {100 V / cm)] ^ ' 



and experimental observations of up to 10 periods have been reported ()205l ) - ()208f ). Since 
one is here concerned with only the classical perio dic ity, the application of variations of 
periodic orbit theory have also proved useful ('2071) - ()209f l. The dynamics of Stark wave 
packets above the field-ionization threshold ff210i) have also been measured experimentally, 
observing that part of the electron wave function returns to the c ore, before escaping over 

the saddle point, with classical periodicities consistent with Eqn. ()308|) . 

Novel Stark wave packets consisting of an — H~ ion pair have been produced (|211[ ) 
in which the energy states are scaled from the simple hydrogen results by 

M„ 

/i ~ me — > /i ~ -— (309) 
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since it is now the nuclear motion which is relevant; the classical perio dicity in Eqn. ()308|) is 
increased by roughly 3 orders of magnitude and data for this system ()211[ ) show good, but 
not perfect, agreement with that prediction. 

The longer-term, revival structure of Stark wave packets, in the case where the time- 
dependent states contain components of both varying n a nd k have been analyzed in 
Ref. ()212l ). The purely Coulomb revival time in Eqn. ()3U3|) is unchanged, but there is a 
non- vanishing cross-revival time given by 

^^^"^ "3[eFao] " ^ ^ 

and the revival (and fractional revival) structure of Stark wave packets differs in interesting 
ways from the free Rydberg case. 



B. Molecular systems 

The vibrational states of diatomic systems, described by anharmonic potentials with 
unequally spaced energy levels, constitute another physical system in which wave-packet 
excitations have been prepared. The subject has been reviewed in Ref. (|37h. where it is 
rightly pointed out that the excitation of localized packets in such molecular systems can 
be easier because of the one-dimensional nature of the vibrational degree of freedom, versus 
the 3D nature of hydrogen-like a toms . Vibr ational wav e packet m otion has been seen in a 
number of systems, including I2 (j213f ). (j214f ). Nal (|215[ ). and ICl (|21fi[ ). often with a large 
number of classical periods apparent. 

For example, the recurrence of the semi-classical periodicity of a vi brat ional wave packet 
at a longer revival time in the sodium dimer Na2 was exhibited in Ref. (j217f ). with a hierarchy 
of classical periodicity versus revival times given by T^i ~ 300 /s and T^ev ~ ^Ips. More 
recent experiments on Br2 (j35h have presented evidence for fractional revivals f of order 
p/g = 1/2 and 1/4) in such vibrational wave packets, as have experiments on I2 ()219l ). 

An especially illustrative method of visualizing the appearance of revivals (and fractional 
revivals) in such systems involves the calculation of a time-windowed Fourier transform 
power spectrum, or spectrogram, S'(u;,r); this is basically a two-dimensional map of fre- 
quency content of the packet versus time delay. The experimentally obtained time delay 
scan signal, s(t), is convoluted with a window function, g{t) to obtain 

S{u,t)= s{t) g{t - t) e^' dt (311) 
Jo 

where the window function is often chosen to be an Gaussian of the form g{t) = exp(— t^/tp); 
it is typical to plot ln(|S'(co', r)p) versus {uj,t). An example of such a plot (taken from 
Ref. (j35h ) is shown in Fig. 021 The dark regions near r = and t ^ 88 ps and / = 95 cm~^ 
correspond to the initial wave packet and the p/q = 1/2 revival, while the features along 
2/ = 190 cm~^ (i.e, at half the classical periodicity) and structures at 3/ show evidence of 
fractional revivals. 

R evivals and fractional revivals in NaK systems for different isotope-selected samples 
(l220f ) have shown evidence for fractional revivals of up to order p/q = 1/6, as well as being 
able to distinguish diff ering classical periods and revival times for different isotopes; isotope- 
selective studies of K2 fl221f ) have also appeared. A novel application of such vibrational wave 
packet revival behavior has been demonstrated in the laser separation of isotopes. Standard 
methods of isotope separation (.222) involve gaseous diffusion and centrifugation which, in 
turn, rely on the differences in isotopic masses. More modern laser separation methods 
(l223f ) make use of the isotopic shifts in various atomic or molecular lines. The excitation 
of vibrational wave packets in diatomic molecules with differing isotopes can yield quantum 
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revival behavior which depends on the detailed structure of the vibrational eigenfrequencies 
(and anharmonicities) which then determine the long-term free evolution and revival times. 
The difference in revival time can yield spat ial sep aration of the t wo sp ecies and this effect 
has been demonstrated experimentally ()224f ). ()225f ) and patented ()226[ ). 

Observations of revival behavior (up to orde r p/q = 1/8) in triatomic molecules have 
been reported p27) . while other proposals (I228[ ) to use differences in revival times in more 
complex molecules, ones with several vibrational degrees of freedom, have also appeared. 

Evidence for coherence in the time-development of rotational wave packets goes back 
to at least 1975 (229) with observations of short-duration birefringence in CS2 vapor with 
periodicity of order 40 ps. The theoretical background and many experimental realizations 
of such molecular structures arising from coherent production of rotational wave packets 
has been nicely reviewed by Felker and Zewail [230). Evidence for reviva l be ha vior in 
the r otational behav ior o f molecular wave packets has also been prese nted ()214D . ()23l| ) - 
()232f ) and reviewed ()233f ) . while more recent experiments ()234f ) . ()235f ) have made use of 
the revivals in such molecular rotational wavepackets to manipulate the form (phase and 
spectral content) of ultrashort laser pulses. The revival time is determined by the difference 
in energy eigenvalues of the rotational states making up the coherent packet. Using a 
standard notation for the rotational eigenvalues, Ej = BJ[J + 1), the relevant differences 
are given by 

AEj = Ej+N -Ej = B[2JN + N{N + 1)] = 2B[JN + N{N + 1) /2] (312) 

where two states differ by N rotational quanta. For molecules with no special symmetries, 
the JN and N{N -|- l)/2 values are integers (neither even nor odd) and the time to return 
to original state will be determined by 

= 2n or Trev = — (313) 

n ID 

with shorter revival times pos sible for molecules of specific symmetries where only certain 
rotational states are allowed (|229f ). Predictions of cross- revivals ()236[ ) due to vibration- 
rotation coupling and discussions of the wavepacket dynamics of rotational quantum states 
of Ceo (237.) have also appeared. 



C. Jaynes-Cummings model 

One of the most frequently discussed fully quantum mechanical models of the interaction 
of a two-level atom with a single mode of the quantized electromagnetic field was proposed 
by Jaynes and Cummings fi238ll . The Hamiltonian for this system is given by 

H = -^^(73 -|- hjjjaa) + hX ((T+a -|- a^o"_) (314) 

where a^, a are the raising and lowering operators for the boson field mode of frequency a;, 
and the 0"+, o"3 are the Pauli matrices representing the two-state system, with hw^ being 
the energy difference between the two le vels. Discussions of this system have routinely 
appeared in reviews of the subject ()239l ) - ()242f l and the model is known to have analogs in 
many other of physics. 

Applied to a two-level atom in a resonant cavity with n photons, the system will undergo 
Rabi oscillations with frequency given by 

= {n\^ + AV4)^/^ with A = ujq-uj (315) 



65 



which simphfies in the de-tuning hmit when uj — ^ 0. For that case, the probabihty that 
the system is in the excited state is given by 

Pe,n(t) = ^ (1 + cos(2v^At)) (316) 

while for a system with a distribution of photons, the solution is averaged over the initial 
probability distribution, p„, to give 

oo 1 1 "'^ 

Pe{t) = J2PnPe,n{t) = 3 + 2 COs(v^At) . (317) 

ra=0 m=0 

For an initial coherent state distribution, one has pn given by the Poisson distribution, 
yielding 

1 e-l"l' lap" 
Pe{t) = 2 + E -jr cos(2v/^At) (318) 

ra=0 

and the average value of n is given by n = |ap, while the spread is An = This 
expression for Pe{t) has many obvious similarities to the autocorrelation function, A{t), 
for wave packets, namely highly localized expansion coefficients (when rI=|ap>>|Q;| = 
An >> 1) and oscillatory terms which are not purely harm onic, so it is not surprising that 
some aspects of the short- and long-term behavior of Eqn. ()318|) have features in common 
with A{t). 

The dynamics of this system exhibit Rabi oscillations with a frequency centered at flr r, but 
with a dephasing given by a Gaussian envelope, exp(— (At)^/2), first derived in Refs. ()243| ) 
and (244), and later improved upon (24l) for A ^ 0. For Ion ger times, using familiar 
physical reasoning (or more formal expansions about n — n), one can see that at the time 
that neighboring terms in Eqn. (|318|) acquire a common 27r phase difference, one expects the 
Jaynes-Cummings summation to return to close to its t = behavior, exhibiting revivals, 
this time due to the quantized nature of the electromagnetic field. The condition for this to 
occur is 

'^^n+lTrev " '^^nTrev = 0n+l ^ (t>n = 27r (319) 

yielding a revival time given by 

27rv/A%+AV4 27rv^ 
Trev = — ^ — — as A = . (320) 

This time scale, given by how long it takes the phases from consecutive frequencies to catch 
up to each other, is more akin to the classical periodicity in Eqn. in the formalism we 
have used so far, but is routinely referred to as the revival time. 

The result in Eqn. ()32()l) wa s de rived in numerical simulations by Eberly, Narozhny, 
and Sanch ez- M ondragon (245), ()246f ). but one should also note the much earlier work by 
J. Frahm ()247l ) who also demonstrated very similar results in this system. Eberly et al. 
also found evidence of subsequent revivals at integral multiples of Trev, as well as accurate 
representations of the decreasing magnitudes at t = kT^ev, and discussed as we ll th e even 
longer-term approach to more irregular behavior, due to overlapping revivals ()248[ ). For 
example, they were able to demonstrate that the peak heights were bounded by the long- 
term limits (when A = 0) 

«W = 5^5(TTA¥7ii)vJ- 

As an example of their solutions, we show in Fig. 0^1 plots of Pe{t) for n = 36 and 
A = 0.01. We note the short-term Rabi oscillations, with Gaussian dephasing factor in 
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Fig. HHTa), as well as the presence of increasingly width revivals at integral multiples of 
in Fig. H3f b). (Plotted as a function of r = At/vr, the revivals predicted by Eqn. ()320|) 

are expected at r = k{2\fn) = 12k for this example.) The reappearance of the 'classical' 
periodicity near the first two revivals ((c) and (d)) is also apparent. We note that there are 
similarities in the spreading time formalism appearing in the envelope functions found by 
them with the formulae by Nauenberg ff40i) : Fleischhauer and Schl eich (249.) have, in fact, 
applied a Poisson summation approach to the evaluation of Eqn. ()H18|1 which makes this 
connection more apparent. The time-development of this s yste m in phase space, in terms 
of the quasiprobability distribution of Cahill and Glauber ()251l ) (or Q- function), has been 
presented by Eiselt and Risken ()252f l. 

The direct observation of these effects (with at least two obvious revivals) was first demon- 
strated by Rempe, Walther, and Klein ()253[ ) using Rydberg atom states interacting with 
single field modes in a su perconducting cavity constituting a one-atom maser ()254f ). More de- 
tailed observations ()255[ ) have been able to extract discrete Fourier components of the Rab i 
oscillation time- dependence proportional to the square root of integers, as in Eqn. ()318p . 
providing very direct evidence of t he qu antization of the electromagnetic field. (It should 
be noted that Wright and Meystre (!25f?) have examined in detail the collapses and revivals 
in micromaser systems and found subtle differences with those in the Jaynes-Cummings 
model.) 

Theoretical work extending these results further to fractional revivals ()257l ). ()258f ) and to 
sub-Poisson photon distributions (258.) has also appeared, as well as suggestions for using 
the revivals in the population inversion to measure novel phenomena (j259); variations on 
these results in an optical (as opposed to a microwave) cavity have also been discussed 
pBDi) . Studies have shown (261.) - (26^1) that at half the revival time it is possible to obtain 
a entangled, Schrodinger-cat-like, atom-field state, similar in some ways to the structure 
in Eqn. ()43|) and Fig. EOl bu t wi th no short-term periodicity (no Rabi oscillations), while 
extensions to A^-level atoms ()264I ) have been observed to give different behavio rs. 

Other similar models which are soluble in closed form have been presented f|265[ ). ()266f ). 
while realizations of similar phenomena in the context of laser-driven excitation of electronic 
transitions in diatomic molecules (268) have also been discussed. 

Using harmonically bound ions, it has proved possible to create non-classical motional 
states of trapped atoms, including thermal, Fock, coherent and squeezed states of motion 
(|269.) . including observations of the classical periodicity. The coupling between the internal 
and motional states is in a regime which can be described by the Jaynes-Cummings model, 
so that the evolution of the atomic state can provide information on the number distribution 
of the motional state. 



D. Revivals in other systems 



1. Atoms in optical lattices and Bose-Einstein condensates 

Optical lattices, formed by the interference of multiple laser beams, can act as local peri- 
odic potentials for atoms. Anharmonicities present in the potential well can cause dispersive 
spreading of wave packets formed in such systems, while dissipation can also arise from spon- 
taneous emission, leading to loss of coherence. Localized coherent-state-like wave packets 
can be formed by suddenly shifting the optical lattice and two groups (p7fll ). ()27lf l have 
examined the interplay between wave packet spreading, dissipation, and tunneling between 
adjacent wells and their impact on wave packet revival times. 

Motivated by the production of Bose-Einstein condensates, several groups analyzed the 
time-dependence of the macroscopic wave function (j272f ) - (j275l ) and found periodic collapses 
and revival behavior. The energy of the external potential experienced by the atoms could 
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be modeled as being due to the interaction energy between the atoms, giving an effective 
Hamiltonian of the form 

H = Mpl) (322) 

where h is the number operator for atoms in the confining potential, and Uq is determined 
by the inter-particle interactions through the s-wave scattering length, a, and the ground 
state wavefunction, w{x), to be Uq = (inh'^a/m) f |w(x)|^(ix. 

In this case, the initial coherent state excitations can be of the form 



\^ " e-^^°"("-i)*/2;.l (323) 



n=0 

where |ap = N and |a| = AA^; the time-dependence has at most a quadratic dependence on 
n. Using by now familiar techniques, we see that since n{n — l)/2 is an integer (neither even 
nor odd), there is a common revival time given by T^ev = 2Tch/Uo at which \ a) +_^_T-^__ = \a)t- 
This effect was confirmed experimentally for a Bose-Einstein condensate (|6G ) confined to 
a three-dimensional optical lattice, where the collapse and (approximate) revivals of the 
number of atoms in the coherent state was monitored. 

We note that the macroscopic matter field, ip = {at\a\a)t, can be written in the form 
(|8f 



%jj{t) = a exp < — |a 



|2 



, 27rt\ f 27rt 

1 — cos — — + i sin 



rev / \ rev 



(324) 



which is very similar in form to the autocorrelation function, A{t), in Eqn. ()125p for the 
minimum-uncertainty, coherent-state like solution of t he harmonic oscillator. The collapse 
time for this syste m can also be derived (|88h . (j272f ) - (l275[ ) from the short-term time- 
dependence of Eqn. fl324j) giving Tcou oc Trev/\ct\ = Treu/AAT 

The dynamical behavior of \a)t can also be examined in the context of fractional revival 
behavior by evaluating the time- dependence (especially the interplay between the n and ra^ 
terms in the exponentials) of each term at p/q multiples of Trev For example, at Trev/2, 
each term in the coherent state expansion contains a term of the form 



.72 ^ ^ 



using the expansion in Eqn. (021) • When used in the evaluation of \a)t one finds 

|a)i=T_/2 = ^ + e+^-Z^-^'^/^^)) (326) 

V 2 

which is another Schrodinger-cat like superposition, with two distinct states rotated by 90° 
from the initial |a) state. The structure of these coherent states are visualized in a standard 
way by defining the overlap of \a)t with an arbitrary coherent state, via 



\a)t = e-l"l'/2 e-l/5|'/2 ^ ^A^e-™(n-iWT_ (327) 



n=0 



and then plotting 

P[Re{f3),Im{(3);t] = P{f3-t) = (328) 

and several examples are shown in Fig. for example, the case of Fig. EI^c) corresponds 
to the state in Eqn. ()326|1 . Related ways of visualizing such coherence phenomena include 
the use of Wigner quasi-probability distributions (as in Ref. (3^ for fractional wave packet 
revivals or Figs 1^ and O0|) and the so-called Q function (,,63.) . (,277,) in quantum optics. 
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2. Revivals and fractional revivals in optical systems 

The quantized structure of important physical properties (energy eigenvalues for quan- 
tum wave packets, or the quantized EM field for Jaynes-Cummings systems) are not limited 
to purely quantum phenomena. The existence of full revivals, mirror revivals (with a ref- 
ormation of the original coherence, but out of phase), and fractional revivals, due to the 
discrete nature of classical wave systems, has been observed experimentally in several opti- 
cal phenomena. In most cases, the observation of revivals in the spatial distributio ns o f light 
is referred to as self-imaging (for an excellent review of this general topic, see Ref. ()278f l) and 
we will focus on two examples, the Talbot effect and self-imaging in waveguides, as they are 
the most analogous to wave packet revivals in their mathematical structure and analyses. 

In 1836 H. F. Talbot f„279i) (a co-inventor of photography, with Daguerre) illuminated 
a diffraction grating with a small (coherent) white light source and examined the resulting 
transmitted light with a magnifying glass. He noted recurring patterns of colored bands, 
repeating themselves as the lens was moved further away from the grating. Rayleigh fii28d[ l 
correctly interpreted these as resulting from the interference of the diffracted beams, in 
what would now be called the near-field regime. Subsequent work showed that the self- 
images of the original periodic structure, illuminated in plane wave approximation, with 
monochromatic light, would appear at multiples of a distance Lt = 2a^/A where a is the 
grating distance and A is the wavelength. The detailed mathematical analysis showed that 
the same pattern would also recur at ^ = Lt/2, but shifted in space by half the spatial 
periodicity (the analog of a mirror revival.) Iligher order fractional revivals, at rational 
multiples of {z = pLx/q) consist of q copies of the original grating, separated by a/q. 
In the context of self-imaging, the revivals (and mirror revivals) are oft en c alled Fourier 
images, while the fractional versions are referred to as Fresnel images ()281[ ). The clear 
physical and mathematical analogies to the structure of wave packet revivals, an d the ir 
i ntric ate dependence on number theoretic identities, has been extensively discussed ()282l ) - 
( 285f ) and the Talbot effect ()286f l. along with corresponding effects in wave packet revivals 
(64fl ()129 ). have been put forward as a novel way to factorize numbers. 

The Talbot effect is not limited to pu rely electromagnetic waves, but has also been dis- 
cussed for matter waves in atom optics ()287f ) and observed experime ntally with relatively 
light atoms ()288f ) - ()29n h as well as with large molecules such as Ceo (j291[ l. Fractional re- 
vivals, i.e., higher order Talbot fringes, of up to 7-th order have been observed (j290), which 
is coincidently quite similar to the fractional revival resolution obtained with the best wave 
packet studies (|33h. 

The formation of images by phase coincidences in optical waveguides is another example 
of such self-imaging processes. The prediction that real optical images could form in planar 
optical waveguides was made by Bryngdahl (292*), was then demonstrated in a series of 
elegant experiments by Ulrich and collaborators (293) and resulted in several patents f.294) . 
For this system, the recurrence length is given by L = AufW^/X where Wa is the effective 
aperture size in this slab geometry (including penetration depth effects) an d tzf is the i ndex 
of refraction of the dielectric slab. Multiple/fractional image formation ()293f ). ()295f ) has 
also been observed, and the phenomenon has also been demonstrated with a;-ray waveguides 
(j296'). The planar waveguide geometry corresponds to the one-dimensional infinite well, 
with W corresponding to the well size, and the length along the waveguide corresponding 
to time. The ray tracing visualization of th e intensity patterns in this geometry leading to 
fractional revivals (see especially Ref. (j293) for a nice example) are therefore quite similar 
to the semi-classical pictures discussed by Born ()101l) . (102.) for trajectories in the infinite 
square well, leading to quantum carpets. 

T he application of this effect for use in novel optical beam splitters has been advocated 
()298[ ). The generalization of t his self- imaging method to two-dimensional waveguides of 
rectangular cross-sections ()297f ) has also been demonstrat ed and fiber optic geometries of 
equilateral triangle and hexagonal cross-sections (as in Sec. IIV.(!|l . have also been discussed. 

Finally, we note that while observations of wavepacket revival behavior has been mainly 
limited to atomic and molecular systems, the possible use of localized low-energy wave packet 
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excitations (I299f l localized at the edge of a 2D electron system (edge magnetoplasmons) 
to probe the dynamical properties of the integer and quantum Hall effects has also been 
proposed. 

The study of classical and quantum mechanical systems subject to external periodic 
forces, focusing on issues related to chaotic behavior, has a rich literature. A number of 
such studies have focused on the appearance of revival- type behaviors f(152i) . f(153) . f SOOf ) - 
()303[ ) in such systems. 



VI. DISCUSSION AND CONCLUSIONS 



The connections between the energy eigenvalue spectrum of a quantum bound state 
system and the classical periodicities of the system have been a standard subject in quantum 
theory since the first discussions of the correspondence principle by Bohr. 

Some semiclassical techniques, such as periodic orbit theory, can connect the quantum 
energy spectrum with classical closed or periodic orbits, but often do so in a way which 
does not exhibit the time-dependence of quantum wave packets. Truly dynamic observa- 
tions of localized quantum wave packets, exhibiting the classical periodicities expected in 
such semi-classical limits, have become possible with observations in atomic, molecular, and 
other systems, with analog behaviors seen in atom-field, BEC, optical, and other systems 
as well. The medium-term collapse of localized quantum states also present in such systems 
is familiar from simple examples of spreading wave packets which were constructed during 
the early days of the development of quantum theory. It can be analyzed quite generally for 
such bound state systems and has been observed in a number of systems as well. 

The truly novel observation that such wave packets can relocalize and once again ex- 
hibit the classical periodicity, first observed in numerical simulations, has now been widely 
confirmed by experiments in a large number of different physical systems including atomic 
(Coulomb and Stark effect) systems using both electronic and nuclear states, and vibra- 
tional and rotational states in molecules. Besides being a fundamental realization of the 
discrete nature of quantum bound states, the simple time- dependence of such eigenstates, 
and interference effects, the phenomenon of quantum packet revivals has been increasingly 
used in the development of modern quantum control experiments, to assist in the shaping 
of specific quantum states, and is likely to remain an important aspect of the production of 
specific target states in the future. Given the relative simplicity of the quantum mechanical 
background which is needed to understand many aspects of these effects, and the connec- 
tions to a wide variety of other revival phenomena (in optics and elsewhere), it is likely that 
this important manifestation of the time-development of quantum mechanical bound states 
could easily find a place in the undergraduate and graduate curricula, and this review can 
be seen as one step towards that end. 
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APPENDIX A: Energy eigenvalues in power law potentials 



by 



The energy eigenvalue spectrum in the family of power-law potentials ()55l ) - ()58l ) given 



(Al) 



can be obtained for large n (which is the situation encountered in wave packet revivals) by 
the use of the WKB approximation. In this case one writes 



2m{Eit^ - V^k){x)) dx = {n + CL + Cr) Hit 



where 



i/fc 



L 



(A2) 



(A3) 



are the classical turning points. The matching coefficients have been discussed in detail in 
Ref. ^5% and are given by C^, Cr = 1/4 for smooth potentials, but with Cl, Cr — ^ 1/2 for 
the case of infinite walls (i.e., in the limit when k — > oo). The integrals can be done using 
standard handbook results to give 



{n + CL + C, 



hir 



R) 



--Vr 



2LV2m 

£^''\n + CL + CRf''/^'+'^ 



i/fc r(l/A; + 3/2) 
r(l/A; + l)r(3/2) 



2fc/(fc+2) 



(A4) 



which does reproduce both the oscillator {k = 2) and infinite well k oo) limits, for large n, 
and even gives the appropriate scaling for the Coulomb {k = —1) case. The same approach 
can also be used for 'half well potentials, wher e an in finite barrier at x = is introduced, 
as for the case of the 'quantum bouncer' in Sec. IIII.Fl for appropriate values of Cl, Cr. 



APPENDIX B: General time scales 



The general expansion in Eqn. (jH} defines three time scales, namely the classical peri- 
odicity, the revival, and superrevival times, all of which we have discussed in detail. In 
this Appendix, we describe two other relevant time scales, as well as another quite general 
approach to visualizing the short-term classical periodicity and subsequent spreading. 

Nauenberg (|40| ) noted that, for the Gaussian coefficients in Eqn. (j24|) . the autocorrelation 
function could be written (in the notation used here) as 



n=0 



-(n-no)2/2An2 

An^y2n 



J[E{no)t+E' (no)(n-no)t+E" {no){n-no)^ /2]/h 



to second order in [n — uq) . Using the Poisson formula. 



oo +00 p „_ 
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he found he could write A{t) in the suggestive form 

piE{no)t/h +~ 



A{t) 



2TTAnJa{t) 
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where 



a{t) 



47]-2 



1 



4:7Tit\ 



T 



(B4) 



(An)2 

and Tci and T^e^ are defined, in the notation used here, by Eqns. © and (fT^ . This form 
exhibits, in a quite general way, the periodicity of the autocorrelation function, with the 
summation of exponentials centered at integral multiples of T^. It also exhibits the spreading 
(due to dispersion) apparent in the time-dependent width given by 



A{t) 



1 



1 + 



1/2 




27rAn 



where 



T 



47r(Ara) 



(B5) 



(B6) 



Th is type of spreading is similar to that due to the pre-factor terms in Eqns. ()82|1 and 
fll02j) . which were written in terms of the familiar spreading times, to? in the free-particle 
(infinite well) or uniform acceleration (quantum bouncer) cases. I n or der to exhibit this 
connection more generally, we use the power-law potential in Eqn. ()A1|) and the results in 
Appendix^ For example, we can write the quantized energies (in the large n limit) for the 
general V(fc)(x) potential in the form 



2m 



which we also equate to the maximum value of momentum, p„ 



4:7!- h 



n 



{k + 2) 



2k{k - 2) 



(B7) 

In these systems, we have 

(B8) 



For an initial Gaussian wave packet of width Ap„, the spreading time can be written in the 
form 



(k) 



mh 



° 2Apl 

and we can use the identification in Eqn. ()B7|) to equate 

k 



Apn = v/2^^n-2/('^+2) 



k + 2 



An 



which gives 
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Thus, for this family of potentials, we have the relation 



Tr 



(k) 



Stt 



k 
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and this ratio is equal to 87r(An)^ for both k = 1 (quantum bouncer, uniform acceleration) 
and A; — >■ oo (infinite well, free particle). We then note that the spreading time, Tg, given by 
the Nauenberg formula for these cases, is calculated to be 



^(fc) 

-L rev 

An (An) 



2tn 



(B13) 
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which agrees with the results in Eqn. ()82|) and p02|) . 

Nauenberg also pointed out that the wavepacket spreads sufficiently that quantum self- 
interference occurs on a time scale of (in our notation) Trev/ marking the approach the 
collapsed state. The notation of a collapse ti me h as also been studied in the context of the 
infinite (i94j well and the quantum bouncer ()149f l. with quite different approaches. For the 
infinite well, the time required for the various ex pecta tion values to approach their semi- 
classical limits was given approximately by Eqn. ()174|) which can be written in the form 

T,,u = (B14) 



while for the quantum bouncer, a direct examination of the phase differences between various 
terms in the eigenstate expansion led to the result 

(where we use our notation for Trev, since the one described in Eqn. (42) of Ref. (|149f ) is more 
rightly associated with the p/q = 1/2 revival in the language of Averbukh and Perelman.) 
A third, quite general measure of the time taken to co llap se to the incoherent state would be 
to use the spreading (dispersive) pre-factors in Eqn. ()B3|1 . or the similar ones in Eqns. (jH^ 
and (jlU2p . and define the collapse time as how long it takes for the envelope of |A(t)P to 
decrease to | AncP, namely when 



, , ^= = \AiT^ou)\ = \Ainc\= 7= (B16) 

T.OU = (4v^An) to = 777^^ (^] (B17) 



which gives 

which gives the same general form as the explicit analyses resulting in Eqns. ()B14|) and 
flB15|) . We note that the earliest time at which a possible fractional revival which might be 
resolvable (over the incoherent background) is given by the condition 
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T 

rev 



11 T T 
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at the time Teariy = = -r — /— (B18) 



q An2^ '""^ q An20F 



so that this time scale is also singled out as one 'comes down' in time from T^,,, . as well as 
from 'going up' in time from t = 0. We note that Fleischhauer and Schleich ()249[) made 
similar use of the Poisson summation f ormu la to obtain improved approximate expressions 
for the Jaynes-Cummings sum in Eqn. ()318|) . 

Finally, the presence of the T^ev term provides the dispersion necessary for the decay to 
the collapsed state. For the harmonic oscillator system, when Tr^ ,, 00, we would expect 
to reproduce something like the result of Eqn. (jl25p . Using Eqn. i ()B3|) . with T^^.^ 00, we 
find that 

+00 

Ait) = e*'^("o)*/'* ^ g-{m-i/Td)'(27rAn)2/2 _ (B19) 
m=— 00 

Near t ^ 0, this reduces to 

A{t) = e«-^{"o)t/?ig-(An)2a;2t2/2 (B20) 

since 2n/Tci = u. The exact result in Eqn. ()125p was derived for P = Pq, but arbitrary 
xo,Po, and using Eqns. ()118|) (to evaluate {E)) and ()122|) (to evaluate An), we find that the 
oscillator result for ut << 1 can be written in the form 

A{t) = e-^mtin^-{i^nn.t?i2 (321) 

which agrees with Eqn. ()B20|) in that limit. More generally, the two formulae agree very 
well for |t — /cTc/| << T^-, i.e., near any multiple of the classical period. 
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APPENDIX C: Fractional {p/q) revivals for even q 



The derivations of Averbukh and Perelman (js^) for the temporal structure of fra ctional 
wave packet revivals at t = pTrevlo. for odd values of q was discussed in Sec. Ill.DI and we 
extend the analysis here. For case of even q values (implying p is odd, since p, q are assumed 
to be relatively prime), the periodicity in / required to satisfy Eqn. is given by 

l = q for g = 2 {mod A) (CI) 

l = q/2 forg = 4(morf4) (C2) 

and we can treat each case separately. 

For the first case, we note that the recursion relation in Eqn. (j53j) will connect 6^. values 
with even and odd r separately. Since q = 2 {mod 4), q/2 will be an odd integer and we can 
rewrite the expression for using the relabeling n = n + g/2, to write 

ho 



since 

because q/2 is an odd integer and p is necessarily odd if q is even. Since 6o vanishes, all of 
the even do as well, and only the q/2 values with odd r are non- vanishing, leading to 
q/2 'clones' or 'mini-packets' near the fractional revival time in these cases. 

An explicit example is for the half-revival time, Tj.ev/2 where p = 1, g = 2, and I = q/2 = 1 
where one can explicitly find that 

fc=0 

and 

fc=0 

and we obtain the result of Eqn. (j38|) . The next case in this series is for multiples of Trev/G, 
which turn out to be similar to the p/q = 1/3, 2/3 cases. 

Finally, for the case of g = A {mod 4), the periodicity in Eqn. is given by / = g/2 
and the analysis proceeds in a similar fashion. For example, for the case of Trev/A, we have 



i ^ g^27rip(n+g/2)Vg 
n 

1 ^-2Tripn^/q ^-2mpn ^-mpq/2 (C3) 



1 ^ g-27ripnVg 
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and e-^^^"'/^ _ Tg— ^21^ = (_i)p (C4) 
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p = 1, q = 4, I = q/2 = 2, and the various expansion coefficients in Eqn. ()48|) are 



2^ 2^ ' 72 

h = -J2 e^^'^'/^-'"/'^ = 1(1 - 1) = (C7) 



1 ^ 1 1 

b2 = i^e2-^(2/c/4-fcV4) ^ 1(1 +z) = ^e+^-/^ 







which explains the resuh of Eqn. (jlgj) in a more 'turnkey' fashion. 

We note that Bluhm and Kostelecky (ffW ) have obtained similar results for the algebra of 
the complex phases arising from the superrevival terms of the form exp{—2iTipk^ / q) . They 
have also extended this formalism of fractional revival analysis to the case of systems with 
two or more quantum numbers (ffoi). 



APPENDIX D: The 'inverted' oscillator 



Fo r the case of the so-called 'inverted' oscillator, the general wave packet soluti on in 
Eqn. (jlllj) . for example, can be directly carried over using the identifications in Eqn. 
to obtain the 'runaway' wavepacket, with probability density given by 
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{x — Xo cosh(u;t) — po sinh(u;t) / muj 



~.\2 
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where 
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Axt 
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B{t) \ = cosh" (cut) + {h/mCo(3ysmh\ut) . 
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Just as for the harmonic oscillator, the expression for A{t) for the general case is cumbersome, 
so we only examine it for one specific case as an example, namely the case where /3 = /3o = 

^yh/muJ. This situation no longer corresponds to a constant width wave packet, since 



Ax. 



cosh" (cut) + sinh"(ci;t) 



(D4) 



increases exponentially, as the individual momentum components comprising the wave 
packet quickly diverge in p-space. For the case of Xq = 0, we have the general expression 



A{t) 
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a/ cosh(cijt) 



exp 



f Pi 
V 2mujh 



cosh(ti;t) — 1 + z sinh(Lo't)(2 cosh(co't) — 1) 
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In the limit when t » l/cu, the hyperbolic functions both approach exp(co't)/2 and we have 
the limiting case 

A{t) — > -^=^==== exp 
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2mujh 



(D6) 



The exponential ('dynamical') suppression once again is seen to 'saturate', as in the free- 
particle case, and for the same reason, namely that both x{t) — Xq and Axt have the same 
large t behavior. The resulting modulus is given by 



2e-'^* exp 



mujh 



(D7) 



which still becomes exponentially small, but now due to the ('dispersive') prefactor. If one 
also has xq 7^ 0, the expression above includes an additional factor of exp(— Xq//?q) (similar 
to that in Eqn. (|HT|. with no cross-term involving Xq times Pq. 



APPENDIX E: The full and annular circular wells: WKB energy eigenvalues, 
classical periods, and closed orbits 

While we focus on the information about wave packet revival times encoded in the energy 
eigenvalue spectrum of quantum systems, it is also interesting to see how the pattern of closed 
(or periodic) orbits supported in a number of simple 2D quantum billiards arise from the 
connections between the classical periods in systems with more than two quantum numbers. 
This is especiall y tru e sinc e most of the experimentally observed data from 2D circular 
billiard systems 1)8041 ) . (JSCS') have involved measurements which are relevant for short-term, 
quasi-classical ballisti c pro pa gatio n. Such closed orbits are also the ones of relevance to 
periodic orbit theory ()306[ ) - ()308[ ) analyses of such billiard systems. 

We have illustrated these connections using the explicit expressions for E{n,m) for both 
the square and equilateral triangle billiard and we now turn to the cases of the circular 
and annular infinite wells. (T he a rguments presented here are adapted and extended from 
Ref. ()177f l. but see also Ref. ()309l ) for discussions of the circular billiard in the context of 
what has come to be known as periodic orbit theory.) For the circular well, we first note that 
the allowed closed orbits are characterized by two integers {p, q) , where q counts the number 
of 'net revolutions' the trajectories make before closing on themselves, while p counts the 
number of 'hits' on the circular walls. For consistency, one must have p > 2q and when 
p, q have common factors, say p = np and q = nq, one simply has an n-fold recurrence of a 
basic (or primitive) closed path. A number of low-lying trajectories of this type are shown 
in Ref. (.309) . The total path length for one classical period in any of these closed orbits is 
given by 

L{p,q) = 2pRsml — ] (El) 



P 



so that the classical period is simply 



T(P.,) = ^ (E2) 

where Vq is again the classical speed of the point particle inside the billiard. The minimum 
distance to the origin for any of these trajectories (distance of closest approach) is given by 

Rmin = R cos . (E3) 

For the quantum case, the two appropriate quantum numbers give rise to classical periods 
given by 

\dE/dnr\ \dE/dm\ 
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and the two periods can beat against each other to produce the classical periodicity iT^"^] 
for closed or periodic orbits if they satisfy 



T, 



(po) 



cl 



rr(m) 



(E5) 



again, with p > 2q for this geometry. 

We can then use this formalism to understand how these conditions can gi ve r ise to the 
classical expressions for the path lengths and minimum radii in Eqns. ()E1|) and ()E3|) . Instead 
of using the approximate expression in Eqn. ()289|) for the (m, nr)-dependent energies, we 
make use of the WKB condition in Eqn. 1)2811) . namely 
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{rir + 3/4)7r 



(E6) 



where the appropriate + Cr = 1/2 + 1/4 = 3/4 factor corresponds to infinite wall 
boundary conditions at r = R and 'linear' boundary conditions at the inner 'turning point'. 



Rm.ir, 
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which defines the useful parameter z. 

We then simply take partial derivatives of both sides with respect to Ur and m respectively. 
We thus obtain the conditions 
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The condition to be satisfied for periodic orbits can then be written as 
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Evaluating the integral and using Rmin = \m\h/y/2jlE, we find that 
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as the condition on periodic orbits, as expected. To find the classical period for such closed 
orbits, we note that 
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where we identify vq = >/2EJJ1 with the classical speed which reproduces the purely geo- 
metric result of Eqn. (jEljl . 
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This type of analysis can be easily generalized to the case of an annular billiard, a circular 
system with infinite walls at both the outer radius R as well as at an inner radius Rq = fR 
(where < / < 1), defined by the potential 



VA(r) 



for i?o < < -R 

oo for r < Ro and r > R 



(E13) 



This system has also been studied in the context of periodic orbit theory (131 
priate quantum wave functions are given by 



The appro- 



(E14) 



where one now includes the 'irregular' (divergent at the origin) Y\m\{kr) since the inner 
wall now guarantees that r > 0. Application of the boundary conditions at r = i? and 
r = Rq = fR yields the corresponding energy eigenvalue condition 



J\^\{kR)Y\^\{kR 
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J\m\ikR. 
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which can be solved numerically just as easily as in the purely circular case. 

In terms of the connection between the classical closed orbits an d th e quantum periods, 
we note that the same periodic orbits described by Eqns. ()E1|) and ()E3|1 are still allowed, so 
long as 
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(E16) 



as shown in Fig. |^ Another set of orbits is allowed, also characterized by the same set of 
integers (p, g), which bounce off the inner wall, an example of which is also shown in Fig. 031 
For t his g eometry, the corresponding path lengths are always at least as large as as those in 
Eqn. ()E1|) . and are given by 



L(p, q) = 2pR^l + P-2f cos(7rg/p) . 



(E17) 



The two classes of periodic orbits and path lengths coalesce in the limit that / fmax = 
cos(7rg/p) and then both disappear from the allowed set of paths. We know that the effective 
distance of closest approach is an important parameter in this geometry, and for the case of 
the annular ring, we can define Rmin = zR as before, but we now have 
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sin(7rg/p) 



^l + /2_2/cos(7rg/p) 
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as illustrated in Fig.HHl We note that 2; < / (as expected) for these trajectories which bounce 
off the inner wall, with the limiting case ('just touching') a rising wh en / = fmax = cos(7rg/p). 

The WKB analysis for this geometry proceeds as in Sec. lIV.D.ll with the lower integration 
limit changed from Rmin to i?o, and the WKB 'matching coefficient' changed to reflect the 
fact that the inner boundary is now also of the 'infinite wall type'. These changes give 
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with Rmin = zR with / < z as noted above. The implicit differentiation proceeds as before 
and the condition for closed orbits becomes 
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The integral can be done exactly and using Rmin = zR and i?o = /-R, we obtain 
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as the condition for closed orbits, or 
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This can be inverted to give 
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which can be solved for z to yield 

z = f 



sin(7rg/p) 



^l + /2_2/cos(7rg/p) 



(E24) 



which is the appropriate condition for these 'inner touching' closed orbits. The corresponding 
classical periods, Tj^°\ which then give the corresponding path lengths, are given by 
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where we have used the relation in Eqn. ()E24|) and the classical connection E = and 
recover the expected result. 

The spherical billiard can also be discussed in the same context, using th ese m ethods. 
The centrifugal barrier term, in a WKB expansion, is obtained from Eqn. ()28H1 by th e 
substitution l{l + 1) ^ l{l + 1) + 1/4 = (/ + 1/2)^, the so-called Langer modification (fl78h - 
(EH. 
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FIG. 1 Plot of the autocorrelation function, over the first 100 classical periods, for the 

model system described by Eqns. fl^i and fl^i and Case A in Table I. The value of at 
multiples of the classical period (squares) and at half a period out of phase (stars) is shown in the 
bottom (b), compared to the 'incoherent' value of |Aj„cP = Xln ~ l/(An2-y/7r) w 0.047 in 
Eqn. I|30|) . shown as the dotted horizontal lines (and indicated by arrows.) Locations of fractional 
revivals are indicated by vertical arrows. 
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FIG. 2 Plot of the autocorrelation function, over the interval {Q,Xrev) for the energy 

spectrum in Eqn. H26|) and the Gaussian distribution of energy eigenstates in Eqn. 1)241) with 
no = 400 and a = 1/800 and three different sets of parameters: (a) An = 3, and /? = 0, (b) 
An = 6, and /3 = 0, and (c) An = 6, and (3 = 2 x W~^. The horizontal lines indicate several 
fractional values of |Ap = 1, namely 1/2 (dashed), 1/3 (dotted), and 1/4 (dot-dash), at rational 
fractions of t/Trev given by 1/4, 1/3, and 1/8 respectively. 
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FIG. 3 Plot of the autocorrelation function, |74(t)p, near various full or fractional revivals, i.e. 
over the intervals {pTrev/l — Tci,pTrev/Q + Td) for various values oi p/q and for the parameters in 
Eqns. (|24() and (|26() and Table I. The values at the bottom are multiplied by 10 to illustrate the 
oscillations about the constant incoherent value |Aj„cP (shown as the dotted horizontal line) for 
times not near resolvable fractional revivals. For comparison, the constant value of \A\'^ = p/q = 
1/37 is indicated by the horizontal dashed line to compare a large q fractional revival versus a 
collapsed or incoherent state. 
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FIG. 4 Numerical calculations by Parker and Stroud (originally from Ref. (3)) for the intensity 
of the ionization probability for atoms excited by short laser pulses, showing numerical evidence 
for revival behavior, as explained in Ref. jsS) ■ The additional notations, identifying the locations 
of the fractional revivals, were added by Averbukh and Perelman. (Reprinted from Ref. (3)-) 
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FIG. 5 Recurrence spectra for the excitation of a Rydberg electron wave packet from Ref. 
showing fractional revivals of up to order T^evH ■ The top plot shows measurements taken using 
rubidium atoms, with no = 46.5 and roughly 6.5 states excited, while the inset shows data for 
TiQ = 53.5 and 9.9 states, showing the improved resolution as one increases An. The bottom 
plot shows the calculated time-dependent transition probabilities for the same system. We note 
that some of the identifications with various fractional revivals differ from the notation used by 



Averbukh and Perelman (;34) . (Reprinted from Ref. 



. fl.) 
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FIG. 6 Plot of the autocorrelation function, |74(t)p, over one classical period, for the 'pulsat- 
ing' Gaussian wave packets solutions in the harmonic oscillator, with xo,Po = 0, described by 
Eqn. H13U|) . for various values of r = 0^ / Note that the plots are invariant under r — > 1/r. For 
r = 1 the solution reduces to the ground state energy eigenstate, with trivial time dependence, and 
\A{t)\ = 1. 



109 




to 



5 10 15 < 




FIG. 7 Plots of the autocorrelation function, over one revival time, for zero- momentum 

Gaussian wave packets in the infinite well, for various values of the initial position, xq. Note the 
special shorter time-scale revivals for cases with additional symmetries, such as for (a) xq = L/2 
(Trev/8) and (b) xq = 2L/3 (r,.e^/3). The corresponding expansion coefficients from Eqn. H153|l 
are shown alongside. 
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FIG. 8 Same as Fig. [3 but for xq = L/2 wave packets with increasing momentum values. Note 
the appearance of the classical periodicity for po >> ^Po for the bottom two cases, where the \an\ 
versus n distribution approaches the expected Gaussian form in Eqn. ()24l) for no >> 1. 
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FIG. 9 Plot of the auto-correlation function (log(|yl(t)p) for the Gaussian wave packet in the 
ID infinite well for the first classical period, for increasing values of the number of eigenstates 
(N = 20,40,80) (with a„ given by Eqn. (|153() 1 used in the expansion (solid curves). The free- 
particle autocorrelation function in Eqn. p64j) as a function of t (dashed) and Td — t (dot-dash) 
are shown for comparison. The values plotted as diamonds are from the generic result given by 
Nauenberg in Eqn. ()B3|) . 
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FIG. 10 Plot of |j4(t)p (solid) and (dashed) over the first 15 classical periods, along with 

the dispersive part of the free-particle autocorrelation function, |j4ent,(t)P (dotted), in Eqn. ()165() . 



The value of the spreading times, to are shown for comparison. 
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FIG. 11 Plot of \A{nTci)\^ (square) and \A{nTci)\'^ (stars) over the first 150 classical periods, 
along with the dispersive part of the free-particle autocorrelation function, |Aeni;(i)P (solid), in 
Eqn. H165|l . The horizontal dotted line corresponds to the incoherent sum of uncorrelated energy 
eigenstates, given by Eqn. ()166() . 
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FIG. 12 Plot of |A(nTc/)| (top) and |74(nrci)| (bottom) over half a revival time (just over T^ev/'^-) 
A large number of possible fractional revivals are denoted by vertical dashed lines, while the value 
of \Ainc\ for an incoherent sum of eigenstates, as in Eqn. (|166() . is also shown for comparison. A 
number of the higher order fractional revivals are obvious. For example, the revivals at T^ev/'^ 
and Treii/6 are obvious in the 'out-of-phase' A{t) anti-correlation function, while that at T^eti/S is 
apparent in A{t); the fractional revival at Trev/^ is equally distinct in both plots, consistent with 
the form in Eqn. (|43|) . 
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FIG. 13 Plots of \A{nTcL)\'^ over one revival time (upper left) and over two classical 

periods, centered about Tj-ev (upper right). The autocorrelation and anti-correlation functions 
near t = Trev/^ (lower right) and the position-space probabilities at t = Tj-ev/^ (dashed) and 
Tj-evf^ + Tci/^O (dotted) are shown, illustrating the mirror revival, the wave packet being reformed 
on the opposite side of the well from the t = packet (solid) centered at xo/L = 2/3. 
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FIG. 14 Same as Fig. llHL except illustrating the fractional revival at T^e^/S. The horizontal dotted 
line corresponds to = 1/3, and xq/L = 1/2 was used. The same fractional revival is also 

visualized using the Wigner distribution in Fig. 1291 
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FIG. 15 Same as Fig. llHL except illustrating the fractional revival at The horizontal dotted 

line corresponds to = 1/2, and xq/L = 2/3 was used. The same fractional revival is also 

visualized using the Wigner distribution in Fig. 1301 
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FIG. 16 Same as Fig. WM except illustrating the behavior of both the probability density (bottom 
left) and A{t),A{t) (bottom right) during the collapsed phase, where P{x,t) = \ip{x,t)\'^ is more 
consistent with the uniform or 'flat' value of P{x) = 1/L (which is indicated by the two horizontal 
arrows.) The values of |j4(t)p, |74(t)p can be compared to the incoherent value in Eqn. (|166|) shown 
as the horizontal dotted line. 
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FIG. 17 Plots of the real (bottom) and imaginary (right) parts of A(t) versus time (over half 
a revival time) as well as a parametric plot (Argand diagram) of the same data. In this case, 
A{t) is evaluated at integral multiples of Tc/, as in Figs. El The solid circle corresponds to 
|^(t)P = 1, while the dashed circle corresponds to |^(t)p = 1/2; the location of a quarter-revival 
{^revl^ is denoted by the diamond, illustrating how at fractional revivals the autocorrelation 
function (sampled at these intervals) is tangent to the circle. 
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FIG. 18 Same as Fig. llTL but for the time period T^i on either side of the T^ev/^ fractional revival to 
be compared to the general wavefunction in Eqn. (|43() and the correlations predicted by Eqn. ()167() . 
The dashed circle corresponds to \A{t)\ = l/\/2. The dotted lines correspond to the directions 

g+i7r/4 and e-''"/'^. 
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FIG. 19 Same as Fig. I17L but for the time period T^i on either side of the T^ev/'^ fractional revival 
and for all times, not just multiples of T^. Note the highly correlated phase relationship consistent 
with the general results of Eqn. (|58|) and in Eqn. H169j) . In this case, the dashed circle corresponds 
to \A{t)\ = 1/^/3. 
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FIG. 20 Same as Fig. I17L but for the time period T^i on either side of a typical time (T*) during 
the collapsed phase, not near any resolvable fractional revival. In this case, the dotted circle 
corresponds to \AinS^ = X] jonl^ = 1/An2y^, typical of the incoherent sum of many eigenstates, 
with little or no phase relationship obvious. 
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FIG. 21 The classical motion of a particle in an infinite well potential. 
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FIG. 22 The expectation value, {x)t (top), and uncertainty, Axt (bottom), for the infinite Gaussian 
wave packet over the first ten classical periods, measured in units of L. The dotted line which forms 



the envelope for the Axt curve is the free-particle spread given by Eqn. (|170a() . (Reprinted from 
Ref. (94).) 
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FIG. 23 Same as for Fig.|221 but for {p)t (top), and Apt (bottom). (Reprinted from Ref. (jM).) 
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FIG. 24 The position-spread, Axt versus t, for various initial width Gaussian states, and the 
approach to the collapsed state, characterized by Axd = L/vT2, evaluated at integral values of 
the classical period. (Reprinted from Ref. (3)-) 
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FIG. 25 The same as Fig. |^ but for {x)t (top), {p)t (middle), and Apt (bottom). (Reprinted from 
Ref. (94).) 
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FIG. 26 The expectation value of position {x)t evaluated at t = (n + l/S)Tci (solid) and t = 
(n + 5/8)Tci (dashed) over an entire revival time. The mirror revival, where the wave packet is 
reversed, at t = T^ev/'^ is apparent. (Reprinted from Ref. l94f).) 



129 



Ax 



free 



1 1 1 1 1 1 1 1 1 1 1 1 1 r 



X 
< 



r 0.2 



X 
< 



r 0.29 




0.27 



_l I L_ 



120 



130 



140 



150 



160 



t/T 



cl 



FIG. 27 Same as for Fig. 1261 but for Axt over one revival time. The return of Axt to its original 
(small) value at T^ev (revival time) and Trevf^ (mirror revival) are familiar, while the small values 
at the Trev/^ and STj-ev/^ fractional revivals are special features of an initial state centered at 
xq = L/2. At those points, the wave packet consists of two 'mini' packets, superimposed and hence 
of small width, but with opposite momenta. (Reprinted from Ref. Q.) 
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FIG. 28 Same as for Fig. I26L but for {p)t and Apt. Note that there are no obvious deviations from 
the collapsed value of Apd = +po at T^ev/^ (200T^/) or 3Tre„/4 (600Tc/), even though Ax^ returns 
to its original (small) size. (Reprinted from Ref. (94).) 
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FIG. 29 Phase-space structure of the T^ev/S fractional revival using the Wigner distribution. The 
parameters of Eqns. (|24() and (|26|) are used, with xq/L = 1/2, but with no = 40 used to make the 
oscillatory 'cross-terms' more obvious. 
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FIG. 30 Same as Fig. I29( but for the fractional revival at t = Trev/4:, where xq/L = 2/3 is used. 
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FIG. 31 Plot of the classical trajectories (or world-lines), shown in the {x,t) plane, of three point 
particles with slightly different speeds, illust ratin g wave packet spreading in the infinite well with 



walls at X = 0, -L, as first discussed by Born (|lOl|), (|l02l ). 
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FIG. 32 Plot of the position-space probability density, P{x, t) = \il!{x, t)\^ , versus x, t in the infinite 
well, for a Gaussian wave packet, over the first classical period. The parameters in Eqn. ()159() . 
with no = 400, are used, so the infinite well is over the interval x G (0, 1). 
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FIG. 33 Same as Fig. 021 but over one classical period, starting at Trei,/3. 
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FIG. 35 Plots of typical worldlines in the (x, t) plane corresponding to the classical contribution 
to the probability density in Eqn. (|195() where n + m >> 1 (a), and two cases from the quantum 
carpet terms in Eqn. ()196(1 corresponding to |n — m| = 1 (b) and |n — m-| =2 (c). 
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FIG. 36 Gaussian wave packet solutions for the quantum bouncer in position-space {\'ip{z,t)\'^ 
versus z, left) and momentum-space t)p versus p, right) for various times during the first 

classical period. The solid curves represent the time-dependent expectation values of position 
{{z)t, left) and momentum {{p)t, right) for these solutions. The similar dashed curves are the 
classical trajectories, z{t) (left) and p{t) (right), superimposed. The wave packet parameters in 
Eqns. (|217j) and (|218j) are used. For the momentum-space figure, the vertical dotted lines represent 



the values p = and the classical extremal values of momentum, ±pm = ±y/2mFzo. (Reprinted 
from Ref. il5lh.1 
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FIG. 37 The top plot (a) shows the position-space probability density, for the initial 

Gaussian wave packet (solid), at at times near the first (dashed) and second (dotted) revival times. 
The bottom plot (b) shows \^p{z,T*)\'^ at a typical time, T*, not close to any fractional revival, 
during the collapsed phase (solid), while the dashed curve is the classical probability density given 
by Eqn. (|22H1 (but note the change in scale.) 
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FIG. 38 Plots of the auto-correlation function, |74(t)p versus for the 2-D square billiard. The 
plots are over a time period equal to ten classical 'back-and-forth' periods, lOr, where r = 2L/vq. 
Plots for different values of the initial angle, tan(0) = poy/Pox are shown. The stars indicate 

the positions of classical closed orbits (and recurrences) as given by Eqn. (|238|) . (Reprinted from 

I I 

Ref. (97).) 
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FIG. 39 Plot of the autocorrelation function, versus t, in units of Tq = 2^1^ /Htt, The 

numerical values of Eqn. (|293j) are used, along with yo = and pox = Poy = 0. The results for 
versus t, as one varies the xq of the initial wave packet away from the center of the circular 
billiard, are shown on horizontal axes. (Reprinted from Ref. (177).) 
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FIG. 40 S ame as Fig. [HHl but with xq = yo = and pox = 0, as one increases poy. (Reprinted from 
Ref. Iml).) 
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FIG. 41 Experimental (top) and theoretical (bottom) photoionization signal from Rydberg atoms 
as a function of the time delay after the initial probe pulse (a) compared with theoretical predictions 
(b). A fractional revival of order p/q = 1/4, with half the classical periodicity, is apparent near 
680 ps. (Reprinted from Ref. \2>li \.) 
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FIG. 42 Log spectrogram contour plot of data from Br2 showing the spectral content of the 
observed signal as a function of time delay, defined by Eqn. (I^llj) . Evidence for full and fractional 



revivals at multiples of / = 95 cm ^ are clear. (Reprinted from Ref. {2^, courtesy of D. Villeneuve.) 
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FIG. 43 Plots of Pe{t) versus scaled time (r = Xt/ir) from the solution in Eqn. (|318)) over (a) the 

first few Rabi cycles and (b) over a number of revivals times, showing the first (c) and second (d) 

revivals in detail. Values of n = 36 and A = 0.01 are used. The dotted curve corresponds to the 

initial Gaussian de-phasing envelope (exp(— (At)^/2), while the dashed curves correspond to the 

long-term suppression given by Eqn. (|321|) . 
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FIG. 44 Plots of P{P;t) = P{Re{P), Im{p);t) versus {Re{[3),Im{(3)) from Eqn. Values 
for (a) t = 0, (b) O.lTret,, (c) T,.et,/2, and (d) T^e^/S are shown. An initial coherent state with 
a = (3, 0) is used for illustration. 
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FIG. 45 An example of two possible closed orbits in the annular billiard, with {p,q) = (5,2), 
corresponding to the path lengths in Eqns. HE1|) (top) and Eqn. HE17|) (bottom), illustrating as 
well the critical value of fmax = cos(7rg/p), beyond which neither orbit can be supported in this 



geometry. 
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